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ABSTRACT 


Design studies are performed and problem areas 
analyzed with respect to the redesign of the MIT Reactor 
in order to accommodate a compact core configuration, 


Mechanical design studies are performed for a D,0 
moderated, cooled, and reflected compact core design. 
Designs are proposed for the core tank and reactor top 
Shielding configuration, the core support mechanism, the 
control rod drive mechanisms, and the refueling mechanism, 


tne mecessary composition and thickness of reactor top 
Shielding are determined and shielding calculations per- 
formed for both a D,0 cooled core and an HO cooled core, 
Shielding above the D0 core consists of an 8 ft. column of 
D,0 plus a shielding enclosure consisting of 6 in, of mason- 
ife, iW emineeOne poral, and | fu. of lead, through which 4 
maximum dose rate of 0.2 mrem/hr is calculated. Shielding 
above the Ho0 cooled core consists of 11 ft. of Ho0 topped 
Dy a 10 in, lead cover, through which a maximum dose rate 
of 0.9 mrem/hr is calculated. 


ticmiemMimilOvuanammesy tLranster problems associated 
with the design of a compact core configuration are analy- 
Zea and it iS determined that the currently installed MITR 
process system is not capable of providing adequate heat 
removal from the compact core, if smooth flat plate fuel 
elements are utilized. An experimental investigation is 
undertaken to determine the effectiveness of using fins as 
a means of augmenting fuel plate heat transfer. Heat 
transfer from a smooth rectangular channel 2.5" x 0.090" x 
25-1/fe" is first investigated and experimental results 
Coniirm the validity of ug ing the modified Colburn corre- 
lation, Nu > Oe Oa Re, Eo, to predict smooth 
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channel heat transfer values, Two finned channel test 
Sections aré inveStigated, one with fins cut into the 
channel transverse to the coolant flow direction, and a 
second with fins milled parallel to the coolant flow | 
Sirection, in each case the fins are milled to be 20 |: 
mils high, 20 mils wide, and spread 40 mils apart. The 
transverse fins are found to augment heat transfer by a 
factor of 2.2, while a factor increase of 1.6 is obtained 
from the longitudinal fins. However, due to safety, 
pressure drop and cost considerations, longitudinally 
finned fuel plates are recommended for use in the High 
Flux MIT Reactor. 


eupervisor: Dr. Theos J. Thompson 


Title; Professor of Nuclear Engineering 
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INTRODUCTION 

Since its completion in the Spring of 1958, the MT 
Reactor has served as a vital research tool for the MIT 
and surrounding scientific community. At the time the 
MIT Reactor was designed, little was known of compact, 
heavy water moderated cores. Consequently, the MIT Reactor 
was designed with fuel element spacing sufficient to pro- 
vide for a high percentage of neutron moderation to take 
place within the core area. Results published in 1964 
upon the completion of the Brookhaven High Flux Beam 
Research Reactor provided the experimental information 
needed on compact, undermoderated core configurations in 
order to accomplish an exaggerated thermal flux peak in 
the reflector. By incorporating a similar compact core 
design into the MIT Reactor, it is reasonable to assume 
that the thermal flux coming out of the beam ports could 
be increased by at least a factor of three, and possibly 
as much as a factor of four or five. Such an increase 
would markedly increase the value of the MIT Reactor to 
research, In addition, if the beam ports were to be 
located below the base of the core, the core gamma ray 
and fast neutron background would be markedly reduced in 
intensity due to the fact that any core gammas and fast 


neutrons which leak out must have been scattered at 





ED 
least once previous to eScaping. Such a design would 
markedly increase the quality of the thermal neutron 
beam emitted from the beam ports. 

it is pihiarily for the above reasons that a MITR 
redesign group was formed under the direction of Dr. 

T.J. Thompson, The design group is mermoued of both 
graduate students and MITR staff personnel. Besides Dr. 
Thompson and the author of this thesis, the members of 
the design group who have been active throughout the re- 
- design process are graduate students: Andrews Addae 
working on core physics, Dan Kennedy working on beam port 
optimization, and Bob Sanders working on a cold neutron 
facility for the thermal column. Mr. Ed Barnett, working 
on all engineering phases of the redesign, represents the 
METR staff. 

Although never stated formally in writing, the design 
criteria outlined orally by Dr. Thompson at the outset of 
the design were the following: 

1) The thermal neutron flux at the beam ports must 
be increased by a factor of 3-5. 

2) The redesigned reactor must be as safe, or 
Safer, than the existing meaccor: 

3) Simplicity of design must be a prerequisite for 


the design of any component. 
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4) There must be no buried mechanisms. 

5) There must always be a means for retreat, To 
quote Dr. Thompson, "We have a going reactor, we don't 
want to replace it with a monument," 

6) The redesigned reactor must be reliable. 

7) The cost of the redesign project must be within 
$500 ,000. 

The redesign project has actually been divided into 
two distinct design phases. The first phase began in 
earnest in the Fall of 1967 and was centered around the 
design of a D,0 moderated, cooled and reflected reactor. 
Phase I concluded in the Fall of 1968 and was succeeded 
by an intermediate phase in which a mixture of H50 and 
D0 was considered as the moderator and coolant for the 
eore, with pure D509 ee Ne test Weewere shel (Betersjorenes (onl 
1968 phase II began. Phase II is centered around the 
design of a light water cooled and moderated core surrounded 
by a heavy water reflector. 

At the outset of phase I, it was apparent that there 
would be many challenges in the area of mechanical design. 
ihe incorporation of a compact core configuration dictated 
that new methods be devised to support not only the fuel 
elements, but the control elements as well. A change in 
core tank design Seemed to be required. This author's 
first responsibility as a member of the redeSign team was, 


then, to attempt to find solutions to some of the many 





Ji 


design problems confronting the redesign project. The 
work done toward this objective is summarized in Chapter 
e Of this thesis. 

As the redesign project continued and a necked cool- 
ant tank was adopted, concern for the requirements of 
reactor top shielding emerged. This author's second area 
of responsibility became, therefore, that of designing the 
shielding to be required for the reactor top, and that of 
performing shielding calculations to justify the proposed 
Shield composition. The results of work done in this area 
are presented in Chapter 3 of this thesis. 

As stated at the beginning of this introduction, one 
of the primary purposes for redesigning the MITR was to 
accommodate a compact core configuration in order to obtain 
a large thermal flux peak in the reflector. However, 
almost by definition the achievement of a pronounced ther- 
mal flux peak in the reflector will be accompanied by a 
Sharply increasing thermal flux gradient across the outer- 
most fuel plates. The outermost fuel plates will conse- 
quently generate a much higher power density than the 
inner fuel plates in the core. From the results presented 
in reference (1), it is reasonable to suspect that due to 
the sharply rising thermal flux gradient in the area of 
the outermost fuel plates, the reactivity coefficient due 


to coolant voiding in the channels adjacent to these fuel 
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plates may be positive. It, therefore, becomes imperative 
that pouiene not be allowed to occur in the area of the 
outermost fuel plates, which is equivalent to Saying that 
porting must not occur anywhere in the core Since the 
outermost fuel plates generate the most power. The maximum 
rate at which heat can be removed from the outermost fuel 
plates without allowing nucleate boiling to occur thereby 
sets an upper limit to the power density allowable in the 
outermost fuel plates, which in turn influence the maximum 
thermal flux peak which can be obtained in the reflector. 
It is this critical area of core heat removal which became 
this author's third and final area of responsibility. The 
results of work done in this area are presented in Chapters 


4 and 5 of this thesis. 
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Chapter 1 


DESCRIPTION OF EXISTING MITR COMPONENTS 

ie General 

Before proceeding with a description of the design 
and calculational work which was performed as a part of 
the High Flux MITR Project, it seems appropriate to first 
present a brief description of those existing MITR compo- 
nents and systems which are the subject of redesign work 
presented in this thesis. The following information was 
taken from the operating manual for the MIT Reactor (2). 
The operating manual should be consulted if more detailed 
descriptions are desired, 
ie Reactor Coolant Tank 

The reactor tank is 48 inches in inside diameter and 
has a wall thickness of 1/4 inch, except at the bottom 
where it is 3/8 inch thick. The overall height from the 
dished bottom shell to the top flange is 7 feet 5-3/8 
inches, The reactor tank is made of special high purity 
aluminum and has a design pressure of 40 psig. During 
normal operation, ROWSE | the pressure in the tank above 
che coolant is only 1 inch son HO above atmospheric pres- 
Sure, The tank is shown in the vertical section through 
the reactor of Figure 1.2.1, and is also shown in Figure 


ee. o. 
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VERTICAL SECTION OF REACTOR 


FIGURE 1.2.1 
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ANNULAR RING 

ce . PPORT BLOCKS 

Lb 
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REACTOR CORE TANK SUPPORT 
FIGURE 1.2.2 





ee 

Heavy water flows into the reactor tank through an 
8 inch diameter aluminum pipe, through the holes in the 
Side of the inlet pipe cap and out into the entrance 
plenum. It then flows up through the fuel elements and 
out into the tank itself. The D,0 from all the fuel ele- 
ments flows out of the tank through a transition section 
which extends a short distance above the plenum head, and 
then through the 8 inch aluminum outlet pipe. The normal 
total coolant flow rate is approximately 2000 gallons per 
minute, 

The reactor tank is attached to the lower shielding 
plug by 24 one inch stainless steel socket head cap screws. 
The cap screwS compress a hollow stainless steel O ring 
between the tank and the shield plug to provide a helium 
seal designed for 80 psig. ConSequently, if a high pres- 
Sure is developed in the tank, the tank will rupture 
downward before the shield plug gives way. 

1.3 Shielding Above the Core 

1.3.1 Lower Annular Ring 

The lower annular ring (shown in Figures 1.2.1 and 
1.2.2) is a concrete-filled steel structure which is 56 
in. in inside diameter, approximately 8 ft. in outside 
diameter, and 4 ft. high. Its lower surfaces are covered 
with 1/4 in. of boral. A 4 in. thick, water-cooled lead 
thermal shield lines the inner surface over the lower 18 


in, There are conduits for thermocouple leads and helium 
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gas pipes, and 6 sleeves are provided for thimbles which 
Slant into the graphite reflector. | 

1.3.2 Lower Shield Plug 

The lower shield plug is a large cylindrical struc- 
ture, the bottom of which is lined with a 1/4 in. sheet 
of boral. Above the boral is a 1-1/8 in. steel plate, 4 
tot stead provided with cooling coils, 29 in. of concrete, 
a steel cover plate, a 6 in. space filled with paraffin and 
steel shot, and a 1-1/8 in. stainless steel top plate. The 
plug is penetrated by numerous sleeves for the concrete- 
filled steel shielding plugs which are part of the fuel ele- 
ment and control rod assemblies. 

1.3.3 Upper Annular Ring 

The upper annular ring, together with the rotary plug, 
complete the biological shielding above the core tank. The 
ieee LOOU thick annulus iS made of dense concrete encased 
in steel, Plugged holes are provided for access to the 
vertical thimbals which protrude through the lower annular 
ring and into the graphite region. The upper annulus is 
spaced 9-1/4 in. from the lower annulus by several support- 
ing blocks bolted to the lower annular ring flanges. 

1.3.4 Rotary Plug 

The heavy concrete rotary plug fits into the center 
of the upper annular ring and rests on a ball-bearing 
thrust ring. It is rotated by a drive motor and friction 


belt. When the plug is latched in its operating position, 
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the control rod racks can move upward into cylindrical 
slots in ee bottom of the plug. To remove a fuel ele- 
ment, the control rods must be fully inserted and the 
plug unlatched and rotated to the appropriate ander oom 
position. The shielding plug directly over the fuel ele- 
ment position can then be removed. The entire rotary 
plug must be removed for access to the control rod drives. 
1.4 Fuel Element Assembly 

A fuel element assembly is shown in Figure 1.4.1 and 
consists of a fuel Section supported by a lower adapter, 
an upper adapter and a shield plug. The 30 positions into 
which fuel assemblies can be inserted are indicated by 
numbers 1 through 30 in Figure 1.4.2. There is a central 
element position surrounded by an inner circle of 6 posi- 
tions, all spaced 6.375 in. center-to-center. These in 
turn are surrounded by an outer circle of 12 positions 
spaced 6.858 in. center-to-center. The 11 additional 
positions are located around the outside edge of the reac- 
tor tank, 

The fuel section consists of 18 curved plates assem- 
bled between two grooved side plates which are 0.187 in. 
thick and 3,168 in. wide. Two nonfuel-bearing outer 
plates are 0.060 in. thick and are curved to a 5-1/2 in. 
meaius. They are spaced 0.117 in. apart to form 17 cool- 


ing water passages. 
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The fuel assembly shielding plug fits into a stepped 
Steel sleeve embedded in the lower Shield plug. The 
weight of the assembly rests on the step in the lower 
Shield plug. This provides a convenient means of removing 
the fuel assembly through the tank top, and when in place 
the shielding properties are equivalent to that of the 
lower shield plug. The shield plug is securely latched 
into place by two Small lugs which slide outward and 
engage a groove in the sleeve. The lugs are moved by turn- 
ing a square projection on a cam which bears against the 
lug. 

1.5 Control Rod Mechanisms 

The MITR has one regulating rod and six shim safey 
rods arranged in the core as shown by the solid black 
circles in Figure 1.4.2. The shim safety rod assembly is 
Shown in Figure 1.5.1 and is described below. 

The absorber section Te hollow, double-walled, 
aluminum cylinder, 28 in. long and 2-1/4 in. in outside 
diameter, in which is sandwiched a 0,040 in. thick cadmium 
ipe 2.15 an. in outside diameter and 26-1/8 in long. The 
absorber cylinder is connected to the magnet armature (a 
nickel-plated soft on disk) and dashpot by means of a 
Spider-like section, thus permitting D50 to flow freely 
through the top of the cylinder as it drops. During a 
reactor scram as the absorber nears its fully inserted 


position, it is decelerated by a cushion of heavy water 
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which leaks off between a washer and the dashpot case 
with the assistance of coil Springs in the dashpot. A 
1/4 in. thick disk of cadmium on the bottom of the dashpot 
Shields the iron armature from neutrons.: 

The electromagnet consists of a silicon-insulated 
copper wire, wound onto a transite bobbin. It is fitted 
into an Armco magnet-iron housing which is enclosed by two 
Stainless steel end plates. The magnet wires are connected 
to a threaded connector, which, when coupled with the racks 
forms a hei ietically sealed unit. Magnet wires are led 
through a tube in the rack rod, then down the outside to 
a take-up real mounted at the top of the casSing assembly. 
The rack is guided by the casing, and the escape of helium 
or D590 around the shaft is prevented by Styrene chevron 
packing. 

The rods are driven by a rack and pinion arrangement. 
Pfesorive shait Trem the pinion gear extends to the motor 
and gear box assembly at the utility shelf (see Figure 
1.5.2). The six drive shafts are arranged radially but 
differ in length and arrangement in order to clear shutters 
maawexperimencal facilities, All drive shafts are enclosed 
in sleeves from the utility Shelf to the between-plug space, 


4} 


and are sealed to limit the amount of A entering the air 


stream, 
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Chapter 2 

PRELIMINARY MECHANICAL DESIGN STUDIES DONE FOR PHASE I 
2, | “General 

In the early stages of the MIT reactor design project 
it was this author's primary responsibility to attempt to 
find solutions to some of the many design problems criti- 
cal tO “ene Comsumuct tem of a High Flux MI@R. In so doing 
many design ideas were proposed, few of which gained 
acceptance by the entire design group. It is the purpose 
of this chapter to relate some of the design proposals 
which can be primarily attributed to this author, though 
at times when working with a design group it is admittedly 
difficult to distinguish which ideas are completely one's 
own from those which may have had their roots in, for 
example, a comment made at a design group meeting. It is 
felt that although a particular design scheme may not have 
appeared best for the MITR redesign, it may at some future 
date prove useful to another reactor design project, and 
therefore is worthy of description in this chapter. In 
addition, the design proposals serve as a kind of chronology 
of mechanical design progress. 
2.2 Core Tank and Reactor Top Shielding Configuration 

The basic necked core tank design as finally adopted 
evolved from the compact core concept and short tank pro- 


posed in the original redesign sketches of Prof. T.J. Thompson 
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and first formally reported in August 1966 in a 22.90 
project report by D.M. Goebel (3). The original purpose 
in designing a necked core tank was to provide for adequene 
wuclaing Of the core during refueling operations. (Prior 
to the introduction of an enclosed refueling structure 
Beop the reactor, 1t was proposed that there be a single 
refueling plug which would be removed during refueling 
operations. Without approximately 11.5 feet of D50 separa- 
bing, Che cere Grom the reactor top, the radiation level 
would be prohibitive to workers in the area of the open 
refueling plug.) In order to accompolish the objective of 
providing adequate shielding during refueling operations, 
put at the same time to minimize the poomtered D50 inventory, 
the concept of two rotating shield plugs was conceived and 
woipanally prepesea=by this author. Tmis design concept, 
outlined in Figure 2.2.1 would utilize a 1-1/2 ft. thick 
lower Shield plug fabricated of lead which would rotate 
PmounG a pracm oi aluminum rollers, Releasable latches 
would prevent the shield plug from beint "shot" upward by 
ea sudden increase in reactor pressure. Inserted in the 
lower shield plug would be three 6 in. diameter removable 
lead plugs, so situated that by rotating the lower shield 
plug one of the smaller plugs would be caused to line up 
atop any desired fuel element. <A tool inserted from the 


reactor top through the wepereshweld plug could be used to 
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unlatch a 6 in. diameter lead plug, lift it out of its 
elie, and Set 1G aside. The upper rotating shield plug 
would consist of 3 ft. of dense concrete, and would rotate 
on stainless steel rollers and contain three 6 in, diameter 
removable plugs similar to those of the lower shield plug. 
An inflatable boot would be located around the upper shield 
which could be deflated to allow the shield plug to rotate 
and inflated to act as a helium seal. (An inflatable boot 
could be installed around the lower shield plug if the 
upper shield plug should be entirely removed in order to 
perform maintenance on the control rod drive mechanisms, 
for example.) Both the upper and lower Shield plugs would 
Be cooled by duplicate D0 cooling coils cast into the 
shield»plugs. When rotated to an indexed position for 
reactor operation, the coolant lines in the shield plugs 
would be caused to line up with their respective supply 
Himes in the annular rings. 

The refueling operation would proceed somewhat as 
follows: One of three 6 in. diameter inserts in the lower 
Shield plug would be unlatched and raised by a tool inser- 
ted from above. The lower shield plug would then be rota- 
ted to an indexed position aligning the opening in the 
lower shield plug over the fuel element to be removed, 

The 6 in. diameter plug could then be rested on the lower 
shield plug for convenience. Next the appropriate 6 in. 


diameter insert from the upper shield plug would be removed. 
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At this time, the insert would not be aligned with the 
opening in the lower shield plug. The refueling flask 
| 


vould be Drought into position ome the reactor top, the 
upper shield plug rotated to an indexed position so seid 
align its opening above the fuel element to be removed, 
and the designated fuel element raised up into the refuel- 
ing flask. Once in the refueling flask, the fuel element 
could either be flipped end for end and reinserted into 
the reactor, or set aside while a new fuel element is 
placed into the core, The above procedure would be 
accomplished in reverse to complete the refueling operation. 
By utilizing the concept of a double, rotating shield 
plug, the reactor D509 tank would remain shaped similar to 
the present D,0 tank (see Figure 1.4). Two important 
differences would be apparent, however. Although 4 ft. in 
diameter, the D0 tank would be shortened from 7-1/2 ft. 
to 6-1/2 ft. in height. The volume of D0 in the system 
would thereby be decreased by 14%. Contained within the 
4 ft. diameter outer tank would be an inner tank of diame- 
ter just sufficient to contain the core. This inner tank 
would extend the full length of the outer D50 tank and 
would be constructed of sufficient thickness of aluminum to 
be able to contain a sudden D50 vaporization resulting in 
4O psi overpressure in the core. AS a consequence, the 4 


ft. diameter D50 tank need only be of sufficient strength 
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to support the coolant contained within it. There are two 
major reasons for utilizing an inner tank to separate the 
core region from the reflector region. First is the 
possibility of maintaining the reflector D350 at a lower 
bemperature than the coolant D0 ; and thus make the reflec- 
cor D509 a better moderator. Second is that by separating 
bae reflector from the core, the entire reflector region 
could be rapidly dumped of its D530; thus creating a very 
imporcan. mechanism for reactor safety. 

Both the upper and lower annular rings would remain 
the same as in the present MITR, and the core tank likewise 
would be supported as it is now. Such mechanical similari- 
ties with the existing reactor would tend to reduce the 
C0st of convertimg to the high flux MITR core configuration. 

iheweOneroimrod mechanisms for the above described 
reactor tank configuration would extend through the lower 
annular ring in a manner similar to the present MITR 
design. The absorber section would protrude into the re- 
flector region adjacent to, but outside of the inner ‘tank 
wall. The control rod drive motors and connecting shaft 
work could remain and be used in their current positions. 
Maintenance could be performed on the control rod mechanisms 
by removing the upper shield plug, thus allowing a man to 


enter the area between the upper and lower annular rings. 
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From the above description it would appear as if the 
above described core tank design with a double, rotating 
shield plug mechanism has no drawbacks and should have been 
adopted. Where the above design fails is that it does not 
sufficiently conform to two of the design criteria mentioned 
in the thesis introduction which have played a major role in 
the selection of final mechanical designs for the High Flux 
MITR. These criteria are that the design selected should be 
Eee ult imace in Simplicity, and Should have no hidden 
‘mechanisms. The mechanics involved in the concept of a 
double rotating shield plug with an inflatable helium seal 
are not what could be described as simple, and the mechanisms 
associated with the lower rotating shield plug could be con- 
eeaered as hidden, 

2.3 Core Support Mechanism 

2.3.1 General 

Cieinajormeemeern Coevhe redesign of the MITR core area 
was the method by which the core is to be supported. It is 
possible that the core receive its primary Support either 
from below or by being hung from above. The core support 
mechanism selected for use in the phase J High Flux MITR 
(designed by Mr. E. nee supports the core from below. 
the overwhelming reason for the choice of this mechanism to 
Support the core from below was safety. In the deSign of a 


reactor in which the sudden removal of a fuel element from 
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foecweOce DPOsit1on would cause a reactivity addition on the 
order of 3% in AK/K, as is the case for the High Flux MITR, 
Lt 1s imperative that each fuel element be supported in an 
absolutely foolproof manner (4). 

There is, however, one very compelling reason for 
wanting to Support the core from above. The aluminum struc- 
iure necessary to support the core from below results ina 
loss on the order of 20% in thermal neutron flux in the 
area below the core (5). It is this area just below the 
core at which the beam ports "look". In addition, having 
the aluminum structure below the core results in an 
increased gamma flux seen by the beam ports due to the ther- 


mal neutron activation of hoot. 


2.3,e¢ @Grid Support 

The first design proposed as a means of supporting the 
gore from above Hs Shown in Figure 2.3.2.1. In this design 
the core is hung from a grid-like structure, the grid being 
Supported either directly by the lower annular ring in the 
case of the double shield plug design, or from the base of 
the necked region of the core tank in the case of the necked 
core tank design. An adaptor which is capable of being 
attached to either end of the fuel element would be connec- 
ted by four aluminum rods to the corners of a rectangular 
plate. The plate would rest on top of the grid support bars 
when the fuel element is inserted into the core. A hinged 


cover would be lowered over the grid structure and latched 





39 


FIGURE Zoe | 
GRID SUPPORT 


ONE ELEMENT OF 
THE GRID SUPPORT 





ADAPTER FOR 
CONNECTING TO 
FUEL ELEMENTS 





ELEMENT 





4O 


in order to hold the fuel element Support plates down. 
There is both a major deSign concept advantage and a major 
disadvantage in this method of core support. The advantage 
menounas With such a megned both the core and the control 
elements are supported from the same structure, the lower 
annular ring. Therefore, should a major catastrophe (such 
aS an airplane crashing into the reactor or an earthquake ) 
eccur in which the reactor structure was disturbed, having 
me wCOre and cOmMmre. elements supported from the same 
Structure would markedly increase the probability of their 
Suaying, together; me., one sheuid not fall without the 
other. The disadvantage to the above method of core sup- 
port is that there must be a space between each fuel element 
of width equal to the uh Ormeaesrrcdmompport bar, The grid 
support bars must be approximately 1/4 in. wide in order to 
provide a sufficient area on which to rest the fuel support 
plates. This spacing would increase the core diameter by 
approximately 1-1/2 in., thereby decreasing the core power 
density, thereby decreasing the thermal flux peaking in the 
fet lector, 

2.3.3 Sliding Bar Support 

in order to Support the core from above, but without 
the spacing between fuel elements required by the method 
of the preceeding section, a design utilizing sliding 


Support bars was developed (see Figure 2.3.3.1). (A 
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design of this general type was proposed earlier by Dr. 
Thompson and Mr. Pee Assuming for the sake of 
{llustration, a 24 element core configuration, a System of 
14 east-west and 14 north-south tracks would be constructed 
around the core. The tracks would be attached at one end | 
to the outer tank wall and supported at the inner end by 
the inner tank wall plus a brace hung from the necked 
region of the reactor tank. On these tracks would slide 
the core support bars, seven bars for each quadrant of the 
core. Special frame type structures with keyhole indenta- 
tions would dock either end of the fuel elements and accept 
the sliding core support bars (figure 2.3.3.2). In order 
to refuel the support bars would be withdrawn in such a 
manner as to allow only one fuel element to be released at 
a time. An interlock mechanism would prevent the with- 
drawal of Support bars in any manner which would cause 
more than one fuel element to be released. The fuel ele- 
Ment Support structures being frame-type Structures with 
Openings on all Sides, the D50 coolant would flow up through 
the core and out the sides of the Support structure. 
Although not a necessity, the top of the support structure 
could be made solid in order to help prevent the D,0 cool- 


ant (with its © 


activity) from being circulated into the 
upper portion of the core tank after leaving the core. As 


Snown in Figure 2.3.3.3, the sliding support bar design 
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“ould require that the control elements be angular shaped 
Penecontorm co the outline of the core, 

The major shortcoming of the sliding bar support 
mechanism is the necessity to design an interlock which 
would be absolutely foolproof in preventing the accidental 
release of a fuel element. It can also be argued that 
this design certainly does not meet the "simplicity" design 
requirement. 

2.3.4 Safety Supports 

From a safety point of view, both the designs proposed 
in Sections 2.3.2 and 2.3.3 would be enhanced if some backup 
means of support were provided below each fuel element. For 
this purpose, a 1/4 in. aluminum rod set into the coolant 
entrance plenum would suffice. On top of the 1/4 in. 
aluminum rod would be a structure resembling a crow's foot 
which would act both to support the fuel element and to 
turbulate the coolant as it approaches the entrance to the 
fuel element. The support bar could be removed and replaced 
my a bar of different length should ieee decmmedero Change 
the height of any element in the core. The aluminum support 
bars serve only one function, that of acting as a backup to 
ensure that no fuel element is able to fall through the 
core. 

2.4 Control Rod Drive Mechanisms 
2.4.1 General 
Onemorecme most difficult design problems encountered 


in the redesign of the MIT Reactor was the location of the 
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Seerrol 1ou Crive mecnanisms. oince it was desired to re- 
tain the present concept of a magnetic coupling between 
the control absorber and withdrawal shaft, it became A 
design premis that the control rod drive mechanism be ee: 
Signed in such a manner as to allow the magnet to be loca- 
ted as far from the core aS posSible, and if possible out- 
Side the D509 tank. By doing so, it was hoped that magnet 
failures due to radiation damage and moisture corrosion 
could be elimineted, or at least reduced markedly in compari- 
son with the rate of failures in the present MITR. The 
magnetic coupling between the control rod drive motor and 
the control absorber was enviSioned to be either the general 
type used in the present MITR or a magnetic clutch arrange- 
ment. Some of the proposed designs and locations for the 
control rod drive mechanisms are discussed in the paragraphs 
which follow: 

2.4.2 Absorber and Drive Mechanisms Outside Inner Core 

Tank 

The first control rod drive mechanism considered was 
merely a modification of the present mechanisms (see Figures 
ime. and 1.5.1). In this design (shown in Figure 2.2.1), 
the absorber would be located in the reflector region just 
Beeside the inner core tank. The absorber would be connec- 
ted by a shaft to the electromagnet coupling in such a 


manner that the magnet is always located inside the concrete 
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shielding of the eee ring and therefore alSo above the 
level of the reflector D50. By locating the magnets away 
from the core in a Shielded area above the D0 level, the 
only weakness in the present control rod deSign has been 
virtually eliminated. What remains is a control rod 
mechanism which has proven very reliable over the past 
eight years of reactor operation, The control rods could 
be removed for repair, or replacement of the absorber sec- 
mon by unlocking ands pulling the entire control*rod assen- 
bly up through the annular ring and out through the top of 
the reactor (after first removing the upper shield plug). 
' Approximately ten control rods would be required to surround 
the core. For a description of the present control rod 
assembly, see Section 1.5 of this thesis, 

2.4.3 Drive Mechanism Through Annular Ring 

Aemnediitcautemeot the comtrol rod mechanism described 
in Section 2.4.2 was designed in order that the control 
absorber could be made to completely surround the core in 
the area immediately adjacent to the fuel elements, inside 
the inner core tank (see Figure 2.4.3.1). As with the 
design of the preceeding section, the Pontrol rod magnet is 
located up within the shielding area of the annular ring. 
However, unlike the preceeding design, it is necesSSary to 
provide a seal in order to prevent D0 from leaking up into 


the magnet area. This seal could take the form of a helium 
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FIGURE 2431 
L-SHAPED DRIVE MECHANISM 
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pressure in and above the magnet area equal to the water 
head below the magnet area. It would be necessary to 
circulate the helium in the magnet area to remove the Dp0 
vapor which would otherwise build up in the helium atmos- 
phere surrounding the magnet area. 

The control rod drive mechanism would be identical 
to the mechanism used in the present MITR as far down as 
the magnet. The assembly below the magnet would consist 
of an "L"-shaped arm connecting the control absorber to 
Die macnes armature. The absorber Section would be pinned 
to the L-shaped arm so that in order to replace the absor- 
ber section it would be merely necessary to unpin the ab- 
Ssorber from the arm and lift it out through the necked 
region of the reactor tank. The elbow of the L-shaped arm 
would be hinged and contain a dual gear arrangement capa- 
ble of raising and lowering the "forearm". With the fore- 
arm in the fully raised position, the entire control rod 
drive assembly could be removed out through the reactor 
foo. It can be noted at this point that should either the 
Glbow of the L-shaped arm or the pin connecting the arm to 
the absorber section fail, the reactor would be shut down, 
As shown in Figure 2.4.3.1, a spring-type dashpot would 
Slow the absorber support arm as it neared the base of its 


travel, and aluminum shields would guide the control ab- 


Perveraceas tLney descend into the core. Eight arced sections 


of absorber would be used to surround the core. 
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2.4.4 Drive Mechanisms Above D,0 Level 

In order to eliminate the complexity involved with 
having the control rod drive mechanisms penetrate through 
the annular ring of shielding, designs were proposed which 
would locate the magnets and drive mechanisms on the reac- 
tor top, entirely above the D50 level in the tank. Three 
possible design arrangements were proposed by this author, 

Tie ims. desien proposed called for horizontally 
placed withdrawal mechanisms installed flush with the 
reactor top (Figure 2.4.4.1). The control absorber would 
be attached to the magnet armature by a cable extending 
vertically to the surface of the D50 column, thence hori- 
zontally to the position of the armature. Eight curved 
absorber plates would form a cylindrical shell completely 
Surrounding the core. The entire withdrawal mechanism 
assembly would be easily accessible through removable 
covers on the reactor top. The cable withdrawal mechanism 
provides excellent protection for the magnet, but suffers 
one serious disadvantage. The absorber will only return 
to the core by the force of gravity. It cannot be driven 
into the core, thus ensuring that the absorber is, in fact, 
eempletely inserted to its furthest point. 

in order to provide for a means whereby the control 
absorber section could be driven into the core, a design 


was proposed whereby the withdrawal mechanisms would be 
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FIGURE 2.4.4 | 
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located directly above the absorber sections. The with- 
drawal mechanisms would be situated vertically atop the 
meadecvor cop shielding and connected to the absorber ee 
tions by solid aluminum rods (Figure 2.4.4.2). As with 
ne previous design, eight arceéd absorber sections would 
eemoletvely encirclestneseeme. The Shaft connecting the 
magnet armature with the control absorber would be compo- 
sed of two sections joined together by a screw-in fitting. 
The two sections could be disconnected at a position 
parallel to the reactor top, thus allowing a cover shield 
fempe Sli1d=into place over the D50 column during fuel 
flipping and fuel changing evolutions. A key at the base 
of the withdrawal mechanism housing and a gear just below 
the armature would mate to allow the upper rod section to 
be unscrewed. The penetrations through the reactor top 
Shielding for the aluminum rods would have to provide 
adequate clearance so as to guarantee that the aluminum 
rod would always be free to travel through the shielding. 
A pressure equalizing line would be required in order to 
eugalize the pressure of the area below the reactor top 
Shielding with that in the area above the magnet. The 
equalizing line would be required to have a cross sectional 
area greater than the cross sectional area of the clearance 
Surrounding the aluminum rod as the rod passes through the 


hegever Lop Shielding, The purpose of the pressure 
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FIGURE 244.2 
WITHDRAWAL MECHANISMS ABOVE SHIELDING 


PRESSURE EQUALIZATION 
LINE 


-. ee eT ie 


mh SSS Short eae) 
ease 


Fae a —————. |= 


fost) Sooo 5 


, oes {Ebon 


ry 
PT LE LE |. PF LAs ve ‘ 





[ ear a 
e 
é 
a 
& 
o 
4 
® 
- 


q* © 


SF nme TYE 
0 
Ls p 
q 
p 
D 


uz 
is 


° 
6 





s 
e e @ 
Pg gL a SS ge I Eg a 


CORE HOLD DOWN 


MECHANISM 


PORT 


«<— BEAM 





D,0 ENTRANCE 








54 


equalizing line would be to ensure that a sudden increase 
ieereactor pressuremwould not exert itself on the bottom 
of the magnet armature, thus using the armature as a pis- 
Bem to drive the control rods upward. With the entire 
control rod withdrawal mechanism located on top of the 
shielded refueling enclosure, all components would be 
readily available for maintenance. 

A third design alternative proposed would employ 
features of both the above designs. Vertical aluminum 
shafts would lead from the eight arced control absorber 
sections to just above the D540 level in the necked tank, 
A horizontal shaft would be coupled to the vertical shaft 
by matching gears. The horizontal shaft would in turn 
iead to a magnetic clutch and thence to the control rod 
drive motor (Figure 2.4.4.3). Upon interruption of elec- 
fiereal current, the magnetic clutch would release the 
horizontal shaft, thus allowing it to turn freely and 
allowing the absorber to drop down into the core. With 
this arrangement, the absorber could also be driven into 
the core should the need arise for such a feature, It is 
a requirement that all of the gears connecting the horizon- 
tal and vertical control rod shafts turn freely at all 
times in order to be effective in a scram condition, Such 
a requirement could be presented as a design flaw. 


However, it was felt that the probability of two or more 
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gear couplings jamming Simultaneously was so small as to 
be nearly nonexistent. Any one of the eight control absor- 
bers would be sufficient to shut down the reactor. 

2.4.5 Drive Mechanisms for Inner Control Elements 

in reviewing the safety requirements for the High 
mux METR, it was decided by #Dr. enc on thauesutficient 
EOncrowmemustumertst CLO maintagan the core in a subcritical 
@Oonaition, even if the core were to be filled with light 
water. in order to satisfy this requirement, Dr. Thompson 
- Suggested that the core should contain an inner ring of 
control elements. The necessity for an incore ring of 
control elements posed a considerable design problem. A 
mechanism which would not interfere with the removal or 
replacement of any fuel element, yet would provide for 
reliable control of the inner absorber elements was re- 
quired. The design proposed by this author is shown in 
Figure 2.4.5.1. The inner control rod withdrawal mechanism 
mela Consist of two wconcentric shafts. The inner shaft 
would be attached at its base to four pivoted support bars 
capable of swinging up or down. The support bars would in 
turn be connected to the outer shaft by four aluminum rods, 
each pivoted at either end. The connection to the outer 
Shaft would be accomplished by a clamp-type collar which 
would ride on threaded grooves cut onto a portion of the 


Pucteteotatte. By rotating the outer shaft, the collar could 
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be caused to move upward or downward, thus raising or 
lowering the support bars. In their lowered (vertical) 
position the support bars would occupy an area above the 
center of the core no larger than the area of the cryostat 
planned for the center of the core. Consequently, the 
mechanism would not hinder the removal of any fuel element. 
In their raised (horizontal) position, the support bars 
would engage slots in the upper Section of the inner con- 
Dol cllenchius, Ihe entire meenanism would be raised in 
order to withdraw the inner control elements from the core. 
(It should be noted that the inner absorber elements would 
require followers in order to prevent the flux peaking 
which would occur if D50 were allowed to fill the space 
vacated by the absorber as the control elements are with- 
drawn.) The withdrawal mechanism would utilize a magnetic 
coupling and be constructed identical to the withdrawal 
mechanisms for each of the absorber sections for the outer 
control elements. 
2.5 Refueling Mechanism 

2.5.1 General 

The initial design concept for refueling and fuel 
flipping was to employ a refueling flask with a self- 
contained flipping mechanism, However, after the concept 


of a permanent, enclosed refueling and fuel flipping 





Dg 
facility to be located directly above the core tank on 


the reactor top was proposed by Mr. Barnett, a new area 
for design was created. This author set about the task 
of developing a refueling and fuel flipping scheme which 
fou1d be used In such anv enclosure. In the process, two 
quite Similar designs were developed, each of which is 
described below, 

Cae enuc Plipger Above Coreslank 

In order to provide a means of rotating a fuel element 
-end for end a rotating table was designed to be mounted 
vertically against a side wall of the refueling enclosure, 
On the rotating table would be a bracket designed to accept 
a fuel element. Directly below the rotating table would be 
another bracket designed to accept the lower adapter unit 
of the fuel element. To one side of the D50 column would 
be a horizontal shute through which fresh fuel would enter, 
and a downward Sloping shute through which spent fuel would 
be discharged into a small D0 filled spent fuel pool. A 
leaded glass window would be provided in the side wall 
opposite the rotating table in order to provide an operator 
with visual access to the refueling cavity during refueling 
or flipping oe A refueling/fuel flipping sequence 
would proceed somewhat as follows (assuming control rod 
drives to be located atop the refueling enclosure): A fuel 
Brabbing tool would be inserted into the reactor tank from 


the reector top, latch onto the desired fuel element, and 
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Peeerc UWOelnce the refueling enclosure. Meanwhile the 
Semurol rod arrive Shafts would be disconnected at the 
surface of the D50 and the upper portion of the drive 
shaft withdrawn into the withdrawal mechanism housing. A 
2 in. thick lead shield would next be slid into position 
above the D,0 column. The lead shield would serve to 
prevent the fuel element from falling back into the reactor 
core area should it accidentally become unlatched. (The 
lead shield would also assist in isolating the refueling 
enclosure from the core while performing maintenance on 
any of the refueling or control rod drive mechanisms. ) 
Mrcm one lead cover in place, a clamping tool would be in- 
serted through the side of the refueling enclosure and 
clamp the fuel element to the rotating table. Once secure, 
fae tool used to raise the fuel element out of the core 
would be disconnected and removed. Another tool would be 
inserted into the refueling enclosure to pull out the pins 
connecting the fuel element to its adapter. Once the pins 
are pulled, the adapter would be allowed to drop down two 
maches, thus giving the fuel element sufficient clearance 
to rotate away. If the fuel element were to be replaced, 
the table would rotate to align the fuel element with the 
chute to the spent fuel pool, the element would be dis- 
charged, and table rotated to receive a new fuel element 


from the horizontal chute. If the fuel element were to be 
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flipped, the table would simply be rotated 180°. The 
remainder of the refueling/flipping operation would be 
the reverse of the above. The lower adapter would be | 
raised and attached to the fuel element, the fuel retrac- 
tion/loading tool attached to the fuel element, the fuel 
element and fuel element adapter clamps removed, the lead 
Shield cover removed from the tank top, the control rod 
drive shafts reconnected, and the fuel element lowered 
ioo its respective position in the core, 

After allowing a spent fuel element sufficient time 
to cool (radiation wise), the spent fuel pool would be 
drained and the fuel element removed through a plug in the 
refueling enclosure Shielding, and thence taken to the 
main spent fuel pool in the reactor building basement. 

2.5.3 Fuel Flipper Away from Core Tank Top 

iMmorecr oom Cleve Lae ivel Pinppine operation from 
the area directly above the core, a modification to the 
design of the preceeding paragraph was proposed. This 
design is shown in Figure 2.5.3.1. The purpose in remov- 
ing the mechanisms involved in fuel flipping/refueling from 
the area immediately above the core to an area which could 
be isolated from the space directly above the core was to 
produce a design whereby maintenance could be performed on 
the fuel flipping/refueling mechanisms while the reactor 


was operating at full power. In this manner it was 
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anticipated that reactor down time necessitated by mainten- 
ance and repairs could be minimized. 

The rotating table would be separated from the area, 
above the core tank by a 6 in. lead shielding wall. A 
removable lead shutter in the center of the shielding wall 
would raise to allow a fuel element and its lower adapter 
to pass from the core tank area to the rotating table. The 
fuel element clamping tool would extend through the center 
of the shaft of the rotating table. The clamping tool for 
the fuel element adapter would penetrate the shielding wall 
just below the rotating table. Instead of being discharged 
into a small spent fuel pool on the reactor top, spent fuel 
would be removed through a refueling door into a small fuel 
element transport flask and transported directly to the 
large spent fuel pool in the reactor building basement. 
Fresh fuel elements would also enter the refueling enclosure 
through the refueling door. Two leaded glass windows would 
be provided to enable visual observation and control of the 
entire refueling/flipping process. One window would be in 
the wall adjacent to the D50 column, facing the rotating 
table. The other would be in the enclosure containing the 
rotating table, facing the refueling door. 

2.6 Underwater Periscope 

While visual access to the refueling enclosure would 

be provided by one or more leaded glass windows placed in 


the walls of the refueling enclosure, it is also desirable 
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that a means be provided for looking down into the core 
itself. Being able to see the core would be a necessity 
during, refueling operations, and in addition would be 
valuable for instructional purposes. It would be possible 
to locate a leaded glass window in the upper corner of the 
refueling enclosure, angled downward in Such a manner as 
to allow one to observe the core. However, with such an 
arrangement one would not get a clear view of the core as 
a result of attempting to look through a gas-liquid inter- 
face at the core tank surface. During reactor operation, 
the liquid (D,0) surface would be turbulent due to the 
rapid circulation of the coolant. The turbulence alone 
would hamper viSion, but in addition it is necessary to 
Shine a light down into the coolant tank in order to see 
the core. Some of this light would be reflected back from 
the coolant surface, thus creating what would be the major 
hindrance to clear observation of the core. 

imporder ee Climninace the difficulties mentioned in 
the previous paragraph associated with observing the core 
by way of the refueling encloSure, a retractable periscope 
was designed. AS shown in Figure 2.5.2.1, the periscope 
is mounted at an angle in such a way that when fully inser- 
ted into the core the outside end of the periscope remains 
mpove the level of the coolant. In this manner the connec- 


mon between the periscope and the vertical viewing tube is 
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Meverequired GO open 2 f1quid interface. When retracted, 
the Cip Of the periscope would be flush with the edge of 
the coolant tank wall, thereby eliminating any interfer- 
ence to the insertion or removal of components in the 
area of the periscope. Although the entire core would be 
visible when looking through the periscope, the mirror 
which extends into the tank would be able to pivot in 
order that one might scan for Some distance up the Sides 
of the core tank wall. The periscope would have its own 
self-contained light source. The light Source would be 
located at the base of the periscope, well removed from 
the intense radiation area, and project its light onto the 
core via a mirror near the periscope tip. One would look 
into the periscope through a vertical viewing tube which 


would extend to a height of 4 ft. above the periscope. 
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Chapter 3 
SHIELDING CALCULATIONS 
3.1 General 

The reactor top shielding design and supporting calcu- 
fae2ons were originally done for the phase one reactor 
fescien. This chapter will first present, therefore, the 
shielding design and calculational work done for a D50 
moderated, cooled, and reflected core in some detail. 

When the core moderator and coolant were changed to light 
water, a new set of shielding calculations was required. 
Since the method and reasoning employed in the phase two 
calculations was exactly the same as that employed for the 
Mmase one calculations, phase two calculations will not be 
reported in detail, but rather a summary of required 
Shielding and calculated dose rates will be given in Sec- 
tions 3.7-3.9. 

In phase one, based on the assumption that the reactor 
top shielding design would consist of an enclosure similar 
to that shown in Figure 3.1.1, it was necessary to deter- 
mine both the composition and thickness of the shielding 
which would be required. As the calculations of this 
chapter will show, no major shielding problems were encoun- 
jmeredq. ouch was nov known in the early stages of the 
design project, however. It was anticipated that shielding 


requirements during reactor operation might be dictated, 
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not by the radiation being emitted from the core but by 

the radiation dose due to gamma emission from te decay 

or the gamma emission accompanying the formation of H? in 
the necked region of the core tank, to mention but one 
example. While the calculations which follow serve to 
justify the chosen composition and thickness of the reac- 
tor top Shielding (shown in Figure 3.1.1), it should be 
noted that the shield composition and thickness were 
originally determined by successive iterations of exactly 
Similar calculations. It is for the sake of brevity and 
clarity that only the final set of calculations is presented 
here. No shielding calculations were performed in the 
rowel Girection since neither reactor power nor radial 
Shielding will undergo any change from the present MITR. 

In fact, radial attenuation will be improved since the core 
boundary will be further from the shielding walls. 


It was required that the DO column combined with the 


a 
Shielding of the refueling encloSure be so designed as to 
be capable of reducing the dose on the reactor top suffici- 
6ntly such that a worker working full time in this area 
would not receive a dose in excess of the maximum whole 
body dose limit established by the U.S. Atomic Energy 
Commission. The Atomic Energy Commission's maximum allow- 


able whole body dose levels are 100 mrem per 40 hour week 


or 2.5 mrem per hour (6). However, it is desirable that 





(0 


the shielding be capable of reducing the dose anywhere on 
the reactor top to well below these maximums. The present 
MITR was designed so that the dose on the reactor top | 
would be less than 1 mrem per hour at 5 megawatts (8). 
Similarly, the goal in the design of the shielding for the 
High Flux MITR was to reduce the dose anywhere on the reac- 
fom woop tO lesaethanei,Ormren per hour when the reactor is 
at a power of 5 Mw. 

The dose rates during both shutdown and operational 
conditions were determined. For the shutdown condition, 
the sources of radiation of most concern are those resulting 
from material activation and fission product decay. For the 
operational condition one must consider, in addition to the 
Shutdown sources, the dose due to gamma ray production in 
fission, fission neutrons, photoneutrons, gamma rays from 
me decay of te | and gamma rays from the production of 
tritium. 

3.2 Calculational Methods 

3.2.1 General 

The calculational methods followed in performing the 
enielding calculations for the High Flux MITR were the 
Same as those used in the design of the shield for the 
Omega West Reactor at Los Alamos Scientific Laboratory 
(see reference 7). This same procedure was also used in 


obtaining the initial shielding estimates for the present 





yes 


Pao. “Radlatyvon measurements recorded subsequent to the 
completion of the Omega West Reactor tend to bear out the 
validity of using the calculational methods employed in 
the shield design. As an example, the calculated dose 
rate due to gamma rays through the lid on the reactor top 
was 114 mrem/hr. The measured value was 107 mrem/hr. 

3.2.2 Gamma Ray Attenuation 

The sources of gamma rays in the High Flux MITR will 
be treated as point sources. The assumption of a point 
source has been estimated by Gage to produce results which 
are at most 20% greater than those obtained from calcula- 
mons PerhermeadsUtckliZingwa distributed source. The error 
involved in assuming a point source of radiation is compen- 
sated for an adjustment to the radiation source strength, 
as discussed in a later section of this chapter. The equa- 
fom for the attenuation of gamma rays from a point source 


nS: 


G = = FBe* y's/em=-sec ’ (3.2.2.1) 
TR 


where, S = source strength (y's/sec), 


R = distance from the source (cm), 

F = fraction of gamma rays in energy range of 
INGerest., 

i = linear absorption coefficient (em7+), 





ce 


xX 


thickness of absorbing material (cm), and 


B= DUOuperactor 


The Omega West Reactor shielding calculations were 
done in terms of decade lengths rather than in terms of 
the Jinear absorption coefficient, and so decade lengths 
will be used here also. A decade length is defined as 
Maat thickness of material, d, required to reduce the 
flux by a factor of ten. Decade length is determined as 


follows : 


1/10 = a 


d 


2.3/u cm = 0.075471 ft. 


The shield thickness in decade lengths, D, is defined as 

— 4 , where x is again defined as the thickness of absorb- 
ine material. It is also necessary to convert the buildup 
fmcecOr, B, from its normally listed value in terms of e- 
folding lengths to decade lengths. Buildup factors in 
terms of decade lengths for water, iron, and lead were 
taken from reference (7) and are given in Figures 3.2.2.1, 
fee. c.c, and eee eerie. It should be noted 
that Figure 3.2.2.1 is for light water. However, since 
Paina ray a2bcenuavion iS most strongly dependent on the 
atomic number of the absorber, and 1,0 and D590 both have 
Pee Same total Z value, Figure 3.2.2.1 is assumed to be 


applicable to heavy water as well. This is a conservative 
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assumption Since D,O is about 10% more dense than H,0. 


2 2 
i should tase be novea that buildup factors for heavy 
concrete were not available, so the buildup factors for 
jron (Figure 3.2.2.2) were used in the calculations 
later in this chapter. This substitution is reasonable 
since heavy concrete is about 75% iron. The result is 
conservative from a shielding point of view because the 


Duildup factor for iron is greater than that of concrete. 


Equation 3.2.2.1 can now be re-written as 


D 2 


y's/em“-sec . (en2aee a) 





Q = Peon 
KarR 

im applying this flux attenuation equation to a shield 
consisting of many layers, it is assumed that the buildup 
mector for the shield is the sum, and not the product, 
of the buildup factors for each layer. Although this may 
not be a conservative assumption, it has been shown by 
Gage to be much more accurate than if the product of the 
buildup factors were used. 

S22 nome Neutron Attenuation 

As with gamma rays, neutrons will be considered as 
Being emitted from a point source. The equation for the 
atvvenuation of fission spectrum neutrons is: 

ge Ae gaa 

AarR 


where N = neutron source strength (neutrons/sec), 





ft 


R 
» 


n 


distance from the source (cm), and 


Hl 


relaxation length for removal of neutrons (cm). 
A decade length for the attenuation of fast neutrons is 


determined as follows: 


-d/d. 


= @ 


Q BIH 


2.3 (cm) = 0.0754 A, (ft). 


Relaxation lengths and decade lengths for some shielding 
materials of interest were taken from reference (7) and 
Bare given in Table 3.2.3.1. 

The attenuation of thermal neutrons is determined by 
Means Of ditfusion theory, The applicable diffusion equa- 
ston for thermal neutrons from an infinite plane source in 


an infinite homogeneous media is: 


pe 2 
C7 oe, neutrons/cm -sec , (ne. 3.2) 


where L = diffusion length (cm), 


D = diffusion coefficient (cm), 
x = thickness of attenuating material (cm), and 
f 
Q, = thermal neutron flux on the core side of the 


oe LC 
3.3 Shielding Calculations, Normal Operation (Phase I 
3.3.1 Fission Gamma Rays 
VAC Pieter yesOuUree Gi radiaviongto be shielded against 


ia che gamma rays produced by fissions in the core. 





ne 





a a a Te 


€19°0 76°8 7° € oe peat 
Gtrt7°O 06°S O°? ed UOIT 
70° T Q°EeT cual a UNUTUNT Y 
Out 2 HL 173°0O L9°T eqtudeiy 
oGL°0 26 °6 66°2 Ol T3YEM 
E09 °O O83. ~ Gas 879T9U0) 
saeqAzeg 
GLi7°O E°9 - My 84V9TOU0N 
UdsABVUYooLg 
9€8°O ‘ail ~ ects a4a10U09 
PUT ALOd 









4 Wd WOZe/SULEQ AYTACLIA 

p ‘ug3ueT YX fUySUeT © “WoT doacessoay 1 jJpoeas 
apeoaq UOTYLCXETAY TeCAOWUSY SATZOAIIY 

DITAODSsOLO IN 


a Let sazey 


T°e* Se STQBL 





tes, 
Assuming 198.22 MEV are released per fission (see Section 
4,2), the source strength of the High Flux MITR operating 


at 5 megawatts would be; 


Sees 10° Watts )(3.1 x fon” fissions/watt-sec) = 1.55 x LOmt 
fissions/sec, 

With each fission event there are, on the average, 7.5 

gamma rays emitted from the fissioning nucleus (9). Thus, 

the source strength of the fission gamma rays from a core 

operating at 5 megawatts is: 


ele 55 x 10+! Bese eas) (75 gam as ea x 1012 gammas 


sec fission sec 


Figure 3.3.1.1, taken from reference (7), shows the 
relative number of fission gamma rays versus energy. The 
straight line extrapolation shown in the figure was used 
to calculate the percentage of gamma rays above a given 
energy. The results obtained are conservative. Gage has 
Shown that for energies above 1.5 MeV the extrapolated 
line predicts about 30 percent more gamma rays than exist. 
The excess in the fraction of gamma rays predicted above 
an energy E in the range O to 1.5 MeV decreases from 30% 
to zero as E decreases from 1.5 to 0 MeV. It is assumed 
that the fission gamma rays can be lumped into three 
eee oe yeoroups. rable 3.3.1.1 gives the calculated frac- 


tions above each of the three energy groups assumed using 
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FIGURE 3.3.1.1 
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on 


the straight line analysis, and also the assumed frac- 
tion at each energy that will be used in the shielding 


calculations. 


aoe oO onl | 


Ca lcilbace Gut ercentace Aaaumed sPereentace 


Energy Above E at E 
O MeV 100% 

1 MeV 33.4% 100% 
3 MeV 3.7% 33% 
6 MeV O. 14% 3% 


Note that the assumed percentage of gamma rays at Bere 
Mave been greatly overestimated. For example, the 
assumed percentage of gamma rays at 6 MeV has been over- 
estimated by about a factor of 20. This overestimation 
of gamma source strength was the method used by Gage to 
compensate for the nonconservative point source assump- 
tion. Gage performed a twelve energy group calculation 
with energy groups from 0-12 MeV in order to check the 
accuracy of the three energy group model. For the water 
and concrete radial shielding of the Omega West Reactor, 
mencound that the three group model as described above 
Was conservative by a factor of 5. Similar conservatism 


could be expected for the D50 and lead shielding of the 
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moon Klux MIiTR. 

The height of the coolant above the core top is 8 ft. 
Since the gamma source is being treated as a point source 
of radiation, it is necessary to add to this the additional 
1 ft. which exists between the center of the core and the 
Gore top. Hence, the total thickness of D50 to be consider- 
ed is 9 ft. The helium-filled space between the D50 surface 
and the bottom of the shielding enclosure is 4-1/2 ft. in 
meiegnc. Ihe refueling enclosure top shielding block con- 
sists of 6 in. of masonite sandwiched between two 1/4 in. 
sheets of boral (neglected for gamma ray calculations), 1 
in. of steel used for support, and 1 ft. of lead. (Note: 
D0 cooling coils will loop through the lead and masonite- 
boral sandwich. ) 


Bore Cclarnivy EQ. 3,2c¢.c.e2 may be rewritten as: 


-(D ptD,+D,) 


at (BytB,+B, )10 (Ceca 


~ lat (R prRytRys ) 

meere the subscripts D, H, I, and US refer to D0; Heliun, 
Iron, Lead, and Upper Shield respectively. For calcula- 
tional purposes, the helium gap will be treated as void. 
mecaces Of attenuation and buildup factors for materials 


of interest as a function of energy are given in Table 


Peo ole. 


_ 
— 
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Seen wire iivomice o.5.1,.1 for the dose rate due 
to 6 MeV core gamma rays on the reactor top directly 


above the core yields: 


18 , 


Nt (94+5+.083+1)x30.48 


De = 7.41 gamma rays /em°-sec. 


For 6 MeV gammas 157.7 eo mea eee mrem/hr., the conver- 
sion from gamma rays/cm ee mrem/hr having been obtained 
from page 7-66 of reference (9). The R.B.E. for gamma rays 
was taken to be 1. Therefore, De = 0,047 mrem/hr. 

Following exactly the same procedure for the 3 MeV group, 


2. = 5,42 gammas /em°-sec = 0.021 mrem/hr, 


and for the 1 MeV group, 


8 el 


5 =A44x 10° gammas /em*-sec = 7.6 x 10 ~~ mrem/hr. 
Summing the above calculated doSe rates, the total dose rate 
Maeve reactor top directly above the core due to fission 
gamma rays is calculated to be: 


IL 0.068 mrem/hr. 


Ye 0,047 x 0.021 + 7.6 x 10° 

3.3.2 Material Activation 
3.3.2.1 General 

Mica priImaryerecacvtOnmWinaicn must be considered is the 


thermal neutron activation of the aluminum contained in the 


mere and Surrounding Structure. The reaction of concern iS; 
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Ree ganna y a Ta, (Co. omen 1) 


fie wvnerna! Chose SeCCtron fOr this reaction is 0.215 barns, 
aize emits an average of two prompt gamma rays for each 
neutron captured by poet (108 

rise Of Concern from a Shielding point of view are 
the gamma rays released in the decay of Nitrogen 16, which 


is formed in the reaction: 


CL ese + p (eres 22) 


Although the activat*on cross section averaged over the 
Molre ticaton= spectrum for this section is only 2O micro 
barns, it is usually the most important reaction associated 
with the reactor coolant. Nitrogen 16 decays with a half- 
mire of 7.35 sec., liberating gamma rays of energy 6.1 
MeV (75%) or 7.1 MeV (7%), plus one prompt proton which has 
no significance from a shielding point of view. (11) 

The final material activation reaction to be consider- 
file oois Section will the formation of tritium, the reac~ 
tion of interest being: 


oer BF) uo y (Stare. Lsah) 


whe thermal neutron cross section for this reaction is 0.57 
milli-barns, and the resulting prompt gamma has an energy 
of approximately 6.26 MeV (9). 

pce oe ny) aac” 


In order to determine the dose rate on the reactor 


28 


top due to the decay of Al”-~, the reactor core area was 





divided into three regions as shown by the dotted lines 
Mab ieure 3.5.c.c,1. "Gamma rays born due to the decay 


of are8 


were assumed to be point sources located at the 
center of each of the three regions. The volume of 
aluminum in each of the three regions was determined based 
fama Dest estimate of utlimate core and supporting struc- 


more Composition: 


Region I = S506 ane 
Region II = 56,700 em= 
Region III Oi OO em? 


The aluminum volumes by region were based on assumed alu- 
minum volume percents within each region of 2% for region 
I above the core, 43% for region II in the core, and 4.3% 
for region III below the core, plus a 1/8 in. thick 100% 
aluminum tank wall forming the boundary of each region. 

ine average thermal flux in each of the three regions was 


conservatively estimated based on the preliminary core 


physics calculations described in Chapter 4 of this thesis: 


Region I G = oe neutrons/em*-sec 

Region II G@=5 x On neutrons /cem~-sec 
i 

Region IIT | G = one neutrons /cm*-sec 


OmiyetnesciecrmatencucCrem induced flux is considered 
mere due to the rapid decrease in the reaction cross sec- 
tion of aluminum at higher energies. The activation cross 


a 
Y 


section for 0.2 MeV neutrons is down to about 10 ~ barns, 





FIGURE 33221 Z 
REGIONS FOR Al*‘, (n,x) Al®” CALCULATIONS 


bases 
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for example. Assuming an irradiation of infinite duration, 


tne number of meulron ceptures per second in acl is 
described by the equation: 
Captures/sec = NQ-, .,9V,, (3, Selec ) 


Proceeding with the calculation for region I, 


N Oat ?Mar = (0.06023) (0.215)(10°) (5500) = 


7.12 x 107° 


captures/sec. 
pam lariy for region II, 


16 


NO 4 PV = 3,66 x 10 captures/sec, 

and for region III, 
Nieieee OV = ol. ol x 1026 captures/sec 
act” Al ‘ 7 ' 


ime umdl LOnmuor aoe 3.3.2.2. 1 was taken from 
references (7) and (9), and lists the number and mean 
energy of the gamma rays emitted per thermal neutron cap- 
ture in ae along with the gamma flux to dose rate con- 
version and the decade lengths in D0 5 iron and lead 
mespectively. 

Equation 393.1.1 was used to obtain the dose rate on 
Pee reactor top directly above the core due to the thermal 
neutron activation of Al“’, Table 3.3.2.2.2 lists the in- 
put data not previously given required to solve the equa- 
tion, and the resulting dose rate on the reactor top due 


to gammas from each region. The calculational procedure was 
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exactly the same as the Sample calculation done in Section 
3.3.1, and so is not repeated here. Summing the dose 
rates listed in Table 3.3.2.1.1 gives a total dose rate on 
the reactor top due to the thermal neutron activation of 
al! to be 0.064 mrem/hr. 


eo eG 


Merch eedeee Came forte restits in hich radi- 


ation levels around exit coolant piping, there is normally 
feet ie reason to consider the cee (ap) nie reeeuLOn fer 
purposes of reactor vessel shielding due to the small 
activation cross section of to, However, the fluid flow 
sememe for=tne fren Fiux MITR calls for the coolant to 


emrcuiave up to very near tne top of the necked portion of 


the D50 tank. It is for this reason that the dose rate due 
iLO yt decay is calculated. 
The threshold energy for the oie fa co ) te reaction is 


very high (approximately 10 MeV). Therefore, it is reason- 
able to assume that te may only be formed on the first 


Pollision of a fission neutron. The probability of a fis- 


ING 


Sion neutron interacting to form N ioe aioky bieeraclioO sof 


16 


the cross section for the formation of N to the total 


memoval croSs Section seen by a fission neutron. Table 1.1 
am Rockwell gives the effective removal cross section for 
fission neutrons in D,0 as 2.8 barns per molecule (12). 


iG 6 


mae probability that an Oo (n,p) yt reaction occurs in 
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iC > .eowao lon OL a Neutron 1s therefore, 


20 x 107° 6 


P(n,p) ~ eee ar en ele ok ue 
where it has been assumed that all neutrons have that 


mie comictonmyitheitner Ho or 0-°. ‘This assumption 


is conservative and leads to a slight overestimation of 


mio). Sines rate of production of yO Pay Cie enc re 


reactor becomes: 


R= 4 fission neutrons x P(n,p) 


# fission neutrons = (5 x 10° watts ) 
NOutiTssLoOns neutrons 
(3.14 x 10 wapeeaee) co iccion 
4F fission neutrons = 3.9 x 10 neutrons/sec. 


-6 12 
) 


R = (3.9 x 10°")(7.15 x 107°) = 2.78 x 107* activations/ 


16 


Sec. it iS now asSumed that all the N atoms formed are 
mormed in the D509 Welen cirentaves wo through the necked 
core tank, This assumption is conservative because some 
N° will be formed in the reflector D0, but the Nuoecee! 
mache reflector will make a negligible contribution to the 
dose rate on the reactor top. The volume of the D350 in tne 
inner coolant tank is about 225 gallons. Therefore, at 
1850 gallons per minute the coolant takes 7.5 seconds to 
pass from the inner tank inlet to exit. The external sys- 
tem D,0 volume is approximately 475 gallons. Therefore, at 


1850 gallons per minute the total cycle type for the D0 
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16 


is approximately 23.3 seconds. Now assuming N”” to be 
formed uniformiy throughout the height if the inner core 
tank (an obviously conservative assumption since few 

fics lon neutrons will undergo their first collisions near 


the top of the necked D,O column), and taking the activity 


2 
of the entire volume of D0 in the inner tank to be that of 


the coolant as it exits the upper portion of the inner tank, 
(again a very conservative assumption) the saturation 


16 


activity due to WN decay becomes: 





- At, 
A p Ome oe 
sat a _ AT, ( ® e e e ) 
j-e 
where \ = decay constant for te = .693 = ,O0943 sec 
4 : ° 7.35 sec : 
tC. — i fpediatvwonmebpime an the inner core tank 
ty = Oval cycle time 
Asat = (2,08 x 10°") Te (0.0943) (23.3) 


a 


A sat = 1,59 x 10°~ disintegrations/sec. From Section 


HO 


3.3.2.1, it is recalled that N”-~ decays by emission of 


either a 6.1 MeV gamma (75% of the time) or a 7.1 MeV 


gamma (7% of the time). It is assumed in the following 
ie 


Salculations that N-~ always decays by emission of a 6 
MeV gamma, (This assumption is conservative since in 


4 


lead a 6 MeV gamma is more penetrating than a 7 MeV gamma. ) 


16 


The dosc rate on the reactor top due to N”” decay is 





Qh 
calculated assuming the total saturated activity to be a 


Dome oource=!Ocaced auewchne Suriace of the DO column. 


e 
Meo SWeStituting into Eq. 3.3.1.1: 


12 
g = —(1-59x10*)(1.0) (1. 2414,6) 107 (-26+6.71) 


Me =(5+,0834+1)30.48 


G = 6,32 y's /em°-sec : 


feererore, 
Dose Rate = 0.04 mrem/hr. 


The dose rate due to nto 


decay calculated above is 
ultra-conservative from a shielding point of view. (The 
dose is overestimated by about a factor of 10.) The 
primary purpose of the calculations was to show that even 
Being Overly conservative assumptions the decay of nto 
would not make a dominant contribution to the total dose 


meee On the reactor top. 


Sucre ha (nu, 7) ao 


since only thermal neutrons are important in the 
formation of tritium, the gamma ray flux resulting from 
tritium formation will vary, as a function of position in 
mre cOre, directly with the thermal neutron flux. Tritium 
decays by emission of a beta particle, which, although 
fel y eieiiiicant Irom a health physics standpoint should 


mube ingested into the body, is insignificant for purposes 





J 


Seer cacvct es UG) ol Cm ine ecalculatvons. for calculational 
Purposes, the D,0 volume was divided into four regions as 
shown in Figure 3.3.2.4.1, and the gamma rays from each 
region were assumed to originate from point sources loca- 
ted as shown in the figure. An average thermal flux was 
assumed for regions I and II based on a conservative 
interpretation of the results shown in Chapter 4 of this 
aceic. tne average thermal flux in regions III and IV 
was calculated by considering the necked core tank to be 
a D0 thermal column in which the calculational method to 
be shown in Section 3.3.4.2 applies. 

The number, n, of gamma rays produced per second in 


each region due to the formation of tritium is given by: 


2 a 
n=N Fact Vp, ‘ (Gee eee an) 


24 


where N = deuterium atom density = 0.06646 x 10 atoms/en*, 


demivacion ¢€ross section for H°(n, y)H° reaction = 


Tact = 
. 5 (es im barns, 
Q@ = average thermal flux, 
Vy ey = coolant volume of interest. 
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REGIONS FOR H* (n.¥) H” CALCULATIONS 


ee eee oe ee eee ee oe 








O7 
Topmemen ooo. 1 


om? 





2 2 
Rezion Volume neutrons /em~-sec rammas/Sec 


i 6.0 x 10° 1.0 x 10! 1.14 x 105 
II 5.4 x 10° oreo 1.02 x 10+ 
ie 1.95 x 10° 1.8 x 101° 6.65 x 101° 
IV 1.95 x 10° Es xo° 1.96 x 10° 





Again the calculational procedures outlined in Section 
3.3.1 are followed in order to obtain the dose rate on the 


3 reaction. Since a 6 MeV 


reactor top due to the ae (ey) a 
gamma iS more penetrating in lead than a 7 MeV gamma, the 
reaction gamma rays are taken to be 6 MeV (a conservative 
assumption). Table 3.3.2,4.2 lists the results of the cal- 


culations for ganma flux and dose rate on the reactor top 


due to the H* (9) H? reaction. 


Vapilemoe sees. ec 


Region (cammas /em“-sec ) Dose Rate (mrem/hr) 


i 0,246 1.56 x 10°° 

TI 0.026 noe ox 107" 
ror 0.0056 nS caer 
IV 0.000036 2.30 x 107f 
Total 0.278 1.76 x 107° 


ee ee ee ee (SRR aR RR SN SRS RE a 





AS shown by the above calculation, the dose rate on the reac- 


fer Lop due co tritium formation is nesligible. 
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3,3. Cue chOnerroduct. Decay 

As one might expect, the dose rate due to fission 
product decay is small compared with the dose rate from 
prompt fission gammas. An infinite irradiation time is 
assumed, thereby yielding a dose rate which contains the 
PenurtvuolOn Of DOLN Short and iong-iived fission products, 
ie results So obtained are conservative Since in practice 
the reactor will operate during the week and be shut down 
on weekends. 

According to Rockwell, the rate of energy emission by 
gamma rays after an infinite operating time at 5 MW is 1.14 


nd 1018 


HeV per se@ond (ea. if the average decay gamma ray 
energy is taken to be 2 MeV, the decay gamma ray source 
strength becomes 5.7 xX oa gammas/sec Ese A point 
source iS assumed at the core center, and the calculational 
pmaeeecdure outlined in Section 3.3.1 is once again used. The 
ieee tace On Lhe reactor top due to fission product decay 
me calculated to be 4,1 x 107" mrem/hr. 

3.3.4 Neutron Attenuation 

3.3.4.1 Fast Fission Neutrons 

Pie wselacliOmmaor tie mabyenvacien of fission Spectrum 
neutrons is fiven by Eq. 3.2.3.1, here repeated in slightly 
Secrerent form: 
S a(n r.) 


= ; (3.3.4. 152) 
MarR 





SIS) 


where S = neutron source strength (neutrons/sec), 


R = distance from the source (cm), 
x = thickness of attenuating materials (cm), and 
r.. = relaxation length for removal of neutrons (cm). 
The relaxation length, Ny is calculated by the relation: 
me (em) , (cece 12) 
where Q7, = micruscopic removal cross section (barns /atom), 
N. = Avagardo's number (6.04 De 10°° molecules/gm- 
mole), 
p = density (ems /em”), and 


A = molecular weight (gms/gm-moie). 


= e 
r (0, 604x10°") (1.105)(2.8 x 107°") 


The relaxation lengths for lead and iron are given in Table 
4.2.3.1 to be 8.94 cm and 5.90 cm respectively. 
Memealculaued Jn Section 3¢3.) there will be 1.55 x 
10°’ fissions/second in the High Flux MITR operating at 5 
MW. Taking the neutron yield per fission to be 2.47, and 
considering that one of these neutrons is needed to sustain 
the chain reaction, there will be a total of (1.47)(1.57 x 
no: ! ) A eas lone Pisstonm, Meucrons per second available for 


Meakase or for absorption in the various reactor :naterials. 
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The conservative approach is to assume that all 2.31 x 
et neutrons leak from the core. 
cup wien into wd. S.3.4.1.1, tmerfast fission 


neutron flux on the reactor top is calculated: 


~ (2X30.48 in .083x30. coer 1x30.48 
Or 


en Ce eo ‘aac 


Ar [(O45+. Coe ee O, 48] - 
@ = 1.69 x Cia neutrons /em@-sec. 


From Figure 2.3 of reference (12), the conversion from fis- 
Pon neutron flux to dose rate is taken to be 8.6 neutrons/ 
om~sec per mrem/hr, Therefore, the dose rate on the reac- 


wer top due to fast fission neutrons is calculated to be, 
1.69 x 107° ae 


= 1,06 x 10 ~ mrem/hnr. AS one might expect, the 
Zoverission neutron flux on the reactor top iS primarily due 
wo une excellent me@derating characteristics of the D50 sie 
mac necked core tank region, Using the above calculational 
mevnod, it was found that the fast neutron flux at the cool- 
ant surface is only 1.96 neutrons/cm@-sec. 

Not included in the above calculation is the moderating 
erfect of the masonite. in practice, most of the fast neu- 
trons reaching the masonite are thermalized in the masonite, 
and subsequently absorbed in the Surrounding boral. ‘To 


neglect the masonite, therefore, yields a conservative 


mesSult. 





one 


3.3.4.2 Thermal Neutrons 

In order to determine the dose rate on the reactor top 
aue to thermal neutrons, it is first necessary to deternime 
the thermal neutron flux at the surface of the necked cosine 
tank. This can be done by assuming the necked coolant tank 
region to be a cylindrical thermal column. The thermal flux 
at the base of the necked coolant tank is conservatively 
estimated, based on the results of the flux calculations in 


le neutrons /om-see, For 


Chapter 4 of this thesis, to be 10 
meee oebloulOn Of Inverest Mere, Tne apolicable relation for 
faieer lux as a LunNectiOn Of axial position in the necked core 


tank is: 


ws ~¥2 
Z(7) = Be (eine 42. 1) 
where Z = axial distance above the necked core tank base (cm), 


(2) 4 (22405) (em™hy , 


emacs vonwrenctn==— LOO cm in D50 ; 


Peete ( L r.. 





il 


4 
ll 


a 
i 


lI 


extrapolated tank radius 


=” 
cr 
hy 
i 


transport mean free path = 3D 


D= diffusion length = 0.80 cm in D,0, and 
r = core tank radius = 25.4 om. 
substituting in numerical values, 
~1 


y, tbo) iat 2-05, = 0,0842 om 


CZ) aca =Ye (0. obs) (2) neutrons/em@-sec. 





C2 


Calculated values of @(Z) are given in the following table: 
Rape os. Settee. L. 


Z (cm above base of 


necked core tank) (2) neutrons /cem*-sec 


Lg 1.8 x 102° 
TAS 5.3 x 10° 
193 (coolant surface) on oe. 10+ 


Because lead is a very poor attenuator of thermal neutrons, 
and because it is desirable to avoid the capture gamma rays 


5 


which are produced in a Fe 6 (n; 7) Fe?! reaction, it was 
decided to line the interior of the refueling enclosure with 
boral. A 1/4 in. thick sheet of boral will result in a 
maximum fractional transmission of thermal neutrons of 5 x 
107? (14). Consequently, if the refueling enclosure is 
lined with two 1/4 inch sheets of boral, the thermal neutron 
flux will be effectively eliminated before entering the 
iron or lead shielding. 

3.3.5 Photodisintegration of Deuterium 

The threshold energy for the ia (ya) ie ee is 
e.c.3 MeV. The reaction cross section increases from zero 
at the threshold energy to a maximum value of 2.24 millibarns 
for a gamma energy of 4 MeV. The cross section then decreases 
to a value of 1.2 millibarns for a 13 MeV gamma (11). For 


calculational purposes, it will be assumed that all gammas 





nO 


of energy greater than 2 MeV are capable of initiating the 
(y,n) reaction and that the average reaction cross section 
m5)1,5 millibarns, 

In order to determine the number of neutrons generated 
Pema tuneviom of pPOSition in the core, it is necessary to 
determine the flux of gamma rays above 2 MeV as a function 
®t POSsition in the reactor tank. For this purpose, the tank 
was again divided into four regions as shown in Figure 
Pes, c.o.i. The gamma flux in a given region is that due to 
gammas born in the region plus that due to gammas entering 
the region from each of the Surrounding regions, Except 
for fission gammas, the number of gamma rays per region of 
- energy greater than 2 MeV was taken from the results of 
PrecwioLli@ns Pperiormed in previous sections of this chapter. 
Table 3 on page 50 of reference (7) shows that 19.3% of all 
fission gammas have energies greater than 2 MeV, and so this 
percentage was used rather than the asSumed percentage given 
in section 3.3.1. The flux entering a region from each of 
its three neighbors is calculated using Eq. 3.3.1.1 where 
the gamma rays born in each region were treated as point 
peurces, and attenuation to neighboring regions is based on 
point source to point source distance. Mean gamma energy 
is taken to be 3.5 MeV (a very conservative estimate). 


Results obtained are listed in the following table: 
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Having determinedthe average gamma ray flux per 
region, the number of neutrons, n, produced per second is 


given by the relation: 


n= No OW neutrons/sec, (ee, 5c 
where 
N = atom density of H° = 6.6 x 10°" atoms /em®, 
We activation cross Section = 1,5 x et om”, 
a = averaze region gamma flux 


Ve= Volunerc: D50 in the region. 


eke Se 


Region Photoneutrons, n, born per second 
r 1,09 x 1014 
eT 1.07 x 10l* 

ita 8.4 » iol 
IV 1.25 x 107° 


The neutron and proton released by the Ho (y,n) H reac- 


tion share equally the difference between the gamma energy 
and the threshold energy. Assuming an average reaction 
gamma energy of 3.5 MeV, the average energy of the neutron 
released is 0,64 MeV. Conservatively, the aaa 
ee Che pnovoneuvtrons born in the coolant tank can be 
approximated as a fission spectrum, The neutron flux on 


Pewrcaever top is therefore calculated by substituting 
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io Loeeeeeniieeeto the results obtained from Eq. 

Bao. ae ARN eluded Meuvrom aveenuation due to dis- 
tance, D,0, iron and lead was multiplied a factor of 107* 
to account for the attenuation due to the 6 in. of mason- 


ite (16). The calculated results are summarized in the 


following table: 


Tap le oro. 5. 


Region Neutrons /em”-sec Dose Rate (mrem/hr) 


I 2.8 x 107!9 oto 
II 6.5 x 107° 7.6 x 107° 
ot 1 OO: 0.004 
IV 0.09 0.010 


Tote L 0.12 0.014 


— 


As the above tabulated values show, the addition of 6 
Pe Onomasonite is Sufficient to reduce the neutron dose 
rate on the reactor top to a very low level. From a dose 
meee POINnt Of yiew, it would be possible to reduce the thick- 
ness of the masonite. However, since masonite is relatively 
inexpensive, the reduction of neutron induced activation of 
reactor top structural materials resulting from the addi- 
tional masonite shielding is well worth the small addition- 
eeeexpense, A lew cemviseys leveleissparticularly desirable 
in planning for the eventuality that one may want at some 
future time to disassemble the reactor top shielding enclo- 


slr e , 





LOWE 


3.4 Shielding Calculations, Reactor Shutdown (Phase TI) 

3.4.1 General 

The predominate radiation sources in the reactor after 
shutdown are: (1) fission product decay gammas, (2) gammas 
Puen une 1MGduUcCea@ractivity Of the structure and the coolant, 
and (3) the photodisintegration of deuterium. Immediately 
after reactor shutdown, radiation due to delayed neutrons 
Should also be considered. 

Quite obviously, the magnitude of the shutdown radiation 
Sources ts Only @ small fraction of that for the reactor 
Operating at power. Consequently, a calculation of the dose 
rate on the reactor top after shutdown would, under normal 
circumstances, yield a minute result which would have little 
meenwme. There are, however, two important situations in 
which the dose rate on the reactor top due to shutdown radi- 
ation may be quite significant: (1) when reactor top 
ee HOin= io removed from the area directly above the core 
tank, and (2) immediately after shutdown when all coolant is 
lost from the necked core tank. 

3.4.2 Reactor Ton Shield Plug Removed 

ohielding directly above the reactor core tank would 
mot normally be removed except for purposes of reactor main- 
menance or repair. Assuming, for instance, that the reactor 
mo shuc down at 5 p.m. on Friday, work in the area of the 
reactor top would not commence until Saturday mornings, a 


cooling time of at least 14 hours. During this time period, 
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virtually all the activity due to Ms (Ty /o = =e eel), 
and nyt6 (y/o = 7.35 sec) decay will have died away. 

From Figure 3 of reference (13), the total gamma energy 
decay rate after infinite operation followed by a 14 hour 


LO 


decay period is approximately 1.5 x 10°~ MeV/watt-sec. 


Following operation at 5 MW this figure becomes 7.5 x 1016 
MeV/sec. Again taking the average gamma energy to be 2 
MeV, the fission product decay Source strength is 3.75 x 
1016 gammas/sec. 
Equation 3.3.1.1 is again used to determine the gamma 
flux on the reactor top, where in this case gamma attenua- 
tion is due only to distance and the coolant in the mecied 
core tank. The calculated result is @ = 1.03 x 10? gammas / 
em--sec, which converts to a dose rate of 303 mrem/hr. 

in order to determine the number of photoneutrons pro- 
Guced by the fission product decay gammas, it is necessary 
to determine the flux of gammas in excess of 2.23 MeV as a 
Bene GlOn Of POSition in the coolant tank. Fisure 3.2b of 
reference (12) shows that gammas in excess of 2 MeV comprise 
roughly 3% of the total decay gamma source strength 14 hours 


after shutdown. The decay gamma source strensgth capable 


oe initiating the Ho Gyan) ih reaction is therefore, 





. iol - 
AO ee 107” MeV/sec) _ 9.8 x 1014 gammas/sec, 


M a /eamma 
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The core tank was divided into three regions (shown in 
Figure 3.4.2.1) and the average gamma flux per region 
calculated using Eq. 3.3.1.1. Knowing the average flux 
capable of inducing the (y,n) reaction, the number of photo- 
neutrons produced in each region was calculated using Eq. 
3.3.5.1, where for fission product decay gammas the average 
activation cross section was taken to be 0,001 barns. 


Results obtained are summarized in the folloving table: 


Wapile ott, e, 1 


Gamma Flux Neutrons Produced 
Region (gammas/cme~sec } (neutrons /sec) 
I CeO. ee non 
II 2.23 x 10° 2,88 x 107 
TIT 7.92 x 10? 1.02 x 10! 


a ne 


Taking the average activation gamma ray energy to be 2.3 
MeV, the average energy of the photoneutron is calculated 
ie be; 


2 OU ee. 
e 


= Sone V 

Neutrons of such low energy will be thermalized very near 
eneir birthplace by the surrounding D0. The photoneutrons 
generated due to decay gammas are therefore considered to 
be thermal neutrons, and the method of Section 3.3.4.2 is 
used to determine the neutron flux on the surface of the 


Pere Cank. Having reached the Surface, the neutrons are 


assumed to radiate as a point source out through the opening 





FIGURE 3.4.2. | Lio 
REGIONS FOR PHOTO-NEUTRON CALCULATIONS 





a 
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Pee wGeae Comoro tein salhe dose rate on the reactor 
top due to photoneutrons is calculated to be 0.03 mrem/nr, 
negligible in comparison to the dose rate due to decay 
gammas. 

3.4.3 Coolant Lost from Necked Tank 

Immediately after reactor shutdown coincident with the 
Meo o) coolant in the necked core tank, the primary sources 
ereroctotlOneWwoula be thesdecay of fission products in the 


28 


core and the decay of Al in the core and surrounding 


28 


merece trom SeCiaen Googe. l, the initial Al rate of 


decay will be OT x 10t© gammas per second, the decay 
gammas having an energy of 1.8 MeV (taken as 2 MeV for cal- 
culations). The total scurce strength immediately after 
shutdown is, therefore, 6.2 x ou gammas per second. The 
pemma tlux on the reactor top is calculated using Eq. 
3.3.1.1, where in this case gamma attenuation is due only 
to distance and the iron and lead shielding on the reactor 
top. The calculated result is 9 = 1.3 x 10° eammas /em”- 
sec, which converts to a dose rate of 383 mrem/nhr. fter 
14 hours, should coolant still be absent from the necked 
mete tank, the dose rate on the reactor top will have 
decayed to approximately 23.2 mrem/hr. 

In Section 3.3.4.1 it was assumed that 2.312 x lov! 
@meuerons per Second leak from the core during normal 


operation. ASSuminz the delayed neutrons fraction to be 
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007 1. OSex oe of these neutrons are delayed neutrons. 
Therefore, at the instant of reactor shutdown the core 


15 neutrons/sec. Using the 


will continue to emit 1,62 x 10 
method of Sections 3.3.4.1 and 3.3.5, the neutron dose on 
the reactor top at the instant of shutdown is calculated, 
where in this case neutron attenuation is provided only by 
distance and the masonite, iron and lead shielding. The 
result is 18.5 rem/or. Fortunately, the half-life for 
delayed neutron emission is small (0,43 - 55.6 seconds) and 
so after about 15 minutes the dose rate on the reactor top 
would be negligible (15). 

Needless to say, should coolant be lost from the 
nected core tank, personnel should stand well clear of the 
area, immediately above the necked core tank, 

3.5 Dose Rate Through Concrete Shielding (Phase 1) 

Thus far in this chapter all shielding calculations 
have been performed for radiation passing up through the 
necked coolant tank and out through the reactor top shield- 
ing enclosure. Although one can reasonably expect that the 
6-1/2 feet of dense concrete which shields the reactor top 
work area surrounding the necked core tank will allow only 
minute radiation levels on the reactor top, one brief and 
approximate calculation is in order. An estimation of th: 
Pequired denaivy for the reactor top shielding concrete 


Will also be made, 
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JMeter omonwercxaens the thermal neutron flux, is 
assumed to be emitted from a point source located at the 
core center. Radiation source levels as a function of 
energy were taken from the preceeding sections of this 


ehapver. 


Table 325.1 


Gamma Ray Source Strength 


Energy (MeV) __{gammas/sec ) oa 
1.74 x 1026 
1.14 x 107° 
4.56 x 10° 
2460 % 1026 
3.87 x 107! 
O550.x lo 


1.17 x 1078 


= no Oo FP HWA N © 





Fast meutron source strength = 2.31 x oe neutrons /sec. 
Thermal neutron flux at the tank wall = ie neutrons / 
em--sec. Gamma ray attenuation through the D5 and concrete 
Pameldime 1s calculated utilizing Eq. 3.3.1.1 where for the 


current case: 


-(D+D,) 


(BytB, )10 ; (CS eeaals) 


Ee S 
Mar (RytRa) © 
monere tne subscripts D and C refer to D0 ena Comer euc 
respectively. Decade lengths as a function of energy were 


taken from reference (7). The calculated results appear in 


Bae Ol rowing tabie: 





cee! 
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Equation 3.3.4.1.1 is used to determine the fast neu- 
fron ilu On tae reactor top. Converting the exponential 


mecenuetion coclficients to decade lengths, 


da =POnCi>4e ON a=40,.806. ft. , 


D50 


ee nate = 0.8 ft. (experimentally measured value for 


O,WeRe concrete shielding, specific 
gravity = 4.6). 
oS ianeintoubq. Gvs.4.d.2, 


2.31 x lot! 


2 cme. Semabcs 
bor [(34+6.5)30.48] © 
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G = 0,242 neutrons /em“-sec : 
™herefore, Dose Rate = 0.028 mrem/hr. 

Thermal neutron attenuation through the concrete 
Shielding was calculated using Eq. 3.2.3.2, where for con- 
crete with specific gravity = 4.5 Gage has determined that 
er cic, , ad ) =O,3/e cm. scubstitucing into Eq, 
Pee. te yields: 


B= 2,58 ge %/1-91 (3.5.2) 


A thermal neutron decade length in concrete is equal to 


0.0754 x (L) = (0.0754)(1.91) = 0.144 ft. Therefore, 


G = 2,58 pO, ee (2 1012) (197%/0- hy 
) = 
ie (ig =) =e year 6.94) 
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In order to reduce the thermal neutron flux to 1 neutron/en<- 


) 
=EGi= = Oita - 


Sec, & concrete thickness of x = 
quired. Thus, nearly all the thermal neutrons from the tank 
Will be absorbed in the first 2 ft. of concrete. Beyond the 
fics 2 tc, ohne thermal neutron flux results from the slowing 
down of the fast flux, and is thus approximately equal to 

toe fast flux at each point. 

PecnouUsmrne CalculabionssoL this section indicate a 
dose rate through the concrete shielding of 0.03 mrem/hr, 
thereby indicating that concrete of specific gravity = 44.5 
is actually more dense than required, it is felt that dose 
Face uncertainties attendant with the various penetrations 
enrough the upper concrete shielding are sufficient to 
warrant the use of the s.g. = 4.5 concrete. 


3.6 Shielding Calculations Summary 





3.6,1 Dose Through Shielding Directly Above the Core 


source Dose Rate (mrem/hr) 


Core Prompt Fission Gammas 0,068 
ae Cin.) azo 0.064 
to Decay 0.04 
H° (n,y) 0.0018 
Fission Product Decay Gammas 0.0004 
West Fission Neutrons OSG 
Thermal Neutrons 07005 
Photoneutrons 0,014 


Low onl: 0.195 mrem/hr. 
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Geen wetocce tL nous Reactor Top Concrete 


Source Dose Rate (mrem/hr) 
All Sources Combined 0.0018 mrem/hr 


3.6.3 Shielding Conclusions 

The maximum dose rate anywhere on the High Flux MITR 
top is calculated to be less than 0.2 mrem/hr. This fig- 
ure is a factor of 5 less than the design goal, and a fac- 
wor @n Ic,50less than the limit set by the ABC. The 
megevor Lop shielding, as described in previous sections 
Beechis chapter, iS therefore considered adequate and safe. 
ber echielding Calculations, Normal Operation (Phase IT) 

Shige Grete ee 

The phase IT design differs from the phase I design, 
for shielding purposes, in three basic ways: (1) the 
water column above the core is light water instead of heavy 
water, (2) the water column above the core is 11 ft. high 
vice 8 ft., and (3) the shielding lid on top of the water 
column is simply a flat cover instead of the "doghouse" 
Py@wesire Shown in Fisure 3.1.1. In going from phase I to 
phase II, the magnitude of the shielding problem was 
greatly reduced. No longer need one be bothered by tritium, 
either from a health physics or a shielding standpoint. 
Light water makes a much better neutron shield than heavy 
Per Cle cOmlus Much elarger neutron absorption cross sec- 


PEO, Wem in addition does not act as a source for 





Aes 

Pievoncutronom™ by iterative calculations, it was deter- 
Poemed tna aewead Shieid lOwinches thick on top of the 
water column would provide adeauate shielding. The 
following calculations serve to justify this choice and 
Peeandicave the expected dose rate on the reactor top. 

3.7.2 Fission tamma Rays 

Table 3.7.2.1 presents the calculations for the leak- 
age of fission gamma rays through the lead shielding lid on 


the tank top. 


TUOve 3, fee. 


1 Mev 3 Mev 6 Mev 
Material Thickness D B D B D B 
0 Heme el) ORO G3) 12,3 44 5.4 
Lead Omi 8.8 Soa 6.4. bon loge 
ion 20,0 mm idowsll;5um 10,7 10.0 ~ 15.6 
G(r) (mo attenuation) 1 x Ope foie. Bex 10.” 
y's/em -sec 
% (attenuated) L aeeee 107° Ton 46.8 
; a 
y's/com™-sec 
Conversion Factor 700 270 158 
Dose Rate (mrem/nr) 2.5 x 10° 0.072 0,292 
Total Dose Rate = 0.364 mrem/hr 
3.{.3 Material Activation 
26 


3.7.3.1 ar?! (n,y) ale? 


Pree section 3.3.2.2 it is seen that the only signifi- 


SEaiGecOlvribu,ioOn FO the reactor top dose rate is made by 
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6 and 8 MeV gammas from regions II and III. Collapsing 
the paige! (new) aie? source strength to a point source lo- 
cated at the core center, a source of 5 x 1026 gaminas/sec 


is obtained. Calculated results are given in the follow- 











ing table. 
Tap emo. (sor. 1 
6 Mev 8 Mev 
Ma Lerial “Thickness D B D B 
H,0 ee AH 5.4 a78 4.0 
Lead eee e5565 On 2 Se On 
LOtals lero SRO CD Loo 
@(r) (no attenuation) 3,42 x 1079 y's /em--sec 
Fraction Ong OrS> 
Q@ (attenuated) y's /om°-sec ie 53 
Conversion Factor 158 eS 
Dose Rate (mrem/hr) O.071 0.415 


Total Dose Rate = 0.485 mrem/hr 
ace tee Mcny ) H 
the production of gamma rays from thermal neutron cap- 
ture in hydrogen was not treated in the phase I calcule- 
tions, and So will be described here. One gamma ray of 


1 (n,¥) H® reaction, 


/ 


energy 2.23 i.eV is emitted in the H 
Biienteacion Creve section being 0,33 barns (7). In order 
to Simplify calculations, the water in the core was con- 


mecerecO tO De a sphere 1 fl. in radius and the average 
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theme leneweron fit Over the Sphere taken to be 5 x 107 


neutrons /em“-sec. The number of neutron captures by hydro- 
gen was calculated by the following equation, which is also 
the number of gamma rays emitted. A point source of 

gamma rayS waS assumed at the core center and attenuation 
of gamma rays through the water and lead shielding was 


calculated as before. 


I! 


neutron captures /em’-sec NO 9 = =O 


neutron captures/sec Va9 = gamma rays/sec 


CO atoms yt 


cm 


N 


I] 


Ore 0 


© 


I! 


V 119,000 ecm 


3(H) = 0.0221 cm™* 
gamma rays/sec = 1.3 x i0rf (assumed to be 2 MeV) 
The results of the dose rate calculation are given in the 


following table: 


Tapes. oe. 


Material Thickness D B 





Water 12 re’. ee 24.0 
Lead i Oigiea. Dyer. roe e 
Totals mer 55 20.2 
ee eS 


G(R) (no attenuation) = 8.9 x 10 yis/om@-sec 


G(R) (attenuated) 0.076 y's/em@-sec 


It 


Dose Rate 


II 


ae me 107+ mren/hr 





jess) 


oc 0. wee” 

In the phase II design the coolant water with its 
yt& activity will not circulate in the upper 4 ft. of 
i@esnecked core Cank. Thus, to calculate the reactivity 
Teeter eacler FOop for the phase TI design, the point 
Pwmree Surenetn calculated in Section 3.3.2.3 will be 
assumed to be located 4 ft. below the surface of the 


water in the necked core tank. Table 3.7.3.3.1 presents 


the calculated results: 


Va enews. 3.5.1 


Material Thickness D Ee. 
Water ie Ila 2.5 
ead ne@ein, 5.6 10 


Totals ZO r 22> 





5.45 x 10° y'g /om@-sec 


5,8 v's /om@-sec 


G(R} (no attenuation) 
O(R) 


Bonversion factor 


(attenuated) 


158 


Hl 


Dose Rate 0.037 mrem/hr 


il 


3.7.4 Neutrons 
Giewatterniac ren equation @iven im Section 3.3.4.1 can 
merwecea to snow thac the fast neutron flux from the core 


meenecarly zero at the reactor top, The equation is: 





eens 
S ry a 
C= ca neutrons /cm*-sec 
{7 

(H,0) = 12 ft. 

(11,0) = 9.98 em. 
& (Pb) = 10 in. 

d,(Pb) = 8.94 om. 


1 
aes oo ¢ neutrons/sec 
~4 wee 

@=5.7 x 10 ~° neutrons/cm’-sec 
fie ,oose Gave cue Go faSt mMeucrons on the reactor top is 
Tous. only 6.6.x ome mrem/hr. 

une Weems race due to thermal neutrons is easily scen 
to be negligible. ‘The diffusion equation for thermal neu- 
fagems tronsaveintiinase plane source in an infinite homo- 


geneous medium is: 


1Ac ¥/T, 


Co 
A aH neutrons /cm ec 


L ami ucLom Lenetun, cm 


Di —Ormmaiiatan coertTrficient; cm 

A= thickness of attenuating material, ecm. 
For H,0, L = 2.88 em, D = 0.183 em, and the column height 
meno0) cm. The thermal neutron flux is thus attenuated by 
about a factor of 2 in the water column above the core. 
3.8 Shielding Cale:zlations, Reactor Shutdown (Phase IT} 


The calculations performed in Section 3.4 were reprated 


Perecne phase 1f reactor design in order to complete the 
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comparison of shielding calculations for the two design 
Piiees her tne Case Of the reactor top Shield plug 
removed, the calculated gamma flux at the coolant tank 
surface after a 14 hour decay period is 7.5 x Tee gammas / 
ei -sec, which converts to a dose rate of 22 mrem/hr. 
Were one fuel element to be raised to a position 8 ft. 
beneath the surface of the water (such as during a fuel 
flipping operation), however, the dose rate at the surface 
of the water would be about 480 mrem/hr. Thus, it is 
apparent that personnel should be protected by shielding 
mem Te OrOucting reiueling or fuel flipping operations. 

the dose rate shortly after shutdown if coolant is 
fest trom the necked core tank is calculated to be le 
rem/nr. After i4 hours the dose rate would decay to about 
700 mrem/hr. 


3.9 Shielding Calculations Summary (Phase II} 


ooo. elooce Through heacpor Top suteld Piug 


Source Dose Rate (mrem/hr) 


Core Fission Gammas 0.364 

Pee! (n,7) al-° 0,485 

uw (n,y) He 0.0002 

nie decay O03 7 

Fast neutrons 0.0006 
Thermal neutrons Ge beet e: 


Tega | 0.89 mrem/hr 
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ee ee oe nOugayneecLvor Top Concrete 
Same as Phase I 0.0018 mrem/hr 

beoecr ote lding Conclusions 

The maximum dose rate anywhere on the reactor top 
is calculated to be 0.89 mrem/hnr. This figure is less 
than the design goal of 1.0 mrem/hr, and thus the 10 in. 
lead shield on top of the necked core tank iS considered 


to provide adequate shielding protection. 
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Chapter 4 
HEAT TRANSFER, PRELIMINARY CALCULATIONS 

4.1 General 

lee was ble pUMDOSe Of Chewprctaminary heat transfer 
Ccmrculavions to determine whether or not the heat from 
tie proposed High Flux MITR core could be adequately re- 
mevea Utrliainge the currenury installed process system. 
Peeemoueh ene Calculations of this chapter were originally 
done for the phase I core design, they are equally appli- 
cable to the phase II design. Section 4.6 will indicate 
the differences in the two phases with respect to heat 
imonlo ley. 

in establish: —~ the limiting temperatnre parameters 
for the operation of the High Flux MITR, the maximum 
Hearen values given in the technical specifications for 
fie presenvynmin were adhered to, These criteria are 
that the temperature at the hottest spot in the hottest 
fuel plate not exceed 450°C (842°F) under any circum- 
stances, and that under normal operating conditions boil- 
mae musc noc be allowed to occur anywhere in th» core. 
ine 450°C temperature limitation was establishel by con- 
emeeration of the structural stability of the aluminum 
fuel elements. While aluminum does not melt un:il it 
reaches approximately 650°C, tbemeevancwco solve sieni fi- 


0 ’ ve ,. 
cantly at about 450°C (18). The requirement that no 
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boiling be allowed to occur in the core was originally 
intended as a safety factor which would act as a backup 
to insure ohnav the fuel plate hot spot would not exceed 
450°C, and also to avoid the problem of flow instability 
caused by boiling. Neglecting the reactivity effects of 
voiding, if the coolant Saturation temperature were ex- 
ceeded sufficiently to cause nucleate boiling at some 
Pernt on the fuel plate, the vapor bubbles formed at the 
iwet piece Suriace Would Cither collapse at the wall or 
would leave the wall and collapse after traveling a short 
(te caiec 8nto the highly suls@¢ooled bulk fluid. The bubble 
formation would cause considerable agitation of the cool- 
ant and thereby lead to a Substantial increase in the 
convective heat transfer coefficient. Heat would thereby 
be conducted away from the fuel plate surface at a greater 
rate, and thus tend to reduce the wall temperature. 
mrevere=as Drei ly mentioned in the introduction to this 
Miescis, Ene reactivity effect Of voiding in the area of 
ee OUvermiost Tucl plates could be posi@rve, aivhough all 
core physics calculations done thus far indicate a nega- 
Give value. Because the possibility of a local positive 
voiding coefficient does exist, however, one is provided 
(ewe ceG INCcenuive co insure boiling does not occur. 
Because the fuel plates are so thin (0.050 in.) 


and because aluminum is a good conductor of heat 





ey, 


(K = 119 BTU/hr-ft-°F for pure aluminum), the maximum 
TonperstUhewc men y DOlLIG Whuninm &@ fuel plate will be at 
most on the order of 3 to 4 degrees F above the fuel 
plate wall temperature at the same point. 
4,2 Distribution of Reactor Heat Load 

The High Flux MITR will generate 5 megawatts of 
thermal ae which must be removed by the reactor cool- 
“hoeovetems, FOr the purposes Of fuel plate hot spot 
Calculations, however, only that portion of the generated 
heat which is transferred to the coolant as it flows up 
tnreuezn the core is of interest. ThiS includes heat 
transferred to the DO coolant both by conduction through 
the faces of the fuel plates and by gamma ray and neutron 
ireeracuLon with the coolant. 

Wolak has determined the total fission energy absorbed 
in the present MITR, operating at 1 megawatt, to be 198,22 
MeV/fission (19). The specific distribution of this absor- 
bed energy is given in Table 4.2.1. 

Although the spectrum from the p-oposed high flux 
MITR will be harder than that for the present MITR, and 
although the amount of aluminum structure may be somewhet 
different, the value of 198.22 MeV/fission for the total 
Mees lLOn enermy absorbed is probably a good estimate f-ir the 


ime Flux MiTR As Shown in Table 4.2.1, not all of this 





ieee «2.2. 1 


Fission Energy Absorption Distribution (Present 





MI'TR) 

Fuel Elements MeV/fission  __% Total 
K,E. Fission Fragments 168,00 
ite beta. Particles eee 
Gamma Energy 4,29 

179.29 90.5 
D,9 
K.E. Fast Neutrons 52 OO 
Gamma Energy 10, 24 

15.24 7.7 
Gam.a Energy Teo 1.8 

Total 198,22 100.0 


energy iS absorbed by the D539 coolant as it passes through 
the fuel elements. 1.8% of the total fission energy, in 
the form of gamma rey captures in the thermal shield, is 
carried away by the thermal shield coolant. The remaining 
98.2% of the fission energy is transferred to the D0 
coolant by conduction through the walls of the fuel plates 
or deposited directly into the D,0. Of the 15.24 Mev/ 
fission deposited directly into the D540 PuulligeLie Teaccor, 
HeVoto has estimated for the present MITR operating with a 


19 element core that 9.03 MeV/fission of the gamma ray 
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energy and 2.5 MeV/fission of the energy, carried by the 
fess Meucronc=are deposited in the D0 outside the fuel 
element boxes (20). Since the ratio of D,0 outside the 
eoolant channels to D50 Piece scooranu channels will be 
evox Mave lL y= one same 10 une High Flux MITR as in the 
present MITR, it seems feasible to assume that the above 
figures will still be reasonably valid. Following 
DeVoto'ts method of estimating the fraction of gamma 
energy deposited in the D540 inside the fuel elements as 
Wiemretctoro1 = -O'Volume Inside the fuel elements to total 


C 
D50 oe one obtains for the High Flux MITR 


(10)* emp 


SE = 0.1644. Therefore, the percentage of total 





(24 \e 
fissicn panes Saers,; deposited in the D5 as 1t flows 
through the fuel elements becomes (0.164)(10.24) + eet. ze 
~ 198,22 


2.1%, Due to the fact that the close packed core will 
iead tO greater particle enerry dissipation before the 
Pros tavo GA Cacao je Ute core than currently exists, the 
above figures for gamma energy deposited in the Ds QO out- 
Side the coolant channels may be a little high and there- 
meeeiOoy COlcety tive. Howcver, the Gifiierence is probably 
Quite small. Hence, it will be assumed that 92.6% of the 
Beta! tission energy will be deposited in the D50 coolant 
as it passes up throusn the fuel elements, and that 90.5% 
Preece TOlLal Tieeion energy will be conducted through the 


walls of the fuel plates and into the D,0 coolant. (Note: 
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approximately 11 MeV per fission is carried off by neu- 
trinos resulting from beta decay. However, since this 
neutrino energy is not deposited anywhere in the reactor 
it is not considered in thermal power calculations. ) 

In the calculations which follow, the heat flux 
mercetrie=tnie Coolant, Lempcravure rise will be reduced to 
92.6% of the total, and the specific heat flux at the hot 
spot will be reduced to 90.5% of its computed value based 
on 5 megawatts. 

4.3 Hot Channel Factors 

POmeCONServative d@sme, 1c ts me@eessary to anticipate 
Me@eerbpclivicam in OpeCratipemcondseions, fwel element fabri- 
Cation, design relations, etc. These uncertainties are 
Spressed in terms of the hot channel factors which are 
elaboraved below. An overall hot channel factor has been 
aematea Or the presens MLTR by Devoto to be 1.25. This 
estimate, however, not onlv includes inaccuracies intro- 
@uced by the manufacturing tolerances, but also includes a 
memo Ceravion Of variations in flux distribution through 
Cach element contained in the core. DeVoto's value is 
somewhat. lower than those normally found in the literature 
due to the fact that he has given credit for several 
eaepects in the design which tend to reduce the value of 
the heat flux at the hot spot. (20) . It is because 


P00 S Fesuics Contain Gonsideracvion of the neutron flux 





dg 


ee tomo Le eprOvooea high Flux MITR differs 
Pe ae Gone eCeDreccniaMl Tt aein this resp@ct, that 
@irect application of DeVoto's results cannot be assumed 
to be valid for the High Flux MITR. 

in determining values of hot channel factors to use 
Pore ule nigh Flux MITR in order Oo eee POP evo 6 1OnS 
in primary design variables caused by manufacturing 
tolerances, it seems appropriate to compare the magnitude 
ac lomo erane coulichearese pool teaasior tiespresent 
MITR, and which will also apply to the High Flux MITR, 
ieee tnose specified forse Brookhaven High Plux Beam 
Imeactor (HFBR) (21). In this magmer the hot charnel 
factors determined for the Brookhaven reactor can be adjus- 
fegescO a5 vO make them applicable to the High Fiux MITR. 
mee contribution tothe overall hot echnanne! factor due to 
the variation in neutron flux along an element will not 
MmewOctcCiInined as a seuarate Suo—laccor, but rather the 
feet. om i meutron flux will be used @Qirectly to deter- 
mene averase power produced per element and hot spot 
Mmocation on the cleiuient, 

oince the calculations to be performed in this chap- 
Wer Will Drimarily we for the hovctest spot in the core, 


aA 


ime POL channel factors of interest wil! be the teat flux 
mo wieay transiter coefiieient hot channel factor:. The 


enous Peuenon Ol Cac i individual item to either o: these 





Te 

Ser NOmmeiomine. taccCOrs are termed sub-factors. The 
Sub-factors are multiplied by Gach other to determine the 
Sveratr OC Matnve! teaccor, A listing and description of 
Lnewhnot charewel sub-factors follows: 

1. Channel Dimensional Tolerance 

PNemeot Spou*tenperavurceritse due to the channel 
eemenslLOnal tolerance wil occur im a channel of minimum 
average dimensions where the local dimension iS a maximum. 
For the Brookhaven High Flux Beam Reactor (HFBR), the 
outermost coolant channels are 0.189 in. thick with 
maximum tolerance of + 0.005 in. The sub-factor deter- 
Peer  1Or tinese GCharacLteristice was 1,09, For the phase 
one Hich Flux MITR, the coolant channels are 0.117 in. 
thick with a proposed maximum tolerance of + 0,005 i 


The sub-factor determined for these characteristics is 











or Snax non) . 221) = 1,11) 
: Oran Gin ioe 


eis Lacu0r 18 applied to the heat transfer coefficient 
Mot channel factor. 


2. Flow Distribution Among the Elements 











Cons19¢ ab20n Of this item was not reported in ref- 
erence (21). For the Iligh Wlux MITR, however, the maxi- 
mum variation in the flow distribution amons the elements 
is estimated to be + 10%, as it is for the present MITR 
(2). It is assuncd that this variation will increase the 


Dee eer COCTTICiCmt nol ch@ne] factor by 1.10. 
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Couetecore jiloy Variagion 

DOmmeene HiGgea wel OCcalevariations in the fuel loading 
are limited to + 5% by the specification of the fuel 
oj meno GhbehnesseateO7O20 +.0,001 in... For the pre- 
sent MITR, however, no specification for the final thick- 
ness of the fuel alloy section is made. Rather, it is 
stated that the overall plate thickness shall not vary by 
here than + 0.003 in., while the variation in the clad 
thickness on either side of the core allcy must be within 
ieee oom eit hs also Stated.2n the MIT specifications 
that there is to be no objectionable segregation of uran- 
Meiiieatoe core alloys Each plate is Subject’ to radio- 
eeaphnic inspection to determine whether such segregation 
Peotovoe. Hence, 1% 1S doubtitul thatethe fuel alloy thick= 
Pesos Olt the spresenv MITER plates exceeds the specification 
as quoted for the HFBR. If the U-235 segregation was 
Suificient to exceed this specification, the fuel plate 
would probably be rejected. 

It is stated in reference (21) that the thickness 
seecifvcaticn made for the HBR can be metewhthout great 
difficulty by the fuel fabricator, except at the ends of 
the plates. Here the fuel alloy develops a local thick- 
ened region due to Gite pdatwe lining process which is called 


mye) Gopbone', Since the svecification for the fuel alloy 
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thickness (excluding the dogbone) as made for the HFBR 
appears sufficiently easy to Satisfy, the Same variation 
will be allowed for the High Flux MI'TR. Although the — 
dogbone thickness in the HFBR is required to be no more 
than 10% greater than nominal thickness, a variation of 
25% was allowed for in their calculations due to (1) the 
fact that the dogbone thickness is a sensitive item in 
the hot spot analysis, and (2) the fact that this is a 
difficult item for the fabricator to control. <A 25% 
variation in the dogbone thickness will also be assumed 
Here toe High NlUx9MEiR. “Thaeewillieeivesan added touch of 
Ses ervauicm vO the not spot calculations performed for 
this design. (In the nhnase II design, the "dogbone" is a 
very sensitive item because the fuel plate hot spot occurs 
in the region of the dogbone at the bottcem of the core. 
werroraved densitometer measurements will probably be 
required to determine the degree of dogboning.) 

inewmie alloy = unrciness 1s Combined with the dimen- 
Prone Saale re.ces Qa,-width and length of the fuel plates 
(and hence heat transfer area). The allowable variation 
in the core heat transfer area for the HFBR is 4.6%. The 


med. vO, eppryOmje Hirsh Flux MIR. 


GQ 


Sane variation is as 
Beer Gme. 2 boc vhe Hiek and the High Flux MITR, the 


heat flux cerm is increased by 1.10 (1.05 x 1.046) 





Cveryvweelre, ex °eoOG ab tne ends of the fuel where it is 
fiieresaseusny) |, clmils25 x 12046), 

4, Reactor Power Measurement 

It is estimated for the HFBR that the power calcula- 
tor will be capable of meaSuring the reactor power to 
within 5%. The same estimate is made for the High Flux 
MITR. Hence, the heat flux hot channel factor is increa- 
Seg 1.05 tovaccount for this estimated variation. 

Demeuest frameeer Cocificient Deviation 

The Phillins Reactor Safeguard Committee recommended 
the use of the modified Colburn correlation in the IIFBR 
with a maximum negative deviation of 20% (21). This same 
correlation is used to determine the heat transfer coef- 
ficient for the High Flux MITR (Section 4.6.3). Conseauently, 
this deviation results in increasing the heat transfer 
Mee icici term by 1.20 for both the HFBR and the High 
Pe MITR. 


6. Summary of the Hot Channel Fectors 





eeu Flux factors, By 


SL OM a OS 


Reactor power measurement error co) 
Fuel allovweartation Le) ( ae) 


Total, F, 1.155 (1.375)* 


*Value at the "deebone" 





irs 


He wiedeeecr ers COoctficicnc. Factors, a 








Channel dimensional tolerance es 
Flow distribution among the elements eaalO: 
Correlation equation ie 

Total, FL 1.446 


Leen ie orewenannelmtaeccor was noe calcu lated. 
instem@d the heat fimeenot chennel factors will We used 
when numerically integrating the power along a channel 
jer bac parpose of calculating weulk temperature rise. 
Tne design of the phase one coolant entrance plenum is 
suc that the flomof coolant up through the outer cool- 
ant channels should be somewhat greater than the flow 
PprOme@meune inner chetaels. ~Conseaucntly, it is argued 
that this effect would overshadow any bulk hot channel 
Sub-factors related to flow velocity, thereby making it 
conservative to simply neglect such sub-factors. It is 
these coolant flow velocity related sub-faectors which 
Peoulnsuish the bulk hot channel factor from the heat 
Mex Oe eniel Taccor. In eddition, by emoloying the heat 
fiwx Nopechannel factors in determinins the bulk tempera- 
reese ticwe = 7571 DULION co tne Dulk coolant tempera- 
ture due to the "dosbone" effect at the fuel plate ends 
fool De taken invo gccount. Tois eiiect woule be missed 


mero meoulk Not «hahnel factor calculavted and used instead, 
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4,4 Power Distribution in the Core 

In order to determine the power distribution in the 
proposed high flux MITR core, the modified twenty grand 
diffusion code was used, and a two dimensional (R,Z) 
plot of the core power resulted. The basic core configu- 
ration and parameters utilized are described here. 

A cylindrical homogenized core 24 centimeters in 
radius and 65 centimeters high was assumed. The core 
and surrounding reflector area were divided into 14 
regions as Shown in Figure eel. The material composi- 
tion of each region is listed in Table 4.4.1. The calcu- 


ijations were done in four neutron energy groups: 


Group 1: co ~ 821 x 10° ev 
Group 2: 821 x 10° - 5.53 x 10° eV 
Group 3: 5.53 z 10° Ol O2a ey 
Group 4: 0.625 - 0.025 eV 


Moe CDiLRermaleeroup Cross Sections for Boron, Hydrogen, 
Oxygen, and Uranium (235) were obtained from the table 


of "Few Group Cross Sections for Metal-H.,0 Systems Fueled 


2 
with u°°>" found on page 142 of ANL-5800. The thermal 
cross sections for the above three elements were taken 
from the 16 group table of "Multigroup Constants for Fast 
Reactor Analysis (LASL)" found on pages 569-574 of ANL- 


5800. The 4 group cross sections for Deuterium and 
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FIGURE 4.4.1. 


LAYOUT OF REGIONS FOR TWENTY 
GRAND CALCULATION. 
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Uranium (238) were obtained by collapsing the 16 group 
cross sections to four energy groups. Fuel plate loading 
was taken to be 1.5 times the fuel loading in the present 
MITR core, or 15 grams of 93% enriched uranium per fuel 
plate. 

It should be noted at this point that the purpose of 
the computer calculation was to obtain a rough estimate 
of the power distribution in the core in order to have some 
foundation upon which to base preliminary heat transfer 
calculations, The assumptions made concerning core size, 
composition, and control rod Peet tion could only be de- 
scribed as educated guesses based on information available 
at the time. Due to the expense involved in running the 
twenty grand program, only one run could be made. More 
detailed and precise core physics calculations were con- 
currently being performed by Andrews Addae. The core power 
distribution is shown plotted in Figure 4.4.2. As expected, 
this plot shows that the outermost fuel plate generates the 
greatest amount of power and therefore it is on the outer- 
most fuel plate that the COre shOUrs DOln Wie soccer. OUNnd, 
Table 4.4.2 lists the power generated (in watts) at the 
various mesh points in the core. (Thermal neutron flux in 
both the core and reflector has been plotted in Figure 
4.4.3 in order to illustrate why the core power distribu- 


tion takes the shape it does.) 
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FIGURE 4.4.2 
CORE POWER DISTRIBUTION (watts/cr’) 
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FIGURE 4.4.3 
FLUX MAP 
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Numerically integrating the power produced by a fuel 
element located on the outermost edge of the core over | 
the height of the core yields a total power of watts ree 
duced in one of these "hot" plates. , 

In order to determine the hottest spot on the hottest 
plate in the High Flux MITR core, the axial power variation 
along the hot channel was investigated. In most reactor 
designs with upward flowing coolant the hot spot, or loca- 
tion of the maximum plate wall temperature, usually occurs 
a few centimeters above the point of maximum power genera- 
m@pon, This is true, for instance, for reactors in which 
mre axial flux distribution is a cosine function which 
Heaks at the f@eometric center of the core. The maximum 
wall temperature occurs beyond the power peak primarily 
due to the fact that the coolant, coming in from the 
bottom at low temperatures, helps keep the bottom portion 
of the fuel element cooler than the top. Such considera- 
tions can also be applied when locating the hot spot in 
the High Flux MITR, and based upon these considerations 
we can expect to find the hot spot about 25 centimeters 
above the bottom of the core for phase one calculations. 
4,5 Wall Temperature Along Hottest Channel 

4.5.1 General 

In arriving at values for the wall temperature 


along the hottest channel, it is assumed that the reactor 





LA5 
is operating at a state in which both the scram condi- 
tions for high coolant outlet temperature and for low 
coolant flow rate have been reached simultaneously. 

4.5.2 Formulation of Applicable Equation 

To determine the magnitude of the wall temperature 
at various points along the hottest channel, one may 
proceed as follows, The temperature difference between 
the maximum wall temperature and the reactor inlet is 
made up of two parts. First is the film temperature 
| drop, which is made up of the differences between the 
wall temperature and the average temperature of the cool- 
ambeac chat point in the channel, Tw - Tc. Second is the 
bulk coolant temperature drop, which is made up of the 
difference between the coolant temperature and the reactor 
inlet temperature, Tc - Ti. Denoting q,/A as the specific 
heat flux at the point of interest and do as the total 
rate of heat transferred to the coolant from the channel 
entrance up to the point of interest, these temperature 


differences may be written as: 
t 


Tw - Tc = q,/hA (4.5.2.1) 


Tc - Ti = 4,/W(ch)Cp(ch) (4.5.2.2) 


The symbols not already defined are: 
h 


convective heat transfer coefficient 


W(ch) = mass flow rate of the coolant in the channel 
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Cp(ch). = specific heat evaluated at the average 


bulk temperature of the D0 in the 


Ti + To | 
2 | 
Adding Equations 4.51 and 4.52, one obtains an expres- 


channel, 


Sion for the maximum wall temperature in terms of reactor 
inlet temperature, the specific heat flux at the point in 
question, and the total heat gained by the coolant from 


the channel entrance up to the point of calculation: 
Tw = Ti + q,/hA + d5/W(ch)Cp + AT, , (4.5.2.3) 


(AT, , is an allowance for error in temperature measure- 
ment, taken to be 5°F.) Following the arguments of 
DeVoto, it iS more desirable that wall temperature calcu- 
lations be based on reactor outlet temperature, which is 
recognized as an incompletely mixed bulk temperature, than 
on reactor inlet temperature as above, This is because 
the outlet temperature recorder is the first monitor en- 
countered by the coolant after it leaves the core and 
hence is the logical place in the circuit to install the 
various alarm and scram circuits. 

By considering the total heat removed by the coolant 
as it passes through the core of the High Flux MITR, an 
expression relating the reactor inlet temperature and the 
incompletely mixed reactor outlet temperature, To, may be 


obtained. Denoting q, as the total heat rate in the 





147 
reactor core, the expression is: 


a = Wy c 


p(r) (TO-T4) ears) 


The symbol W, refers to the total coolant mass flow rate 


ic 
in the reactor core, and Cp (r) is the specific heat of 


Ti + To 


the D,O evaluated at the reactor bulk temperature ——;—— : 


2 
Solving this expression for Ti, and substituting the 


resulting expression into Equation (4.5.3) yields the 


Result s 
°F, q q 
u 2 g 
hA W(ch)Cp(ch) Wear) tm 


(4.5.2.5) 
4.5.3 Heat Transfer Coefficient 
The heat transfer coefficient, h, for turbulent flow 
in confined rectangular passages is calculated by means 
of the modified Colburn correlation (see references (22) 
and (23) for experimental verification of the use of the 


Colburn equation): 


a _ ono ant ome 
= 0.023 [Pete] : eg" : (4.5.3.1) 
where 

h = heat transfer coefficient, BTU/hr-ft°OF 


De = equivalent diameter, 0.0188 ft. ae 


va 
ll 


thermal conductivity of the coolant, BTU/hr -ft°R, 


Fig. (4.5.3.1) 
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V = coolant velocity, 11530 ft/hr for flow rate of 
1850 gpm. 
p = coolant density, lb/hr-ft, Fig. (4.5.3.2) 
lb = coolant viscosity, lb/ft, Baee (4.5.3.1) 
Cp = specific heat, BTU/1b-°F, Fig. (4.5.3.2) 
The subscript "f" means that the properties are evaluated 
at the film temperature. The film temperature is taken 


to be the average of the heat transfer surface temperature 
_ Tw + Tb 

ee aes = 
‘previously defined, and Tb is the bulk temperature of the 


and the bulk water temperature, Tp Tw has been 


reactor coolant at the point of calculation. Equation 


(4.5.1.1) can be rewritten as: 


0.023K, [W,De 0.8 fig }o3 
h = 7 ie jae , Ke], (Ub. 3.2) 





Implicit in the above expression is the assumption that 
the coolant flow rate is the same for all channels, i.e. 
W,/A, = Wch/Ach. Hence, the Reynolds number will be the 
same whether defined for one interplate fuel element or 
mae entire core. Heat transfer coefficient, h, is shown 
plotted versus film temperature in Figure (4.5.3.3) for a 
constant velocity through the core of 11530 ft/hr, which 
corresponds to a flow rate of 1850 gallons per minute at 


a bulk water temperature of 1O4°F, Taking a horizontal 
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cut through the core, coolant flow area was taken to be 
the entire core cross sectional area times the ratio of 
coolant channel cross sectional area to fuel plate Suoee 
sectional area. Taken as such coolant flow area is a. 
conservative factor because flow area would be somewhat 
smaller than this in the actual core, thereby increasing 
flow velocity and consequently increasing heat transfer 
coefficient for the actual core. 

4.5.4 Limiting Parameters 

4.5.4.1 Reactor Outlet Temperature 

Following the reasoning of Section III.3 of the 
Technical Specifications for the present MITR, the maxi- 
mum outlet temperature, To, is taken to be 55°C. (The 
technical specifications quote a maximum outlet tempera- 
ture of 58°C which is the audit and review level. The 
Scram point was set at 55°c.) 

4.5.4.2 Maximum Allowable Wall Temperature 

As stated in Section 3.1, in order to insure that 
the hottest spot on the hottest fuel plate does not reach 
450°C under any circumstances, the criterion established 
is that cooling Shali be sufficient to prohibit boiling 
of any type within the coolant channels. Therefore, the 
difference between the plate wall temperature at the hot 
spot and the saturation temperature of the coolant at the 


Same location, TW. - Dee must not be permitted to 


Ax 
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reach the point of inception of nucleate boiling. 
Several correlations are available which are designed to 
predict the amount of superheat necessary for the incep- 
tion of nucleate (24,25) boiling. (See heat transfer 
texts by L.S. Tong or W.M. Rohsenow.) Three of these 
correlations are used in Section Hach on However, such 
correlations are by no means exact and so for the purposes 
of determining the maximum allowable wall temperature it 
is prudent that we attempt to design our system such that 
_we simply rely on the fact that boiling will not take 
place as long as the fuel plate wall temperature does not 
exceed the coolant saturation temperature at any point in 
the core. The top of the phase one design core is under 
Bert. of D0; at this pressure the saturation temperature 
of the 1,0 is 223.5°F; 
4.5.4.3 Coolant Flow Rate 

Following arguments of Section III.2 of the techni- 
cal specifications for the present MITR, the minimum 
allowable flow rate through the core of the High Flux 
MITR will be taken as 1850 gallons per minute. 

4.5.5 Calculations 

Utilizing Equation 4.5.2.5, here repeated, wall 
temperature is calculated as a function of axial height 
above the core bottom for the hottest (outermost) fuel 


plate. 


154 


eee, + q,/hA + Go/Wonln ch - q./w,C + AT 


Den m 


(4.5.5.1) 


The first calculation will be done in detail as a sample. 


For a point on the outermost fuel plate 10 centimeters 


above the 


wo 


ol 


bottom of the core: 

55°c = 131°F, a limiting condition by definition 
FyGqyP(2=10)V 

heat flux hot channel factor = 1.155 
fraction of core power conducted through the 
fuel plates into the coolant = 0.905 

power produced per unit core volume = 113.1 
watts /em? (from Table 4.5.5.1) 

volume of 1/2 fuel plate plus 1/2 coolant 
channel per unit axial dimension = 1.60 em? 
(see Figure 4.5.5.1) 

heat transfer area for one Side of one plate per 
unit axial length = 7.61 em* (see Figure 4.5.5.1) 
1265 BIU/hr-ft°-°F from Figure 4.5.3.5 for a 


film temperature of 168°F; F, = 1.46 from Sec- 


h 
tion 4.3.6 
10 
eG = P(Z)Vv(Z) 
2 do Zia) 


fraction of core power carried off by the cool- 


ant as it passes through the core = 0,926 





FIGURE 4.5.51 e 
CORE VOLUME ELEMENT 





| COOLANT 
CHANNEL 
0.050" 0.117" 





P(Z) = 


V(Z) = 


ch 


p,ch 


C = 
p,r 


Substituting 


q,/nA = 


ISIS 


power produced per unit core volume at 
axial position Z in the hot channel; anes 
taken from Table 4.5.5.1 
volume element of core producing power | 
P(Z) = V x mesh spacing 

one-half the coolant mass flow rate through 


a Vachanne ly) _ 
one coolant channel = v/aw, Uenannsi) . 


965 lb/hr 


= average specific heat in the hot channel up 


to the point of calculation = 0.9995 from 
Figure 4.5.3.4 

total power produced by the reactor which is 
renovea Dy Ghewcoolant aS a6 passes up 
through the core = (0.926) (5x10°) =a C2 ax 


10° 6 


watts = 14.5 x 10~- BTU/hr 

total mass flow rate of coolant through the 
core = 1,020,000 lb/hr (1850 gpm at 104°F) 
average specific heat in the core = 1.0005 
for a core bulk temperature of 122g 


in the numerical values and solving: 


(1.46) (1.155) (0.905) (3.412 BIU/hr-watt) (113.1 watts /cm?) 5 


(1265 BTU/hr-ft~-°F) (7.61 cm“) 


(1.60 cm?)(30.48 em/ft)* = 91°F 
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Io/Wony, ch ~ 


(1.375) (0.926) (3.412 BTU/hr-watt)(1015.0 watts/cem~) (1.60 em“) 
(965 lb/hr)(0.9995 BTU/1b-°F) . 


= 7.3°F 
q,/W,C _ (0.9237 5x10° watts )(3.412 BTU/hr-watt 
ee (1.0005 BIU/1b-°F)(1.02x10~ 1b/hr) 
= 15.5°F. 
Hence, 


ioe] + 78S = 15.5 + 5 amenocear , 


Repeating the above procedure for the other values 
of axial distance above the bottom of the core resulted in 
the values plotted on Figure 4.5.5.2. As shown on this 
figure, the maximum fuel plate wall temperature is 231-F 
and occurs about 23 centimeters above the base of the 
core, At this position in the core, the saturation tempera- 
ture of the D,0 coolant is 225,4°F, Although the above 
calculated hot spot value is slightly above the coolant 
Saturation temperature at that point, the amount of sSuper- 
heat is insufficient to cause boiling by calculation of 
any of the three correlations for the inception of nucleate 
boiling given in Tong, Boiling Heat Transfer and Two Phase 
Flow (25). 

For the calculated values of heat flux, pressure, and 


Saturation temperature at the position of the hot spot, the 
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wall temperature necessary to cause nucleate boiling to 
commence is calculated as follows: 


Correlation (1), 


2.3/_0.0234 
! 15.60 pr 156 p 


Ge eat 


q" = 124837.8 BTU/hnr-ft* 


O 
me = 225,68 
p = 19.17 psi 
= O 
Ty =9eo9, 2h 


Correlation (2), 


= ~_ alt 
Trp = Meat + ST yu, - a/b 
a _ P/900 1! 6 eet 
— oO 
To. = 254 4p 


Correlation (3), 


= 0.015 x 10° (q"/10°)1/" 


6 apDY/ att 


He OC = Cea ial ae )s10 © 
262.0°R 


by 
ll 


5 : 

Consequently, even if one were to abide by the re- 
quirement of Section 4.5.3.2 that Tamrac not exceed T,.,. 
only minor modifications would be required to the fluid 
flow system. Since the main coolant pumps are capable of 


supplying approximately 2100 gpm of D0 through the core 


one could, for example, Simply redefine the Scram point 





iol) 


flow rate to be 2000 gpm. In addition, one might re- 
define the outlet temperature scram point to be 53°C. 
The combined effect of redefining scram point flow rate 


and outlet temperature would be to reduce T by 


w,max 


approximately 8°F, more than enough to reduce T 
w,max 


below T Preliminary calculations, therefore, have 


sat* 
indicated that fluid flow and heat transfer would not 
pose a major design problem. 

4.5.6 Revised Data 

As Andrews Addae began to obtain reliable results 
from his phaSe I core physics calculations, one fact 
became very clear: the maximum core power density would 
be much higher than the values used in the eAleuvations 
of the preceeding section. Best estimates indicated a 
maximum homogenized core power density of approximately 
320 watts/cm°, a factor of 2.5 greater than the maximum 
power density computed in Section 3.4. It should be 
noted that the computer code used by Mr. Addae in his 
calculations is based on diffusion theory, and as Such 
values of power density near the edge of the core where 
the power density gradient is steep are only approximate. 
Since the film temperature rise tern, q,/nA, is the 
dominant term in the calculation of fuel plate wall 


temperature (Equation 4.5.2.5), it was clear that to a 
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first approximation one is faced with the necessity of 
increasing the magnitude of the product "hA" by a factor 
of approximately 2.5 for the Phase I design. 

4.6 Phase I vs Phase II 

With respect to heat transfer, the major difference 
between the phase I and phase II core configurations is 
that the phase II configuration is much smaller. The 
power distribution in the two cores is very similar. For 
a constant mass flow rate, the flow velocity through the 
phase II core is a factor of 2.1 greater than that through 
the phase I core, which increases the heat transfer 
coefficient by a factor of 1.81. Partially offsetting 
this advantage is the fact that the maximum power density 
in the phase II core is a factor of dls Drenter than that 
for the phase I core. The combined effect is that for the 
light water cooled core, instead of the factor of Ae 5 
augmentation in heat transfer required by Section 4.5.6, 


Ten eas 


1.8. Maximum power density for the phase II core occurs 


one must augment heat transfer by a factor of 


along the bottom edge of the core. 

In addition to being smaller, the column of water 
above the phase II core is 3 ft. higher than that above 
the phase I core, The additional head of water increases 
the coolant saturation temperature by about 5°R, thus mak- 
ing it conservative to use estimates based on phase I cal- 


culations to develop a phase II heat removal analysis. 
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° Chapter 5 
AUGMENTATION OF FUEL PLATE HEAT REMOVAL 


| 


Results of core physics calculations being performed by 


5.1 General 


Mr. A. Addae (21), together with the results of the prelimin- 
ary heat transfer calculations performed in the previous 
chapter of this thesis, have indicated that the product "ha" 
must be increased by a factor of approximately 2.5 (Phase I) 
or 1.8 (Phase II) in order to provide adequate heat removal 
from the High Flux MITR core. It is the purpose of the 
paragraphs which follow to explore the various means by 
which the product "hA" may be increased, and to experimen- 
tally verify any conclusions reached. 
5.2 Methods of Increasing the Product "ha" 

Pach DeCreased Muel Plate spacing 

The fuel plates in the present MITR are spaced So as to 
allow a coolant channel width of 0.117 inches. A Survey of 
research reactors utilizing flat plate MITR-type fuel ele- 
ments has revealed no parameter which would prevent plate 
Spacing from being decreased. The narrowest flow channel 
discovered was the 0,050 inch channel between the fuel 
plates of the Oak Ridge High Flux Isotope Reactor (22). 
For a rectangular channel with high aspect ratio, decreas- 
ing the fuel plate spacing causes a proportional decrease 
in the effective channel diameter and increase in flow 


velocity. Taking the McAdams correlation for heat transfer 
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coefficient as typical of the several available, 


0.8 0.4 
) =e), 3 ae os) 


hD pVD 
e e 
y— = 0.023( 
b Hy 





it is seen that for a constant mass flow rate the decrease 
in equivalent channel diameter and corresponding increase 
in flow velocity cancel each other, thus leaving the 
Reynolds number unchanged. AS a result the Nusselt number 
remains unchanged, which in turn requires that the magni- 
tude of the heat transfer coefficient h increase in direct 
proportion to the decrease in flow channel equivalent 
diameter. (Equivalent diameter is calculated for rectangu- 
lar channels using the standard equation De = 


4 X cross sectional area 
wetted perimeter 


decrease the coolant channel width from the present 0.117 


-) Therefore, if one were to 


inch to 0.090 inch, one could expect the heat transfer 
coefficient to be increased by a factor of 1.3. (Note: A 
0.090" coolant channel width was assumed for all Phase II 
calculations. ) 

5.2.2 Increased Mass Flow Rate 

Increasing the mass flow rate of coolant through the 
core results in a proportional increase in flow velocity 
through the core. Referring to Eq. 5.2.1.1 it is seen 
that the heat transfer coefficient is proportional to 


velocity to the 0.8 power. If one were to add a third main 
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coolant pump, thereby increasing the flow velocity by 


approximately 50%, one could expect the heat transfer 
coefficient to be increased by a factor of 1.38. 

When an increase in mass flow rate is contemplated, 
one must also investigate the associated increase in pres- 
sure loss around the primary coolant circuit. 


AP = 2f 5 > ft. of H,0 (Ss2K2m) 


As shown by Eq. 5.2.2.1, pressure loss varies with the 
Square of the flow velocity. Consequently, for a 50% 
increase in flow velocity one could expect pressure losses 
to increase by a factor of 2.25. In the present MITR the 
pressure loss through the core is on the order of 0.3 feet 
of H,0 (2). From Eq. 5.2.2.1, pressure loss through the 
8-inch piping is calculated to be approximately 0.07 feet 
oy HAO per foot of pipe. Since each main coolant pump 

has a capacity for providing a pressure head of 46 psi 
while pumping 1000 gpm, it is seen that neither the pres- 
sure loss through the core nor the pressure loss through 
the 8-inch piping would prevent a 50% increase in flow 
velocity. However, the pressure loss through the heat ex- 
changers is approximately 15 psi. Consequently, if one 
were to add an additional main coolant pump to increase 
flow velocity, an additional heat exchanger may be required 


to limit the pressure loss across the heat exchangers, 





166 


In summary, the extensive equipment modifications and 
additions necessary to obtain a 50% increase in flow 
Vevloctamaweuld beva highwpricemtospay for a factor of 
1.38 increase in heat transfer coefficient. 

5.2.3 Modifications to Fuel Plates 

The various methods currently available to augment 
convective heat transfer have been summarized in a paper 
by Dr. A.E. Bergles (27). Dr. Bergles has grouped the 
augmentive techniques available into the categories of 
surface promoters, displaced promoters, vortex flow, 
heated surface vibration, fluid vibration, electrostatic 
fields, and additives. Of these categories, all except 
surface promoters and vortex flow can be immediately dis- 
carded as not being feasible due to the nature of the MITR 
core heat transfer problem. The idea of using some type 
of twisted tape insert in order to generate swirl flow 
(vortex flow) in the coolant channel was seriously con- 
Sidered. However, the concept was finally discarded due 
to serious doubt as to the effectiveness of a swirl 
generator in a narrow rectangular channel of high aspect 
ratio, and due to the possible danger of a tape breaking 
loose and blocking the coolant flow through a channel. The 
only applicable method by which the product "hA" can be 
increased is, therefore, by modification of the fuel plate 


surface, 
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Unfortunately, reference (1) compares the various 
means of augmenting heat transfer as found in the litera- 
ture with the corresponding unaugmented case based on main- 
taining constant pumping power. For the case of coolant 
flowing through the MITR core, pressure loss through the 
core appears to be of minor concern. Consequently, a 
comparison of heat transfer data based on constant mass 
flow rate is desired. The results of a literature search 
for data reporting the effect of various types of surface 
promoters is presented in the following section. 

5.3 Literature Search 

5.3.1 General 

Table 5.3.1.1 was taken primarily from reference (1) 
and lists the various experimental investigations of arti- 
ficial roughness. One fact stands out clearly; there is 
no data for water cooled thin rectangular plates. Results 
of investigations which most nearly approximate the situ- 
ation in the MITR core are presented below. No attempt at 
completeness is made here. The original papers should be 
consulted if detailed information is desired. 

5.3.2 Brouillette, Mifflin, and Myers 

The authors! purpose was to determine the effect of 
a type of internal fin on the heat transmission and pres- 
sure drop characteristics of water flowing through copper 


tubing. The fins were produced by cutting 60-degree, 
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ROHL e  oterOeteet 


EXPERIMENTAL INVESTIGATIONS OF ARTIFICIAL ROUGHNESS - 


NONBOILING HEAT TRANSFER AND FRICTION 


Investigators 


Cope 
Sams 


Kolar 


' Lancet 


Dipprey and Sabersky 


Nunner 


Koch 


Nagaoka and Watanabe 


Sams 


Kreith and Margolis 


Boelter et al. 


Edwards and Sheriff 


Kemeny and Cyphers 


Description 


Cooling of water; internally 
knurled tubes 


Heating of air at high AT; square 
threads inside tubes 


Heating of water and air; internally 
threaded tubes 


Heating of air; rectangular channel 
with machined roughness 


Heating of water; Sand-grain type 
roughness in tubes 


Heating of air; rings of various 
cross section inserted inside tubes 


Heating of air; axially supported 
rings inserted inside tube 


Heating of water; wire coils inside 
tubes 


Heating of air; wire coils inside 
tubes 


Heating of water and air; tubes 
with wire-coil inserts 


Heating of air; rectangular duct, 
one surface heated, with transverse 
strips 


Heating of air; rectangular channel, 
one surface heated lined with wire 


Heating of water; annulus, inner 
heated surface with spiral grooves 
or protrusions 
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aD Aes Saco lee 


(Continued) 


Investigators 


Brauer 

Bennett and Kearsey 
Brouilllette, Mitilin, 
and Myers 


bairant et al, 


Wilkie 


Sheriff and Gumley 
Kjellstrom and Larsson 


Draycott and Lawther 


Kattchee and Mackewicz 


Sutherland 


Description 


Heating of water; annulus; inner 
heated surface with grooves 


Heating of Superheated steam; annu- 
lus, inner tube heated with grooves 


Heating of water; tubes with inter- 
nal V-shaped fins 


Heating of water; annulus, inner 
heated tube with diamond knurls 


Heating of air; annulus, inner tube 
with Square or rectangular protru- 
sions 


Heating of air; annulus, inner heated 
tube with coiled wires 


Heating of air; annulus, inner heated 
tube with various protrusions 


Heating of air; annuli and tube 
bundles with threads, knurls, and 
coiled wires 


Heating of nitrogen (water for fric- 
tion tests); wire coils outside rods 
in a bundle 


Heating of air; wire coils outside 
rods in a bundle 





lifes 


Nunner correlated his data for Reynolds numbers ranging 


3 


from 10~ to 10° with the result: 


| 
Nu = £/8 RePr ; (Sea. oe) 
ae f 
1 + —y7g—yelPr ee - 3) 
Re Pi O 
: AP d 
where f is the Darcy friction factor (f, =2—- ) 
2 pv 
for the rough surface and f is the smooth tube friction 
factor. Equivalent diameter for all rough tube calcula- 
tions is taken to be the volumetric mean diameter, d= 
fi 
Th * 

Nunner's data indicated that the larger the €/d 
ratio the greater was the increase in heat transfer coeffici- 
ent. His data also indicated that the optimum value of 
the ratio €/L, where € is the height of the ring insert and 
L is the center-to-center spacing of the insert in the tube, 
is approximately 10. However, since the manually installed 
ring inserts in Nunner's experiment may not have been per- 
fectly seated against the tube wall, the effect of increas- 
ing the number of ring inserts per inch of tube length may 
have been to decrease (or at least fail to increase) the 
effective area available for heat transfer. Consequently, 
it is not possible to predict the optimum fin spacing dis- 
tance based on the results of Nunner's work, Nunner's 


work did indicate, however, that ring inserts were capable 
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V-shaped notches on the inner surface of five specimens of 
tubing. The number of notches or fins was varied from 8 

| 
to 32 per inch of length and the ratio of fin height to 


tube diameter ranged from 0.009 to 0.050 inch. Data was 


taken for Reynolds numbers ranging from 1.5 x 10+ to 
ex 10°. Table 5.3.2.1 below summarizes the results obtained, 
In addition, the authors stated that fin height has a 


greater effect in increasing the heat transfer coefficient 


than the number of fins, 


Tape. 5235.2... 


Tube € hnA( fin A( fin h( fin 
No. ss hA(smooth) _A( smooth) h( smooth 
1 0.009 eee lace 0.89 
2 0.013 1.08 Heple 0.91 
3 0,017 WS) 1.48 1.04 
4 0.027 1.90 1.44 Le 
D 0.050 Bt 1.48 1.68 


Sec. 3 SNunner 

Nunner investigated the effect of various sizes, 
Shapes, and densities of ring inserts on heat transfer 
coefficient and friction factor for air flowing through a 
50 mm diameter brass tube. Both round and rectangular 
springy split rings were manually inserted into the tube 


at equal distances apart to serve as roughness elements, 
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of increasing the Nusselt number for air flowing ina 
Svindrical Quceu by Neally a Lacuor Of Chree. 

5.3.4 Kolar 

Kolar analyzed the mechanism of heat transfer for 
both air and water flowing through smooth and rough tubes. 
Internal tube roughness was formed a machining a 60° 
triangular thread into the tube. The smooth tube was 26 
mm in diameter and the rough tube thread depth varied from 
0.5 to 1.5 mm. Kolar's smooth tube data fit the standard 


. Dittus-Boelter correlation; 


0.85.04 


Nu = 0.023 Re (So alias) 


In order to correlate his rough tube data Kolar proposed 
a relation in the same form as the Dittus-Boelter equation 


for smooth surfaces: 


Nu = 0.0517 Re,Pr°*? , (5.3.4.2) 
where Re» is the Reynolds number defined by the friction 
velocity, 

Up = U, Vivo, 


wherein Un is the mean velocity and Ur is the friction 

velocity. Kolar shows that his data is in general agree- 
ment with that of Nunner and also with that of Cope (30). 
(Cope's data was for water flowing in circular tubes with 


the inside surface roughened by knurling) Reynolds number 





ye 

was varied from 4500 to 145,000 and Prandtl number from 
0.71 to 5.52. Unfortunately, Kolar makes no mention of 
the fin effect and what influence this effect might have 
on the rate of heat transfer. 

D=3.> ances 

Lancet investigated the effect of a knurled surface on 
the heat transfer coefficient for air flowing through a 
small rectangular duct. The rough tube was 4.5 inches long 
and the maximum width and height were 0.145 in. and 0.045 
in., respectively. The diamond-shaped protrusions were 
mepmeominavlclyO,O10 an. x O201O0 an. x 0.010 in. and were 
Spaced approximately 0.010 in. apart leaving an open duct 
area of 0.025 in. x 0.125 in. For a constant Reynolds 
number in the range 3110 <¢ Re< 26900, Lancet found that the 
Nusselt number was increased by about a factor of two, 
Lancet's data is fit quite well by Nunner's relation for 
Nusselt number, Eq. 2.3.1, but even better by Nu = 0.042 
Re? 8 pri/3 | The fin effect is not mentioned in Lancet's 
results. 

Lancet's results for a Smooth rectangular channel were 


found to fit the Colburn equation for smooth round ducts: 


Nu = 0.023 Re?:Spy1/3 


; (Sack oad: 
5.3.6 Conclusions 
From the literature survey, it is apparent that the 


increase in heat transfer from a surface roughened duct is 
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proportional toechne increase in the duct roughness ratio, 
€/d. When correlating cylindrical tube data, the equiva- 
lent diameter, d, is taken to be a volumetric mean saee Lee 
defined as d = ae When attempting to apply the eee 
of experiments performed on cylindrical tubes in order to 
predict the effect of roughness on rectangular ducts of 
high aspect ratio, it is not clear how equivalent diameter 
should be defined. (Aspect ratio is the ratio of width to 
height of a rectangular channel.) When applying smooth 


tube correlations such as those of McAdams or Colburn to 


rectangular ducts equivalent diameter is taken to be d = 
4 x cross sectional area 


LA 
mETSCer perincter (=>). Synan, however, has shown 


that the rough rectangular channel data of Lancet is best 
correlated with the data of other investigators by using 
the volumetric mean diameter, d = (ae as the equivalent 
diameter (35). It should be noted, though, that the 
aspect ratio for Lancet's rectangular ducts was only 5: l. 
For a MITR coolant channel, the aspect ratio is approxi- 
mately 28 : 1, a considerable difference. 


For the proposed MITR coolant channel shown in Figure 


>. bee sks 
P i i d 
dn = TI, — 0.496 Tk, 5 — 0.04 
| a 
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Since it has been concluded that the increase in heat trans- 
fer is proportional to the roughness ratio, €/d, utiliza- 
tion of the correct value of d is imperative if one ae to 
predict rectangular duct heat transfer coefficients based 
on cylindrical tube experiments. From reference (35) one 
would be tempted to make the rectangular-cylindrical 
comparison based on d = V/7L . However, using the 
volumetric equivalent diameter implies that the increase in 
heat transfer due to roughness will vary with the width of 
the rectangular channel. For example, for a given roughness, 
€, a 3-inch wide thin rectangular channel would have a 
smaller €/d ratio than a 2 1/2 inch wide rectangular chan- 
nel, and thus a lower Nusselt number, This is obviously 
not correct. It would appear, therefore, that for an arti- 
ficially roughened thin rectangular channel with high aspect 
ratio, the correct equivalent diameter lies Somewhere 
between 4A/P and fe . Using the pessimistic, oo 
equivalent diameter, the heat transfer coefficient for the 
channel of Figure 5.4.2.1 can be estimated (using the same 


velocity and temperature conditions as in Section 4.5.3): 


Area 
Increase 
d C App ny hy Due aK 
in. d Dittus-Boelter Kolar Nunner to fins DB 


0.496 0.04 1360 2520 2270 leroy 1.83 





176 


On the basis of these estimated values it can be concluded 
that the finned channel of Figure 5.4.2.1 will provide at 
least a factor of 1.5 increase in heat transfer over an 
equivalent smooth channel. The heat transfer augmentation 
may, however, be greater than a factor of 1.5, possibly as 
much as a factor of 3 or 4. | 
5.4 Experimental Program 

5.4.1 General 

Based on the literature survey of the previous section, 
it was apparent that heat transfer from artificially 
roughened cylindrical tubes is at least partially understood 
by correlation of roughness ratios, €/d. However, it was 
not clear how one should apply the results of rough tube 
experiments to predict rough rectangular channel heat trans- 
fer. An experiment which would duplicate as nearly as possi- 
ble the conditions which will exist in the High Flux MITR 
core was clearly desirable. 

5.4.2 Coolant Channels 

After making a brief survey of the methods available 
to heat the walls of a test section coolant channel, it 
was decided to heat the coolant channel by passing an elec- 
trical current through the channel walls. As described in 
Section 5.4.3.2 the heat transfer laboratory power supply 
is capable of supplying 3000 amperes at 24 volts. 
Consequently, if the test section is to take full advantage 


of the power supply available, it should provide a resistance . 





iL 


of 8 x 107° ohms. Since, 
_ Pl 
R= £ eee 
where R = resistance (ohms) 
L = channel length (cm), and 
A = cross sectional area for current flow (cem“), 


p = resistivity (ohm-cm) 
there are three parameters which can be adjusted in order 
to obtain a correct resistance value. A stainless steel 
- channel wall 2.5 inches wide by 0.020 inches thick was 
decided upon. For stainless steel of these dimensions, a 
coolant channel wall electrical resistance of 6.8 x 107° 
ohms was calculated. Maximum power input into the coolant 


channel was therefore estimated to be 61,200 watts, which 


will provide a maximum heat flux into the coolant flowing 
BTU 


hr-ft 
In order to provide a basis for future comparison, 


through the channel of 2.5 x 10° 


the first channel to be constructed and tested was a 

smooth surface rectangular duct. A flat piece of 0.020 
inch thick stainless steel was folded three times to form 

a rectangular channel with inside dimensions of 2.5" x 
0.090". The Aneta closed by a Single weld which 

ran the entire length of the coolant channel. It was 
expected that test results obtained from this channel would 


agree fairly well with existing smooth tube correlations, 


- 
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such as those of Colburn or McAdams. 

The second channel constructed was designed to provide 
a "proof of principle" for internally roughened rectangular 
ducts. It was decided to mill fins into the stainless 
Steel channel wall in preference to other means of wall 
roughening (such as knurling) for the rough wall run 
because, (1) the increase in heat transfer due to heat 
transfer surface area versus that due to increased turbu- 
lence could be accurately separated, and (2) the process 
was one which was within the immediate capability of shop 
equipment and personnel. 

Fin height was chosen to be 0.020 inches based on the 
philosophy that in the High Flux MITR no point of two fuel 
plates should be closer than the smallest existing separa- 
tion of reactor fuel plates, the 0.050 inch separation of 
the Oak Ridge High Flux Isotope Reactor fuel plates. The 
projection of opposing 0.020 inch fins into a coolant 
channel whose unfinned width was 0.090 inches results in 
a minimum spacing of 0.050 inches. Fin width was taken 
to be 0.020 inches in order to provide the fins with 
structural stability and to prevent the machining problems 
which would accompany the milling of very narrow fins. 

Fin spacing waS a compromise, and in fact, somewhat arbi- 
trary. From Nunner's work (Section 5.3.3) the optimum fin 


Spacing to fin height ratio appears to be 10 : 1, which 





aie? 
would require a fin spacing of 0.200 inches for a 0.020 
inch fin. However, from the work of Brouillette, Mifflin, 
and Myers it appears as though the optimum arrangement is 
to have fins spaced as close together as possible. It is 
this author's opinion that the manually inserted rings of 
Nunner's experiment are not capable of providing data on 
which to base conclusions concerning optimum fin spacing. 
It is felt that heat transfer will be increased by spacing 
fins close together due to the greater heat transfer area 
increase over widely spaced fins. However, there must 
exist some minimum spacing beyond which fluid circulation 
between fins will be hindered and boiling in the inter- 
fin spaces will be caused to occur. With the above infor- 
mation in mind, a fin spacing of 0.040 inches was chosen, 
The optimization of fin spacing might well be the object 
of future experiments. 

Rectangular grooves 0.040 inches wide, 0.020 inches 
deep, and spaced 0.020 inches apart were milled cross-wise 
three at a time into a flat piece of stainless steel. Two 
grooves 0,100 inch nate were then milled lengthwise into 
the stainless steel plate at positions which were to 
become the channel edges. The stainless steel plate was 
next folded into a rectangular channel and welded as with 


the smooth channel, 
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5.4.3 Description of Apparatus 
5.4.3.1 Hydraulic System 


The experimental facility used was the low-pressure 
test loop located in the MIT heat transfer eta ae 
The basic apparatus was designed and constructed in 1961. 
A schematic Seine loop is presented in Figure 5.4.3.1.1. 
The pipes and fittings are made of brass and stainless 
steel for corrosion nesetiece. Rayon reinforced rubber 
hose is used where flexible connections are required. 
Distilled water is circulated by a bronze, two-stage, 
regenerative pump capable of providing a discharge pres- 
sure of 260 psig at 3.6 gpm. The pump is driven through a 
flexible coupling by a 3-hp Allis-Chalmers induction motor. 
A Fulflo filter is installed at the pump inlet. Pressure 
fluctuations at the outlet of the pump are damped out by 
means of a 2.5 gallon Greer accumulator charged with 
Piiero2en, The accumulator contains a flexible bladder- 
type separator which prevents the nitrogen from being absor- 
bed by the system water. After the accumulator, the flow 
splits into a bypass line and a test-section line. 

In the test-section line, fluid flows through a 
Fischer-Porter flowrator followed by a preheater, thence 
through a Hoke metering valve and the test section, after 
which it merges with fluid from the bypass line. The flow 
then goes through the heat exchanger and returns to the 


pump. The preheater consists of four Chromalox heaters of 
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approximately 6 Kw each. Three of the heaters are con- 
trolled simply with "on-off" switches while the fourth 
can provide a continuous range of O to 6 Kw by means of 
a bank of two variacs mounted on the test bench. Quick- 
action Jamesbury ball valves are installed before the in- 
let to the flowrator and after the exit from the test 
section. The ball valves permit quick isolation of the 
test section to minimize fluid loss when conducting burn- 
out tests. The exit valve can also be used to adjust the 
test Section pressure. 

Flow through the bypass line is controlled by a ball 
valve on each side of which there is a 300 psig pressure 
gage. Pump operating pressure, and hence the pressure up- 
stream of the test section, iS controlled by this valve. 

The heat exchanger is a counterflow type with system 
water flowing in the inner tube and city water in the 
outer annulus. The heat exchanger iS capable of maintain- 
ing a constant pump inlet temperature. A Fulflo filter 
is installed on the city water line to reduce scale forma- 
tion in the exchanger. 

The distilled water was de-ionized continuously by 
passing a portion of the flow through four mixed-bed 
resin demineralizer cartridges installed in a Barnstead- 
"Bantam" demineralixer unit. A 4.7 gallon degassing tank 
is provided with five electrical heaters (3-200 Vac and 


2-110 Vac). The degassing tank also serves as a surge tank. 
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A 15-gallon stainless steel storage tank for filling the 
System is mounted directly above the degassing tank and 
can be filled with distilled water from standard ac 


bottles with a small Hypro pump. | 


5.4.3.2 Power Supply 
Test section electrical power is provided by two 
36-KW dc generators connected in series. Each generator 
is rated at l2 volts and S,0CO amperes. The power control 
console permits coarse or fine control from O to 24 volts. 
A water-cooled shunt installed in parallel with the test 
section protects the generators against the shock of the 
sudden open circuit which would occur at burnout. Power 
is transmitted from the main bus to the test section by 
water-cooled power leads. In the present tests, the down- 
Stream or exit power lead was at ground potential. Rubber 
hose connected both inlet and exit chamber plenums to the 
main loop in order to electrically isolate the test sec- 
tion, 
5.4.3.3 Instrumentation 
All temperatures were meaSured by copper-constantan 
thermocouples made from 30-gage Leeds and Northrup duplex 
wire. The test section inlet and outlet temperatures were 
measured by thermocouples directly in the fluid stream. 


The thermocouples were introduced into the piping 
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immediately before and after the test Section through 
Conax Lrtcines. 

Coolant channel inlet and exit presSures were measured 
with Helicoid 8-1/2 inch gages with full scale readings of 
200 psig and 100 psig, respectively. Both are specified 
to an accuracy of + 0.25% of full scale. Pressure loss 
across the test section was alSo measured by a manometer 
located adjacent to the test stand. 

A float tube calibrated for flows of 0 to 8 gallons 
per minute was installed in the basic Fischer-Porter 
flowmeter housing and provided measurement of the test 
section flow. 

The voltage drop across the test Section was read 
directly on a Digatec multiple range dc voltmeter with a 
specified accuracy of + 1/2 percent. The current flow was 
determined by using a Wilson-Maeulen multi-volt meter 
(having 0-30 and 0-60 millivolt scales) to measure the 
voltage drop across a calibrated shunt (60.17 amp/mV) in 
Series with the test section. Thermocouple voltages were 
read from a HONS wel Brow continuous reading chart 
ieec Order. | 
5.4.3.4 Test Section 

The test section consisted of a Stainless steel cool- 
ant channel (described in Section 5.4.2) which protruded 


into two stainless steel end blocks. The end blocks 
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served,to act aS an entrance plenum for the coolant water, 
to provide a point of attachment for the electrodes of the 
power supply, and to serve as a point of attachment for 

the pressure taps. Each end block was fabricated in three 
pieces; a base piece into which a 1/2 inch NPT was machined 
in order to provide a point of attachment for the coolant 
hose, and two entrance plenum sections which were machined 
so as to taper Peete. channel from a 1/2 inch diameter 
cylinder to a 0.130 inch by 2.5 inch rectangle. Figure 
5.4.3.4.1 shows both a horizontal and a vertical cross sec- 
tion through an end block, Figures 5.4.3.4.2 and 5.4.3.4.3 
Show the test section in various stages of construction and 
assembly. The rectangular coolant channels protruded into 
a tight-fitting slot in the end blocks to a distance of 

3/4 inch. Once inserted, the coolant channel and end 
blocks were bonded together with a stainless steel solder, 
and pressure taps were drilled into the coolant channel at 
the positions indicated in Figure 5.4.3.4.1. During the 
finned channel experimental runs, it was discovered that 
water was seeping around behind the coolant channel inside 
the end blocks, Consequently, the pressure taps, which 
protruded into the end blocks, were reading the pressures 
of the inlet and exit plenums rather than the pressures 
inside the inlet and exit of the coolant channel itself, 


as intended. In order to eliminate the inlet and exit 
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A. Three-piece end block before B. Channel end block showing 
assembly. Coolant channel in coolant channel slot. 
foreground. 





wa} 





C. Finned coolant channel inter- D. Assembled coolant channel 
ior shown before bending and and end blocks. 
welding plate into rectangular 


channel, 


FIGURE 5.4.3.4.2 
EXPERIMENTAL TEST SECTION 
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A. Coolant channel inserted into B. Test section assembled. 
Support block. 





Meerest Section in place in test Ds Instrumentation stand, 
Poe lity. 


FIGURE 5.4.3.4.3 
EXPERIMENTAL TEST SECTION 
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effects and thus to measure only the pressure drop across 
the coolant channel, pressure sensing tubes were inserted 
through the endmost thermocouple holes in the front 
fiberglass block and through the coolant channel wall at 
those points. The existing pressure tap holes in the end 
blocks were plugged. | 

The coolant channel was physically supported by, and 
was electrically and thermally insulated by two 2" x 6" 
Nema grade G-5 fiberglass blocks. A rectangular channel 
2 1/2" x 0.125" was milled into one of the fiberglass 
blocks in order to accept the coolant channel. The two 
fiberglass blocks were held together by 16 1/8" bolts 
Spaced along either side of the coolant channel. Holes 
were drilled in the fiberglass blocks as shown in Figure 
5.4.3.4.1 in order to provide access for the thermocouples 
to the coolant channel. 

Teflon tape was used to electrically insulate the 
coolant channel from the thermocouples. A total of seven- 
teen thermocouples were mounted on the coolant channel 
(see Figure 5.4.3.4.4 for arrangement). The thermocouples 
were pressed into place and held onto the coolant channel 
by the fiberglass support blocks. Four thermocouples were 
mounted near the surface of the fiberglass blocks in order 
to monitor the temperature gradient across these insulating/ 


Support blocks. in order to maintain a near zero temperature 
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THERMOCOUPLE ARRANGEMENT ON 
THE’ SURFACE OF THE COOLANT CHANNELS 
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gradierit, three separate flexible tape heaters were wrapped 
around the fiberglass blocks, Heat output from these 
three heaters was adjusted using three variacs to maintain 
a nearly zero temperature difference between the outside 
channel wall and the fiberglass block surface, thus creating 
an adiabatic wall condition at the channel outer wall. The 
fiberglass blocks and heater strips were blanketed with 
fiberglass insulation in order to reduce the required 
heater output and to cause the strip heaters to give a more 
| uniform heating effect. 

Four steel Support arms were bolted between the fiber- 
glass blocks and the stainless steel end blocks in order 
to join the entire test section more firmly. By so doing, 
the probability that the solder joints between the coolant 
channel and end blocks would crack during handling was 
reduced. 

The method used to change coolant channels in the 
test section was to heat the end blocks (not the coolant 
channel) sufficiently to melt the solder permitting removal 
of the channel. Once removed and allowed to cool, the end 
blocks were emissed in a crock of nitric acid. The nitric 
acid successfully removed all the remaining solder from 
the end blocks without any damage to the stainless steel, 
thus providing a clean smooth surface for installation of 


the new channel. 
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5.4.4 Operating Procedure 

With all system components installed and ready for 
operation, and with both the test section inlet and exit 
valves closed, the coolant pump was started and the bypass 
valve adjusted to provide 20 psig backpressure to the test 
section. Next the test section exit valve followed by the 
test section inlet valves were opened and adjusted to pro- 
vide the desired flow through the test section. All high 
point bleeder valves were opened and trapped gasses 
allowed to escape. At this point, the guard heaters were 
turned on and adjusted to provide a insullating block sur- 
face temperature equivalent to the expected coolant channel 
wall temperature. The bus bar cooling water was turned on 
as well as the city water to the heat exchanger. The 
final stage in the startup process was accomplished by 
Synchronizing and starting the generators. After allowing 
the generators to warm up, they were brought up to power 
and allowed to stabilize. The flow of city water to the 
heat exchanger was adjusted until a constant inlet coolant 
water temperature to the test section was obtained. With 
all conditions stable, data was taken. This data included 
coolant channel inlet and exit pressure, voltage drop 
across and amperage flowing through the test section, and 


all thermocouple readings. 
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valves closed, the coolant pump was started and the bypass 
valve adjusted to provide 20 psig backpressure to the test 
section. Next the test section exit valve followed by the 
test section inlet valves were opened and adjusted to pro- 
vide the desired flow through the test section. All high 
point bleeder valves were opened and trapped gasses 
allowed to eScape. At this point, the guard heaters were 
turned on and adjusted to provide a insullating block sur- 
face temperature equivalent to the expected coolant channel 
wall temperature. The bus bar cooling water was turned on 
aS well as the city water to the heat exchanger. The 
final stage in the startup process was accomplished by 
synchronizing and starting the generators. After allowing 
the generators to warm up, they were brought up to power 
and allowed to stabilize. The flow of city water to the 
heat exchanger was adjusted until a constant inlet coolant 
water temperature to the test section was obtained. With 
all conditions stable, data was taken. This data included 
coolant channel inlet and exit pressure, voltage drop 
across and amperage flowing through the test section, and 


all thermocouple readings. 
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5,4.5 Assumptions for Data Analysis 


In analyzing coolant channel heat transfer data, the 


Colburn correlation Nu - = 0,023 Re 27° pees 


was 
used as the basis for comparison. From the basic equation 
for heat conduction Q = hA(T,-T,) , Q, T,, and T, are 
experimentally measured, and A is the constructed internal 
Surface area of the coolant channel. The heat transfer 
coefficient can thereby be calculated directly. For pur- 


poses of comparing smooth channel experimental data with 


the Colburn correlation, equivalent diameter was calculated 
eae 

= - 

In order to compare heat transfer from a smooth chan- 


using the relation, D. 


nel with that from a finned channel, it is necessary to 
establish certain constants for comparison. For the pur- 
poses of reactor design, it is both desirable and reasonable 
that both the center-to-center spacing of the fuel meat in 
the core, and the mass flow rate through an individual 
coolant channel remain constant. Consequently, in compar- 
ing the smooth channel with the finned channel, it was 
assumed that the two channels had both the same equivalent 
diameter and the same cross sectional area for coolant flow. 
Thus, for a given mass flow rate Reynolds number (smooth) 
equals Reynolds number (finned). As a result, a relative 
comparison of heat transfer ability was obtained. 

The only temperature correction incorporated into the 


experimental data analysis was the correction for temperature 





7 


re 


drop agross the coolant channel wall. Coolant channel 
outside wall temperature is measured whereas inner wall 
temperature is desired. The correction was made using 


the relation, 


Twi = Two 7 ORR : 
where Q/A = heat flux 
K = thermal conductivity of the channel wall, and 
t = channel wall thickness. 


Although to a first approximation the heat flux along the 
coolant channel was uniform, a small correction was required 
to the heat flux as a function of position along the coolant 
channel due to the temperature dependence of resistivity. 
Since resistivity increases with increasing temperature 
according to the relation p = p,(l+aT), more heat would be 
generated near the exit of the electrically heated coolant 
channel than near the entrance. The temperature coeffici- 
ent of resistivity for a piece of the stainless steel 

used in the fabrication of the coolant channels was measured 
over a range of temperature between that of liquid nitrogen 
and 100°F water and found to be 0.002/°F, Wall temperature 
rise between channel initet and exit was typically on the 
order of 50°F, Therefore, the actual input heat flux to 

the channel varied typically from -5% at the entrance to 


+5% at the exit from the channel average. Since for a given 
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measured temperature drop between the channel wall and 
bulk fluid heat transfer coefficient is directly propor- 
tional to input heat flux, a corresponding correction was 
made to the experimentally measured heat transfer coeffici- 
ents. 

5.4.6 Experimental Results 

Smooth channel data is presented in Figure 5.4.6.1 
and is shown to be about 12% above the Colburn correlation. 
Such increased heat transfer was expected because the 
smooth channel tested was not perfectly smooth. When the 
thermocouples were pressed onto the channel outside wall 
small indentations were formed which protruded through the 
wall and into the coolant flow path. These thermocouple 
"nipples" acted as roughness elements in promoting heat 
transfer. The smooth channel data taken is best fit by 
oe cee . (It should be 


noted that a leak developed at the entrance solder joint 


the correlation Nu p = 0.027 Re 


between the test section end block and coolant channel. 
It is felt that this leak did not affect experimental 
results, however. ) 

Finned channel data is presented in Figures 5.4.6.2, 
5.4.6.3, 5.4.6.4 and 5.4.6.5. The upper curve of Figure 
5.4.6.2 was drawn through points representing the average 
value of calculations from 10 centerline thermocouple 


readings for each experimental run. The curve is fit by 
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0.8 O13 
f Pro ; 


along the lower curve of Figure 5.4.6.2 represent average 


the relation Nu OOo a= Ke The points 


values for each smooth channel experimental run. In | 
calculating the heat transfer coefficients for the finned 
channel experimental runs, it was observed that the heat 
transfer coefficient was highest near the entrance of the 
coolant channel and decreased toward the exit of the 
channel. This effect did not show up in the smooth channel 
runs. Heat transfer coefficient averaged over the seven 
experimental runs is shown plotted versus thermocouple 
position along the coolant channel in Figure 5.4.6.3. 

Since one would normally expect the heat transfer coeffici- 
ent to increase slightly toward the exit of the coolant 
channel, the validity of the finned channel heat transfer 
data was subject to investigation. It was suspected that 
some aSymmetry existed either in the thermocouple readings 
or in the channel itself. In order to investigate this 
possibility the test section was inverted. Everything 
remained as it was except, now, upside-down, and an entire 
new set of experimental runs was made. Average values of 


O53 


Nu/Pr are shown plotted in Figure 5.4.6.4 and are fit 


0.8 5, 0.3 


by the relation Nu, = 0.050 Re p r » almost exactly 


f 
the same as the relation which fit the original data, 

Heat transfer coefficient versus position along the coolant 
channel is plotted in Figure 5.4.6.5. Although the slope 


of this curve is not as great as that of Figure 5.4.6.3, 
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the same trend of decreasing heat transfer coefficient with 
distance along the coolant channel is present, and the 
Shape of the curves is identical. At present there 
appears to be no good explanation for this effect. 
However, it is felt that the average .value curves of Fig- 
ures 5.4.6.2 and 5.4.6.4 can be used with confidence in 
heat transfer calculations for the proposed High Flux MIT 
Reactor. The average values are conservative when one 
considers that due to control rod insertion into the reac- 
tor core, only the lower 12-16 inches of the fuel plates 
will be generating a Significant amount of power, and it 
is in this region near the coolant channel entrance that 
the heat transfer coefficient was measured to be above 
average. 

In computing average values for the finned channel 
data, only the readings from thermocouples along the chan- 
nel centerline were used. Heat transfer near the edge of 
the finned channel was observed to be much greater (by 
about a factor of 2) than heat transfer along the channel 
centerline. The most probable reason for this effect is 
that there existed a small gap between the channel side 
walls and the ends of the transverse fins in which 
secondary (additional) turbulence was generated. 

It is recommended that the modified Colburn correla- 


Halon, Nur —meRO23 Re 0*® Pr,'* » be used for smooth 


eter} 
' ° 
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channel? calculations and that the correlation Nu, =O O50 


O20 0,3 
PY p 


Re be used for calculations for finned channels 


a 
Similar to that of this experiment. Heat transfer from 
the finned channel is thus a factor of 2.2 greater than 
heat transfer from a Smooth channel. - The heat transfer 
surface area of the finned channel was a factor of 1.67 
greater than that of the smooth channel. Since fin 
efficiency in this case was nearly 100%, the factor of 2.2 
increase in heat transfer can be taken as the product of 
a factor of 1.67 increase in heat transfer area multiplied 
by a factor of 1.3 due to increased turbulence caused by 
the fins. To be more precise, the product "hA" was in- 
creased by a factor of 2.2, while "nh" alone was increased 
Dyea factor of 1.3. 

Pressure drop data for both the smooth and finned 
channels is presented in Figure 5.4.6.6. The smooth 
channel pressure drop was much larger than that normally 
associated with a smooth channel, but, like smooth channel 
heat transfer, this can be attributed to channel roughness 
caused by thermocouple indentations. Channel entrance and 
exit pressure losses are included in the plotted value of 
Smooth channel pressure drop. 

Finned channel pressure drop was very large, with AP 
varying roughly proportional to yi 6, This large pressure 


loss is a limiting factor on the coolant mass flow rate 
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which ean be pumped through the core using the existing 
pumps. | 
5.4.7 Analysis of Results 


D,0 pressure measurements around the primary coolant 


2 
loop of the existing MITR show that there exists a 23 

psig pressure drop between the primary heat exchanger 

exits and the core inlet piping (36). The pressure loss 

in the piping between these two points Should be no more 
than 1-2 psi. The 23 psig pressure drop can be explained 
by the fact that between the heat exchanger exits and the 
core inlet are valves which are kept partly closed in 
order to regulate the primary coolant flow, Were these 
valves completely opened, it is reasonable to assume that 
an additional pressure head on the order of 21 psig would 
be available. Consequently, core pressure drop could be 
increased by approximately 21 psig while maintaining a 

flow rate of 2000 gpm with the existing primary coolant 
pumps. However, a pressure drop of 21 psig through a 
finned channel corresponds to a flow rate of only 4 gallons 
per minute. Were 10.95 gpm to be pumped through a smooth 
channel, equivalent heat transfer would be realized. 
Implicit in the above information is the fact that there 
exists a certain cutoff core size, assuming no additional 


pumping power is to be installed. 





eon 


Smooth coolant channels will give better heat in 
cores of smaller size, whereas it is advantageous to Nee 
the finned coolant channels in cores of larger size, 
Table 5.4.7.1 below gives a brief analysis of two coolant 


channel sizes of possible interest. 


Hamden Sy 47a 


2.5"x0.090" 2 EO ROO 
channel channel 


Max. finned channel 
flow rate* 


channel flow rate 10.95 gpm 9.20 gpm 


4 gpm 3.36 gpm 





218 channels, or 
15-15 channel 
elements 


183 channels, or 
11-17 channel 
elements 





Break even point 






Actual heat transfer 


hA( finned 
Increase eS 


*Core AP ax = ct psi 





1.84 for a 26 
element core 


Lec ©Orean lc 
element core 






Thus, for purposes of the MITR redesign, the benefit to be 
derived from fuel plate finning varies with the size of 
the core under consideration. 

One core configuration considered for the phase II 
High Flux MITR contains 26 fuel elements, with each element 
containing 15-2.1" x 0.090" coolant channels. Comparing 
Smooth channels to finned channels for this core configura- 


tion shows that the pressure drop limited finned channel 
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will provide a factor of 1.84 better heat transfer than a 
smooth channel. However, since the heat transfer area of 
the finned channel is a factor of 1.67 greater than the 
smooth channel, the heat transfer coefficient is in fact 
increased by only a factor of 1.1. More important, 
though, is the fact that the factor increase of 1.84 is 
greater than the factor of 1.8 required for the light water 
core. While the optimum plate roughness condition for 
heat transfer augmentation has not yet been achieved, it 
has been demonstrated that the currently installed process 
System is capable of removing the heat from the hottest 
fuel element of the proposed High Flux MITR core. 

For a flow of 2000 gpm through a 26 element core, 
5.15 gpm would flow through each coolant channel. From 
Figure 5.4.6.6 this flow rate would cause a pressure drop 
across the core on the order of 33 psi, or 12 psi above 
the existing system capability. However, if one were to 
either install new main coolant pumps or install booster 
pumps downstream of the heat exchangers capable of provid- 
ing an additional le psi of pressure head, the finned 
channel pressure drop limitation could be eliminated, thus 
allowing one to fully realize the 2.2 factor increase in 
heat transfer for a finned channel over a corresponding 
Smooth channel. In addition, if new coolant pumps were to 


be purchased, pumps of not only higher head capacity but 
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also greater flow rate capacity could be installed, thus 
adding even a little more to the degree of heat transfer 
augmentation possible. 

5.4.8 Further Investigations 

At the suggestion and with the assistance of this 
author, the heat transfer and pressure drop characteris- 
tics of a longitudinally finned coolant channel have been 
investigated by Mr. Donald Uhl. A longitudinally finned 
(fins running parallel to the coolant flow path) fuel 
plate enjoys three advantages over a transverse finned 
(fins perpendicular to the coolant flow path) fuel plate 
such as investigated in the previous sections: (1) the 
longitudinal fins provide little additional turbulence and 
thereby create less pressure drop across the core, (2) 
fabrication of the fins is less time consuming because 
fewer milling cuts are required, and (3) fuel plate buck- 
ling cannot completely block a coolant channel because 
flow can always occur in the grooves between fins. 

It was originally intended that the coolant channel 
plate be milled so as to produce fins 20 mils wide, 20 
mils in height, and spaced 20 mils apart, thus increasing 
surface area over that for a smooth plate by a factor of 
2. However, due to difficulties experienced in milling 
the fins into the stainless steel plate (difficulties 


which should not be encountered when milling fins into 
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the actual aluminum fuel plates), the spacing between 
fins was increased from 20 mils to 40 mils, thus providing 
an area increase of 1.67. 

PresSure drop data versus flow rate for the longitu- 
dinally finned coolant channel is shown on Figure 5.4.6.6. 
Although larger than that for a eth channel, pressure 
drop is shown to be considerably less than that for a 
transverse finned channel and will not impose a limiting 
condition on primary coolant flow rate. 

Experimental heat transfer data is presented in 
Figure 5.4.8.1. Each point shown is an average or the 
data from each of the 10 thermocouple readings recorded 


for each run. The experimental data is fit by the rela- 


0.8 0.3 
f Pr, 5 


finned channel augments heat transfer by a factor of 1.6 


tion Nu, =20.03/ Re Thus the longitudinally 
over that from a smooth plate. Since the area ratio, 
finned to smooth, was 1.67, it is seen that the increase 
in heat transfer was 0.958 times the increase in surface 
area, Thus, it is reasonable to expect that if the 
longitudinal fins were to be spaced 20 mils apart, thereby 
doubling the heat transfer area, one could expect heat 


transfer to be increased by a factor of 1.92. 
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Chapter 6 
CONCLUSIONS 

It is not intended that the work which has been 
presented represent a final design or a final design 
analysis, but rather that the preceeding work make some 
contribution toward that final design. Chapter 2 has 
Shown that while much work remains to be done in the 
area of mechanical design, there seem to be no prohibi- 
tively difficult or unsurmountable design problems. 
While the shield design and supporting calculations of 
Chapter 3 must certainly be refined based on final 
Shielding requirements, it has been shown that shielding 
will not pose a difficult design problem. Chapter 4 has 
indicated the need to find a means of augmenting the con- 
vective heat transfer from the High Flux MITR fuel plates. 
In response, the experimental work of Chapter 5 has demon- 
strated that by modification of the fuel plate surface the 
required heat transfer augmentation can be achieved. 

In surveying the work pertinent to the contents of 
this thesis which must be accomplished prior to setting 
a final High Flux MITR design, it is apparent that mechani- 
cal design work is being thoroughly covered by Mr. E. 
Barnett. Final shield design can be obtained by simple 
extrapolation of the work of Chapter 3. It is in the 


area of core heat removal that the most intriguing area 





a 


for future work lies, that of optimizing the means by 
which fuel plate heat transfer is to be augmented. In 
order to make a tentative decision based on information 
available at this time, the following Summary analysis 
based on heat transfer coefficient, h, is presented for 
the phase II core design. 

For a maximum homogeneous power density of 1.3 x io 
fissions/cm°-sec along the base of the core (equivalent 
to that assumed in previous analyses), the maximum heat 
flux through fuel plates spaced 0.090" apart will be 
2, 1o°x 10° BTU/hr-ft*. At the coolant entrance, the 
Saturation temperature of H,0 is 228.5°F. Assuming the 
reactor were just about to scram due to high coolant out- 
let temperature, bulk inlet coolant temperature would be 
114,.5°F, Thus, since fuel plate wall temperature is the 
sum of bulk coolant temperature plus the film temperature 
rise from the bulk fluid to the fuel plate wall (T, = T, + 
Q/hA), in order to prevent boiling the quantity Q/hA must 
be less than 114°F. Employing hot channel factors and 
solving for heat transfer coefficient yields a maximum 


a 


required value of h = 3750 BTU/hr-ft<-°F, Were hot 


channel factors not employed, the required value of "h" 
would have been 1870 BTU/hr-ft°-°F, From the experimental 


results of Chapter 5 at a flow rate of 2000 gpm, "h" is 
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calculated to be 2430 BTU/hr-ft°-°F for a smooth channel, 


3880 BTU/nhr-ft°-°F for a longitudinally finned channel 
(fins 40 mils apart), and 5300 BIU/hr-ft*-°F for a trans- 


verse finned channel. Repeating in tabular form: 


hyequired (with hot channel factors) 3750 
beeanired (without hot channel factors) 1870 
Asmooth channel 2430 
eee elie (spaced 40 mils apart) 3880 
Dee ences eee (spaced 20 mils apart) 4660 
hi nansverse fing (requires booster pump) 5300 


It is apparent that if hot channel factors are not employ- 
ed even the smooth channel is capable of providing adequate 


heat removal. The ratio 


Dsmooth ehamiied = 

Nrequired Wathout Hac. . - 

However, a safety factor of 1.3 is not considered adequate 
due to the uncertainties involved in this "first-of-a-kind" 
design. 

On the basis of the above analysis, it is recommended 
that longitudinal fins be milled into the fuel plates of 
the High Flux MITR. The fins should be 20 mils in height, 
4O mils in width, and should be spaced 20 mils apart. 
Should milling difficulty require that the fin spacing be 


increased, the spacing may be increased up to but not 


beyond 40 mils. 
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Appendix A 
EXPERIMENTAL DATA 
In the experimental data which follows, the run 
numbers beginning with "A" are the smooth channel runs, 
run numbers beginning with "B" are the initial transverse 
finned channel runs, and run numbers beginning with "c" 
are the transverse finned channel runs after the channel 


0.3 values 


was inverted. Average Re and average Nu/Pr 
are an average of calculations based on centerline 


thermocouple readings only for the finned channel runs. 





Data Sheet 
Run No. A-l 
Flow Rate = 1.66 gpm 


Electrical Power = 43,900 BTU/hr 


i 


eel 


Q/A = 52,800 BTU/hr-ft* 

Inlet Pressure = 6.1 psi . 

Exit Pressure = 4,6 psi a 

Inlet Temperature = 98°F 

Exit Temperature = 151.5°F 

Taele Two Twi tb TP h Re Nu/Pr?: 9 
2 201 197 +145 #172 +«211010 8430 27.9 
QR 193 189 145 167 #+«®\11200 8150 S2n0 

192 188 138 163 £41050 #47920 28.5 

4 190 186 131 £159 960 7760 25.8 
AR 186, 182 131  @i67 ~1e205. “7590 Peo 
5 173. 169 125 147 #«2»41200 7642 S15 
5R 166 162 125 144 =14e0 6920 Sis 
6 170 166 125 146 1285 #4246950 6c2G 
Gi Lfdie  WiOs. M25) Wilelon meson maac® 8046 
7R 169 165 '125 145 1320 #6930 34.5 
8 165 161 =$118 =$139 #1240 ~=#=6650 eral 
8R 168 164 °-118 141 1145 6660 29,8 
9 7 eases ld 860 6680 223 
10 15% 152 105 geo Seiieoe ers 28.5 
10R 166 162 105 134 920 6390 Eel 

Average Re = 7230 

Average Nu/Pr?: > = 29.7 
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Data Sheet 2 


Run no. A-2 


Flow Rate = 2.5 gpm 
Electrical Power = 60,460 BTU/hr 
Q/A = 72,550 BIU/hr-ft* 


Inlet Pressure = 9.0 psi 


MER 2 psi 

Exit Pressure = 7.0 psi . 

Inlet Temperature = 91°F 

Exit Temperature = 139.5°R 

ona Tio Twi ‘» Te h Re Nu/Pr?:? 
2 189 183 134 158 1475 11400 39.4 
2R 176 170 134 152 2010 11000 eee 

17or 173 les 150m 1626 10900 42.8 

4 186 180 122 150 £1240 10900 S270 
AR 176 170 122 146 = 1510 10600 39.4 
5 17 165 15 Mom 1470 10050 SCS, 
5R 166 160 115 138 # £1630 9820 41.9 
6 171 165 115 #4140 ~~ 1480 10050 38.2 
7 169) 168) 115 £l14@m 1526 10000 39.4 
7R 175 169 115 14e 1345 10200 34.7 
8 160 = 154 109 132 1640 9400 41.8 


8R 164 158 109 134 1490 9400 38.1 

9 161 155 103 129 1400 9300 35.4 

10 153147 97 122 1465 8650 So 

10R 163° «157 oO] 127 en 9050 30.6 
Average Re = 10,050 


Average Nu/Pr?>? = 38,9 
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Data Sheet 3 


Run no. A-3 

Flow Rate = 3.4 gpm 
Electrical Power = 53,860 BTU/hr 
Q/A = 64,630 BTU/hr-ft* 


inlet Pressure = 15.2 psi 


AP ele (as 1 
Pate Pressure = 13558psi 
inlet Temperature = 101°F 
Exit Temperature = 133°F 
een Two Twi Tb Tp h Re Nu/Pr?: 3 
2 165 159 129 i144 2170 14700 Doe 
2R 161 156 129 143 2460 #13650 64.0 
157 152 125 139 2430 13200 62.5 
4 a58 9253 21 9137 -Be¢6o e200 5G 
AR M56 (5 t We es eee sooo 58.4 
5 tO a 7 30 2350 ens 200 G2n2 
5R 151 eS G7 “esl esc e540 58.8 
6 150 44 117 #+19131 2400 13200 oO 
7 1445 140 117 #+%1128 2940 #12500 74.7 
7R oa 56 “lay a7 “167.0 ers2c0 coer 
8 142 86137 | 113 125 2740 #11800 69.0 
8R 154 = =148 —2113°—=«131 1840 12500 46.6 
9 142 136 109 123 2410 11670 60.5 
10 235 G30 GOS Melly We2e20 WesCo 65. 1 
10R W448 143 #105 3ehk 1705 #£«£411760 43.0 


Average Re = 12,740 


Average Nu/pr?: 3 = 58,6 
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Data Sheet 4 224 


Run no. A-4 
Flow Rate = 3.5 gpm 
Electrical Power = 73,090 BTU/hr 
Q/A = 87,770 BTU/nr-ft* 
Inlet Pressure = 14 psi 
AP = Ell WS 1. 
Exit Pressure = 11.6 psi 
Inlet Temperature = 91°F 


Exit Temperature = 133.5°R 


Average Re = 13,100 


Average Nu/pr?:* = 61.6 


Souris wo wi 7» tr h Re Nu/pr? 
2 170 163 125 01468) 25a) 15700 655 
2R 165 158 128 143 #+.2920 ~=614400 (Ox0 
3 162 155 123 139 2740 38214000 71.4 
4 162, I56GR 11eme 13 Coeem I seCe 60.6 
AR 15g) 5a 1ISmeet3s9" 278m 13970 (On 
5 158 iGy 112) 120mN 25008 (12650 65a 
5R 154. 147 +1112 +130 2510 #212900 63.6 
6 154 147 #+2.+112 #+|11380 2510 #212900 63.6 
7 156, HO 112) 1339) 2370. 12900 60.6 
7R 162 155 ‘Lie 134. 2040 13250 pe 
8 148 141 #107 +124 2580 #12300 65.0 
8R 157 150 “107 129 2040 #12850 5 Lae 
9 146 139 102 #2121 «2370 =#& 11800 59.2 
LO 143 +136 96 116 2190 11340 54.4 
10R 152 39145 96 121 1790 #11800 AM 7 
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Data Sheet 5 


Run no, A-5 

Flow Rate = 4.0 gpm 

Electrical Power = 62,530 BTU/hr 
Q/A = 75,040 BTU/hr-ft* 


Inlet Pressure = 16.7 psi 


Average Re = 15,080 


Average Nu/pr?:? = 62,1 


AP = 3.0 psi 

pxivy Pressure = 13.7 psi 

Inlet Temperature = 100°F 

Exit Temperature = 131.5°F 

| BOhale Two Twi ‘> Tr h Re Nu/pr?-? 

2 166 160 128 i144 2345 #16640 Gare 
ER 162 156 128 142 #2680 16360 69.7 
3 160 154 124 139 2500 #16080 64.7 
4 159 gi53 ) G20 essence | 15700 5G 
AR 155 149 120 135 #2590 #15460 Gime 
5 150 144 116 130 2680 14870 68.1 
5R 151 145 116 131 4.2590 14870 65.7 
6 150 144 116 130 2730 #14870 69.3 
ih 150 gs 116" “1e0m) (2m) 141870 69.3 
7R 159 153 116 9912555 2020) 5915460 52.0 
8 1446 140 112 #126 2680 14530 67.4 
8R 153 J47@ 112 1309) 2i4o  Mise70 54.5 
9 1446 j14o0 108 124 2345 #14020 DOR 
10 142 #136 06104 €6120 © «6©2340)6=—- 113620 on 
10R 150 144 104 124 1900 #14020 47.9 
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Data Sheet 6 


Run no’, A-6 — 

Flow Rate = 4.2 gpm 

Electrical Power = 55,500 BTU/hr 
Q/A = 66,590 BIU/hr-ft* 
Inlet Pressure = 15.8 psi ; 

AP = 1.8 psi 
Exit Pressure = 14 psi 

Inlet Temperature = 119°F 


Exit Temperature = 146°R 


Average Re 


18,990 


Average Nu/pr?:° = 65.5 


7 couic Two Twi "bp Ts h Re Nu/pr°- ? 
2 177 #172 #4143 «157 2230 @8§=©20500 59.8 
2R 168 163 143 152 #3300 #£#+19700 oven 

172 «©1670 =: 139 s-d153)S—2390 =~) —:19750 Ger 
4 ie “ley sou tea) 2s ae ESGOO 56.4 
AR N63 “W5S 136595175) sono Seseoo (oro 
5 165 160 133 146 2480 19200 64.8 
5R 157 “53 133) Baiee ee seios “Tesco 86,2 
6 165 “60 933 Wiig 2420 9200 64.8 
7 165 “860 133° 146" )2550 ~ 16200 64,8 
7R 169 164: 133 148 #2160 #19000 bone 
8 165 160 129 144 #$%§.2210 19400 et 
8R 163 158 129 143 2360 18300 61.4 
9 162 57 1265 ee eee 136200 56.3 
10 156 “151 “W230 W365 "2550 1700 60.5 
10R 154. 149 123 #136 2490 ~=617850 63.9 
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Data Sheet 7 


Run no. A-7 

Flow Rate = 4,95 gpm 

Electrical Power = 86,460 BTU/hr 
Q/A = 103,750 BTU/hr-ft* 
ivet Pressure = 23.2 psi . 
AP = 3.7 psi 
Exit Pressure = 19.5 psi 

Inlet Temperature = 103°F 


Exit Temperature = 138,5°F 


Average Re = 19,750 


Average Nu/pr?** = 73.2 


Tasic Two Twi Tp Te h Re Nu/pr?*? 
2 180 171 134 $$153 2780 22050 fe 
2R 176 168 134 #151 #£®.3110 £21600 82,2 

173 164 130 147 #22990 21050 ont 
4 173. 164 #125 145 £2650 #20700 69,2 
AR 168 160 125 142 $%§$}3030 #19900 Ton 
5 161 152 121 136 3340 #19300 86.0 
5R 173 164 #121 142 8 2410 #19900 62.6 
6 163 155 121 =%.138 $3060 19450 79.4 
i 162 153 121 #138 #3230 £19400 Ogee 
7R 173, 164 , 121 i14e e410 8619900 62.6 
8 160 152 116 134 ($3000 #19000 75.9 
8R 169 160 £116 138 $=+:.2390 #19500 61.8 
9 160 152 212 23255 2630) 16700 67.1 
10 151 143 #108 $125 #+.2960 38 17800 Wo 
10R 159 150 108 129 2420 8618350 62m) 
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Run no, B- 


il 


Data Sheet 8 


Flow Rate = 1.38 gpm 


Electrical Power = 34,800 BTU/hr 


Q/A = 41,800 BIU/hr-ft* 


Inlet Pressure = 


Exit Press 


wc 


—- 
— 


9.5 psi 
4.5 psi 


Inlet Temperature = 


Exit Temperature 


Thermo- 
couple 


a 
en 


10 
1OR 


Average Re 


Two 
ino 
D 
169 
158 
si) 
139 
Or 
148 
142 
138 
143 
io 
126 


ee 


Average Nu/Pr 


39D 
Ons 


= 152°R 


Twi 
ee 
172 
166 
> 
154 
135 
134 
145 
139 
134 
139 
133 
123 
124 


20 


Ton r 


- 


133 
133 
Hae 
127 
127 
127 
Tey 
120 
114 
107 


107 


Taro 


b 
146 


146 
139 


162 
ie 
153 
144 
143 
131 
130 
iG 


133 


130 
130 
Tas 
it 
116 


1300 
eS 
1570 
1940 
2020 
4920 
Doe 
2300 
3480 
5430 
2200 
2235 
2770 
2560 


5 psi 


Re 
6650 


6500 
6190 
5760 
5710 
5170 
5140 
5400 
5290 
5140 
5140 
4820 
4490 
4530 


228 


Nu/Pr?° 


35.1 
42.8 
41.7 
50.8 
ee 


125.4 
152.0 


5 oe 
89,2 


138.2 


56.0 
Sor 
68.5 
Gime 


S 
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Data Sheet 9 
Runa no, b=<2 
Flow ore = 2,0520m 
Electrical Power = 48,400 BTU/hr 
Q/A = 58,100 BTU/hr-ft~ 


Inlet Pressure = 13.8 psi 


Average Re = 7,640 


Average Nu/Pr?- ? = 67.6 


Exit Pressure = 6.5 psi *_— 

Inlet Temperature = 94,5°F 

Exit Temperature = 143°F 

ie Tyo Twi to | Tp h Re Nu/pr?: 3 
2 io “Sig als 7 sel 1720 9140 45.6 
2R 173 168 137 42.152 41900 8980 Does 

165 160 131 145 2040 8520 53.4 

4 152 147 #421125 136 2620 #47900 Cis 
4R mo HS “le5 “ier 2 Secoo 63.8 
5 i338 eS “Wr | 123 6455 7050 162 
5R 133 128 $119 123 £6455 #£«7050 162 
6 146 141 #119 134 $=$.}2665 #7800 68,4 
i We7 “832 Glo 125) 42805 s72ce 108 
7R 134 129 #119 124 5480 7140 WS 
8 iss 133 “113 23> 2690. 7660 72.5 
9 130 “25 S107 ile “sil6O Seago Tome 
10 1260 115 01) 108s V4eo Giles 97.5 
10R m4 719 “101 1230 3280 £6250 (Sce 





Data Sheet 10 230 


Run no’. B-3 
Flow Rate = 2.6 gpm 


Electrical Power = 80,000 BTU/hr 


Q/A = 96,100 BTU/hr-ft* 
Inlet Pressure = 19.5 psi 

BEe On ps2 
Exit Pressure = 9 psi 


Inlet Temperature = 95°F 


Exit Temperature = 157°F 

ie Two Twi th Tr h Re Nu/pr?*? 
2 197 189 149 169 2440 12900 66.6 
QR 194 186 149 163 2590 #12350 70.7 

185 177 142 159 +2690 12000 eno 

4 ifs 368, agh ato “3266 Jaie® S5—8 
4R 171 163 134 #$$J49 3260 = 811100 Como 
5 te 1368 #126 @22  7elo 9650 200 
5R 1448 140 126 133 7120 9780 182 
6 158 150 i226 138 #4055 # £16150 105 
7 152 Ws 126 135 5340 9930 Sy. 
7R 148 140 126 133 =~ @«=©6720 9780 172 
8 153 145 #118 $132 #3560 9650 90.9 
9 143) 135 dl 123-5 3936 8900 100 
10 1205 lee 103 112 “5065 8050 123 
10R iia elaie = LOG) Lio a OL® 8280 99.5 

Average Re = 10,440 

Average Nu/Pr?: 3 = 90.0 
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Data Sheet 11 


Flow no. B-4 


Flow Rate = 3.0 gpm 
Electrical Power = 79,500 BTU/hr ! 


Q/A = 95,400 BTU/nr-£t* 
imies Pressure = 25 psi 

AR = Gigs eens. 
Exit Pressure = 11.8 psi 


Inlet Temperature = 99°R 


Exit Temperature = 152°F 


fotnie wo twi tp Te h Re -Nu/pr°?? 
2 191 183 145 169 2540 #15300 69,2 
2R 185 lt i145 acim 3620 14350 Silas 
179 171 #+$139 155 2940 #13800 78.4 
4 167 159 132 146 3530 #«412850 93.7 
AR 165 157 132 145 #£=3820 #«%12700 99.8 
5 1446 138 126 132 7630 ~=@«11450 195 
5R 1 6 1386—C:tCiG—(esiad COOSA 232 
6 157 149 126 137 «©4150 #11950 107 
7 150 i14e 126 134 #%§$§5780 #£«£411650 148 
7R My 86 13890CtCi«C«CdGsC‘C( TCO ESO 181 
8 15d eS 486119 «6H 63986) «~“aeGo 101 
9 143° #135 #4112 =+$§4123 4260 #10550 107 
10 1383 125 106 1155 5066 9800 124 
10R i 26 106 1167 4636 9900 116 
Average Re = 12,250 


Average Nu/Pr?:° = Ove 
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Run no. B-5 


Flow Rate = 3.5 gpm 


Data Sheet 12 


Electrical Power = 80,000 BTU/hr 


aoe 


Q/A = 96,100 BTU/nr-ft* 

Inlet Pressure = 32 psi 

Exit Pressure = 14.8 psi — “ae 

Inlet Temperature = 94.5°F 

Exit Temperature = 140.5°F 

Sate wo wi tp Fe h Re Nu/Pr?:? 
2 177 169 135 152 2870 15650 fone 
2R 173 65a) 135m 150m 3280 15400 84.6 

167 159 129 i144 3170 ~# 14700 3320 

4 156 149 123 136 3770 #13750 Oa 
LR 157 l14e> 123) 136) 3760 13750 bee 
5 139 131 #2118 =$dek 7120 #12450 179 
5R W4o 132 $4118 $=#$4125 6720 # £12500 169 
6 150 142 118 130 3920 #13100 99.9 
7 142 134 $118 #126 »5820 #12600 147 
7R 140 132 644118 4125 6630 #12500 167 
8 143 135 112 124 =j;Allo 12450 103 
9 135 127 106 117 #4=«®944580 11650 114 
10 126 118 100 109 5490 #10800 134 
1OR 130 2 100m line Who 169560 1S 

Average Re = 13,220 

Average Nu/Pr?:? = 100.4 
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Data Sheet 13 


Run no, B-6 


Flow Rate = 3.9 gpm 


Electrical Power = 79,400 BTU/hr 


Q/A = 95,300 BTU/hr-ft* 


Inlet Pressure = 38,5 psi 


Average Re = 14,860 


On 


Average Nu/Pr = 111.1 


233 


AP 21 psi 
Exit Pressure = 17,5 psi 
imtet Temperature = 101°F 
Exit Temperature = 141°F 
mane wo Twi th Te h Re —-Nu/Pr 
2 i799 “171 “Sls6 SierS see e700) 17350 T3508 
2R 7] “S6S So tS sooo S6g00 94.7 
168 | 160 131 145 £3310 #16350 S6.6 
4 156 NS lG 37 ee sso etas50 ate 
AR 156 148 126 137 4410 #815350 114 
5 140 132 121 #127 «#8360 814100 212 
5R 40 132 121 126 8910 #13900 223 
6 150 14e 121 132 j%&4540 #14700 116 
7 32 WSs le) les woe 0 200 173 
7R 141 133 121 127 8080 14100 204 
8 115 ls? 9116927 ete o S00 115 
9 139 LS 1 Selly Serer So aes 300 122 
10 129 lel 106 114 #6230 #412500 154 
1OR 133 125 106 116 5020 #12700 124 


jim 
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Run no. B-7 


Data Sheet 14 


Flow Rate = 4.3 gpm 


Electrical Power = 80,000 BIU/hr 


234 


Q/A = 96,200 BTU/hr-ft* 
Inlet Pressure = 44 psi 
AP = 24 psi 

Exit Pressure = 20 psi 

Inlet Temperature = 97 .5°F 

Exit Temperature = 135°R 

nonin ; Tow wil *b Te h Re Nu/ pr?:9 
2 163 155 130 142 #£«:(.3850 #17600 100 
2R 160 152 130 141 4370 17500 114 

159 151 126 138 #3850 #17100 99.4 

h 150 ie 121 131 4510 16100 115 
4R 49 141 #121 #9131 @®«2°44810 ~=# 16100 123 
5 137% 129 116 1227 7580=— 14650 190 
5R 137 129 $116 123 £47400 #815000 185 
6 145 137 #116 #4126 «4740 ~—=615400 120 
if 139 131 116 123 6420 8 14950 161 
7R 137) 1297" 114 ese (5 (Oe ose 190 
8 Wo 8132 46112 «©8122 €64690-~—614850 118 
9 134 126 107 #=\1116 5200 14100 129 
10 126 118 102 =$110 %|}6000 13250 147 

Average Re = 15,780 

Average Nu/Pr?:° = 118,2 
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2 
Data Sheet 15 = 


Run no. Cw] 

Flow Rate = 1.5 gpm 

Electrical Power = 77,500 BTU/hr | 
Q/A = 93,100 BTU/hr-ft* | 
iniet, Pressure = 987 psi 


AP = 5.2 psi 
Exit Pressure = 4,5 psi 
Inlet Temperature = 97 OF 


Exit Temperature = 199°F 


i. Tou Eden Te h Re Nu/Pr°?:° 
2 250 ehe 186 214 #1660 10100 An 8 
eR e5e e4h 186 215 1605 #10100 43.4 

235 227 +%174 +#+.200 1760 £9330 47.5 
4 212 204 161 183 #2160 #8400 58.4 
5 mS SaifO S48 ~SI59uFes0  S7650 114 
5R 170 162 148 #=+§$155 £46650 #6800 179 
6 190 182 148 165 e740 7400 74.0 
7 170 162 148 $155 6650 #6800 179 
7R 172 164 148 156 5800 6900 156 
8 182 174 135 155 2390 #6800 64.6 
8R iS Wi7O 6135 @elS53 B2660 G00 eS 
9 169 161 123 14e 2450 ~~ 866150 66.2 
10 1448 J40 110 125 #3150 5330 Obae 
10R 12 144 110 127 ~ #«©2780 = 8§65430 (520 


Average Re = 7570 
Average Nu/Pr?:? = 63.1 


Peer erg 
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Data Sheet 16 


Run no. C-2 

Flow Rate = 2.0 gpm 

Electrical Power = 78,900 BTU/hr 
Q/A = 94,800 BrU/nr-ft° 


Inlet Pressure = 13.8 psi 


Exit Pressure = 6.5 psi ee 

Inlet Temperature = 98°F 

Exit Temperature = 176°F 

| oe tyo (iwi 4p Mie h Re Nu/pr?9 

2 220 212 166 189 2060 #411600 54.5 
2R 220 22 dG6 @isS Wie2950 daiGoo 54.5 
3 208 200 157 #+%9178 2180 #10800 See 
4 192 184 #147 #421165 #2600 9860 68.9 
5 164. 156 137 #42®147 #4990 8540 130 
5R 158 150 137 144 #«7300 8350 193 
6 M76 @68 237 153 —S3050 8950 80.8 
7 158 150 137 4144 7300 8350 193 
7R 160° 252 a8? G45 Seeec0 8400 167 
8 M69 “16 227 Wael ee 7so 8350 73.9 
8R 165 WibyY 1227 Glo Gei6o 8230 S37 
9 159 151 /118 134 2830 7700 10 
10 we dey 108 S421 Wie 6850 98.3 
10R 16 268 205 G23 Mestico 6970 83.7 

Average Re = 9090 

Average Nu/Pr??? =— ice 
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Data Sheet 17 


Run no, C-3 

Flow Rate = 2.5 gpm | 
Electrical Power = 78,900 BTU/hr 

Q/A = 94,800 BrU/nr-rt* 


Iniet Pressure = 19 psi 


hPa Oep os 
Exit Pressure = 9 psi 
Inlet Temperature = 98°F 
Exit Temperature = 161°F 
ele wo twi Tp Te h Re Nu/pr?*° 
2 197 189 153 171 2630 12800 68.4 
2R 198 190 153 Giga Geasce 12800 CEn6 
188 BOO 1s “86s S27Ie) 1200 (O25 
4 lj 166 wy 152 270) 128/00 85.0 
155 Wer Wo Wee ite 9960 141 
5R 15 143 i1g0- dB6G "o20 9750 183 
6 166 158 130 144 $%§$}3390 #10400 88.1 
a 151 143 130 136 # 7020 9750 183 
7R 152 We a0 Me]? e566 9850 170 
8 NG@m 152 122 pkey cco 9850 82.1 
8R 15g iGo 122 @e6 Bega 9750 one 
9 151 243 124 129 £3270 9200 85.0 
10 7 elem 106) 116 426 8330 107 
HOR 141133 106 120 #& 3510 8500 92.1 


Average Re = 10,480 


Average Nu/Pr??? = 83.6 
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Data Sheet 18 


Fun no, c-4 
Flow Rate = 3.0 gpm | 
Electrical Power = 78,900 BTU/hr | 
Q/A = 94,800 BTU/hr-ft* 
imerets Pressure = 2o5ejecu 

Bee eames 1 
Exit Pressure = 11.8 psi 


Thermo- 


couple Tyo Twi Tb Ts h Re Nu/pr?? 

2 183 175 143 #9159 2960 14100 (aD 
2R 183. 175 143 #159 2960 #14100 25 

173 «61650 SY ear S200 86.2 
4 165 157 130 144 3510 12500 89.5 
5 147 139 «#124 131 46110 #11200 156 
5R 44 #136 124 $4130 ©7580 8 11100 193 
6 156 148 124 #136 3870 #211700 98.6 
7 4 136 124 130 # #97580 #&©11100 193 
7R 1446136 0=«©6124 = =6130)=0—s 7580 ~—S 11100 193 
8 1is26 Te iy IS salon 1si2e6 89.5 
8R 1449 141 #117 #29129 3950 11050 101 
9 W445 137 +4110 124 + +%}3510 #10600 89.5 
10 1s) 25> IOs SE 4oil6 9700 eS 
10R 136 128 $104 $=$116 #3950 9800 Aoi! 


Average Re = 11,800 


Average Nu/Pr?: 3 = Oem 
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Run no. C-5 


Flow Rate = 3.5 gpm 
Electrical Power = 78,900 BTU/hr 


Q/A = 94,800 BTU/hr-ft@ 

Inlet Pressure = 32 psi 

Exit Pressure = 14.8 psi oo 

Inlet Temperature = 98°F 

Exit Temperature 140°R 

eels wo iwi ‘pn te h Re‘ Nu/Pr 
2 Ive lee 137 Asi eezse, 154600 Som 
OR 174 166 137 «+151 3260 15400 # 83.1 

166 158 131 145 £3505 #14700 89.4 

4 157 149 126 137 4030 13850 103 
5 145 137 #4120 129 +5560 #12900 Lhe 
ok We 13 120 Wey eica se7e6 We 
6 152 144 120 132 3950 13250 101 
i 42 8134 ©6120 © «66127 0 8§«€©6©6760)6=—s:« 12900 2 
7R iy ise, i2o ey 72keny dieses 185 
8 Tee Iie i104 “Wey 36ten icc 92.8 
8R 145 137 #124 #125 4120 #12450 105 
9 ise? 134% 109 32 3796) I2eo 96.6 
10 131 123 104 #113 £94850 ~~ 11100 124 
10R 134 =$126 104 #=$+$115 #=;;}4421l0 811350 107 


Average Re = 13,250 


Data Sheet 19 


Average Nu/Pr°?:? = 98,4 
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Run no. C-6 


Flow Rate = 3.5 gpm 
Electrical Power = 137,300 BTU/hr 


Data Sheet 20 


Q/A = 165,000 BIU/nr-ft* 


240 


imiet Pressure = 32 gosu 

Exit Pressure = 14.8 psi e > ailesenh- 

Inlet Temperature = 93°F 

Exit Temperature = 170°F 

same Two Twi Tp Tre h Re Nu/pr?*? 
2 219 205 160 183 3670 19600 96.5 
2R 218 204 160 182 3750 #19500 98.5 

205. 190) 51 eo. 750 109 

4 195 181 2141 161 £=4130 #&#£416700 109 
5 168 154 132 143 £97340 14450 193 
5R 163 149 132 1400 «69430 ~=614150 248 
6 179 165 132 148 $4920 #13750 eo 
7 164. #150 132 141 #48930 14250 235 
7R 163 149 132 lho 9hho 8614150 248 
8 174 160 122 141 4350 14250 115 
8R 170 156 122 139 4860 14000 128 
9 164 150 112 #=+<;9131 £44350 813100 115 
10 148 134 4103 #4119 5330 #11750 140 
10R 15) 237 103 120 4660 11900 128 


Average Re = 15,250 


Average Nu/Pr 


0, 


= 


= 116.7 


p 4 
SF 








ee -*. i. 
eee @4 8 O8CE8 
ee SIP reese 


Perart |i 
eeteeceegese 
stoaeeseaeese 


seetetiat aay 
ese00ee@eeegse 

















Run na, 


Con 


Flow Rate = 3.5 gpm 


Data Sheet el 


Electrical Power = 35,450 BTU/hr 


Q/A = 42,550 BIU/nr-ft* 
Inlet Pressure = 32 psi 
AP = 17.2 psi 

Exit Pressure = 14.8 psi 

Iniet Temperature = 95°F 

Exit Temperature = 115.5°F 

apie wo -wi ‘md te h Re Nu/pr?:? 
2 ioe eee eee 2430 12100 59.5 
2R 12h Jon seeteomen c4cOu 12100 59.5 

130) 12 lO me 1LISeee 26HG. «©6600 64.8 

4 125 i222 10S) 11 Sie ane5o | Oman 
5 121 #2118 #6105 #=+$9111 #243490 #10900 85.4 
5R 109) 1169105) 1105 4ifee Teseo 102 
6 125m 12255 105e 113 2636) inloo 64.3 
7 116% l16QI0O5e) 110m 417em 1L0e0e 102 
7R 118 115 #105 $110 #«=;\}44620 10800 nae 
8 123 120 #103 #112 «2192550 #10900 Geet 
8R 199 116 103 Be 10Gm 3 350m 16700 82.0 
9 178) 115 100 = 107m 2955. 110500 72.4 
10 nee 16 98 104 3580 #10200 Ofa5 
HOR 115 112 98 105 3060 #10300 150 


Average Re = 11,075 


Average Nu/Pr?: 3 = 70.4 
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Run no, c= 


Flow Rate = 4.0 gpm 


Data Sheet 22 


Electrical Power = 78,900 BTU/hr 


eho 


Q/A = 94,800 BTU/hr-ft* 

Inlet Pressure = 39.5 psi 

Exit Pressure = 18 psi ‘=: 

Inlet Temperature = 98°F 

Exit Temperature = 137°F 

ey Mee Twi Os Te h Re Nu/Pr?:? 
2 165e5 N5e.5 132 145 3710 16800 93.5 
2R 165.5 157.5 132 145 3710 16800 93.5 

158 150 127 139 4120 16000 104 

4 151.5 143.5 122.5 133 4500 15300 IIe 
5 140 132 117.5 125 6530 14200 164 
5R 138 130 11755 Wet) 75comn1 Go 190 
6 146 138 117.5 128 4610 14550 116 
7 138 130 117.5 124 7560 14100 190 
7R 137 129 117.5 123 8220 13950 207 
8 Vises 165.5 112.5 iweb Annie 14100 HOSS 
8R 141 133 11285 a2) 4enomn 3950 116 
9 138 130 107 119 4110 13450 103.5 
10 129 121 103 112 5250 12550 132 
10R 132 124 103 Jie 4510NN12800 114 

Average Re = 14,630 

0.3 


Average Nu/Pr = 108.9 
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Data Sheet 23 


Rumeno., C—O } | 
Flow Rate = 4.5 gpm 
Electrical Power 78,900 BTU/hr | 
Q/A = 94,800 BrU/nr-ft~ 
Inlet Pressure = 48.5 psi 

AP = 26.5 psi 
Etro Presgmure = 225 0S1 
Inlet Temperature = 97°F 


Exit Temperature = 132°R 


Seis | Tyo Twi oe. el ea Re _Nr/Pus° 
2 161 153 128 141- 3790 18300 94.7 
2R 160 152 128 140 83950 18150 98.7 
3 154 146 123 135 (4120 #17500 103 
4 147 89139 119 129 ©4730 = 8=16550 108 
5 138 130 114.5 122 6110 #15600 Iss 
5R 135 127 #+.4124.5 121 #7570 #15400 189 
6 44 136 114.5 125 <Ahio 16000 110 
7 16 1265 Wi “el 7 ero wb oo 6 
7R 135 127 #+%\9114.5 121 7560 £15400 189 
8 4100-133 110 121 =+%(4120 ~=©15400 103 
BR 137 129 110 120 4990 15300 125 
9 136 §©128 106 117 4310 14800 108 
10 127. =—-:1119 101 110 5260 13900 131 
WOR 1292512145 101 1112 #£x\44610 14000 115 

Average Re = 15,990 

Average Nu/pr?: 3 = 109.6 
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